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A review  is  provided  of  recent  (December  1976-June  1978) 


developments  in  problems  associated  with  estimating  the  potential 
effects  that  aircraft  emissions  at  cruise  altitudes  may  have  on 
the  earth's  protective  ozone  shield  and/or  surface  climate. 
Background  information  is  provided  as  deemed  necessary  for  con- 
text. The  review  shows  that  one  recently  measured  key  reaction 
rate  (HO,  + NO  + NO^  + HO)  has  had  dramatic  effects  on  the  ozone 
question.  Computed  effects  on  the  ozone  column  of  nitrogen  ox- 
dies  (which,  unless  in  very  large  quantity,  now  cause  an  increase) 
and  water  vapor  emissions  (which,  unless  thermal  feedback  effects 
are  included,  cause  a decrease)  from  supersonic  transports  (at 
17-20  km)  are  now  small  and,  on  balance,  apparently  positive,  at 
ieast  for  moderate  fleets  (several  hundred  aircraft).  Subsonic 
aircraft  also  appear  ,to  cause  small  increases  in  the  ozone  column, 
however,  no  new  modeling  results  are  available.  The  new  HC>2  + NO 
rate  has  affected  model  duplication  of  the  natural  atmosphere 
adversely;  another  important  new  rate  (HO^  + 0^  HO  + 2 O2), 
not  yet  incorporated  in  available  results,  may  reduce  this  dif- 
ficulty. Second-order  effects  (thermal  feedback)  have  become 
important,  particularly  in  modeling  water  effects;  the  modeling 
of  water  transport  processes,  however,  involves  many  uncertain- 
ties. Additional  modeling  studies  are  needed.  Progress 'in 
emissions  measurement  uncertainties  and  brief  comments  on  pos- 
sibly important  climatic  aspects  are  also  included. 
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SUMMARY 


This  report  summarizes  recent  (December  1976-June  1978*) 
developments  in  the  estimation  of  potential  effects  of  aircraft 
emissions  at  cruise  altitude  on  the  protective  ozone  layer  and 
on  surface  climate.  Background  information  is  provided  as 
deemed  necessary  for  context.  The  subjects  covered  include 
developments  in  aircraft  emissions,  ozone  modeling  (atmospheric 
dynamics  and  chemistry),  and  briefly,  climate. 

S.l  BACKGROUND  AND  STATUS  IN  BRIEF 

The  subjects  discussed  herein  have  been  under  continuing 
study  since  about  1970  when  the  advisability  of  continuing 
development  on  the  United  States  supersonic  transport  (SST)  was 
being  intensively  debated.  The  original  environmental  argument 
related  to  exhaust  water  vapor  reduction  of  the  stratospheric 
ozone  layer;  in  1971  a much  larger  effect  on  ozone  was  predi- 
cated due  to  the  nitrogen  oxides  (NO  ) in  the  exhaust.  Climatic 
effects  due  to  water  vapor  and  sulfur  dioxide  (from  sulfur  in 
the  fuel)  were  also  forecast.  Although  the  U.S.  SST  was  can- 
celled, the  issues  remained;  also,  other  SSTs  were  in  develop- 
ment. Accordingly,  to  pursue  the  questions  raised,  a large, 
several-year,  multimillion-dollar  study,  the  Climatic  Impact 
Assessment  Program  (CIAP)  of  the  Department  of  Transportation, 
was  undertaken  in  late  1971.  Results  of  this  study  (Grobecker 
et  ai.,  197^) » as  well  as  of  a concurrent  National  Academy  of 
* 

The  principal  cutoff  date  for  this  work  is  June  1,  1978. 

Where  feasible  and  critical  however,  certain  information 
reported  at  the  World  Meteorological  Organization  (WMO) 
meeting  (Toronto,  June  26-30,  197b)  has  been  included. 
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Sciences  study  (NAS,  1975),  were  released  In  197^-7  , . 
confirming  early  predictions  (but  narrowin  ' the  r : . . 
stratospheric  ozone  depletion  due  to  I’C  , with  wave: 
ozone  ambiguous,  and  also  predicting  climate  chan 
cooling)  due  to  sulfur  dioxide  and  water  emissions; 
adverse  effects  were  detailed  at  L<  n ;th.  Furl  h r 
the  market  situation  for  SSTs  declined  during  )IAP, 
extension  of  the  stratospheric  ozone  depletion  problem  t 
subsonic  aircraft  was  added.  Many  unanswered  ques  S 

In  1975,  to  implement  and  extend  the  findings  cf  , 

High  Altitude  Pollution  Program  of  the  Federal  Aviatic: 
tration  was  established.  At  about  the  same  time,  • - 

recognized  threat  to  the  ozone  layer  from  halocarbcr  s ear. 
intensive  investigation  by  other  groups.  During  the  c • 
these  various  efforts,  new  chemical  reaction  rate  mru/  r.t 
measurements  have  been  undertaken  and  new  chemistry 
included,  all  of  which  have  caused  dramatic  revisi  ns  Li 
puted  effects  of  aircraft  emissions  on  the  ozone  la;. 
comprehensive  review  of  these  various  developments  : r- 

as  of  December  1976  and  distributed  in  1977  (Oliver 
1977)-  The  1977  review  showed  computed  SST  depleal 
on  .zone  (at  constant  N0x  input)  to  have  declined 
(although  uncertainties  increased),  and  that  NO  f . 
would  result  in  an  increase  in  the  ozone  column,  rat;.- 
the  decrease  reported  in  CIAP. 

Changes  since  December  1976  have  continued  • h- 

CIAP  trends  in  NO  effects  on  ozone,  but  have  ir.tr; 

x 

new  difficulties,  in  that  modeling  results  reported  ; 

from  models  which  do  net  adequately  match  the  "n  r ur  '"  atv 

sphere.  Matching  of  the  nature  atmosphere  Li 

sidered  to  be  a necessary  (but  clearly  not  a stiff  !.  ; - • • 

tion  to  credibly  predict  effects  of  perturbation: . 

very  late  chemistry  changes,  which  have  not  yet  t-  ■ r ' 
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should  ameliorate  this  problem.  Available  modeling  results  are 
thus  somewhat  unsatisfactory;  nevertheless,  modeling  efforts  with 
rate  constants  treated  parametrically  permit  certain  observations 
as  follows:  Currently  (see  S.3  below  and  Section  3.5.5  in  the 
train  text),  it  appears  that  the  NO  constituent  of  engine  ex- 
haust  from  both  SSTS  at  17-20  km  (except  in  very  large  numbers) 
ana  from  subsonics  at  9-14  km,  will  contribute  to  an  increase  in 
'he  ozone  column ; water  vapor  emissions,  on  the  other  hand,  at 
least  from  SSTs,  (will)  tend  to  decrease  ozone  unless,  as  seme 
m ■ de ling  results  indicate,  radiative  feedback  effects  compt  nsatt  . 
The  net  effect  on  the  ozone  column,  of  combined  water  and  NO 
■missions , with  current  NO emission  indices , is  apparently 
si  ' ght  y positive  for  both  classes  of  aircraft  in  foreseeable 
numbers . 

Quantitative  estimates  for  the  effects  of  projected  fleets, 
combining  revised  forecasts  with  revised  chemistry,  are  not  yet 
available.  It  seems  quite  probable,  however,  that  the  effects 
over  the  next  20  years  or  so  from  subsonics  will  be  greater  than 
from  supersonies.  Subsonic  aircraft  tend  to  increase  the  ozone 
column  by  increasing  the  ozone  concentration  in  the  upper  tropo- 
sphere. Conceivably,  this  potential  increase  may  become  a mat- 
ter of  concern. 

As  discussed  in  our  previous  report  (Oliver  et  al.,  1977), 
climatic  effects  of  both  subsonics  and  supersonies  are  doubtful, 
but  may  involve  warming;  water  vapor  is  probably  the  emitted 
ingredient  of  greatest  climatic  importance  for  either  type  of 
aircraft.  More  study  is  needed  on  effects  of  altered  vertical 
orefiles  of  radiatively  significant  species  such  as  ozone  and 
water  vapor. 


As  in  the  past,  many  uncertainties  remain.  : enl  ievelop- 
n.ents  are  described  in  the  following  sections,  with  emphasi:  n 
■L-i'  -r:  which  are  felt  to  be  of  current  interest. 


S . 2 AIRCRAFT  EMISSIONS 


5.2.1  NOx  Emissions  Measurements 

Large  uncertainties--factors  of  2 to  6 — still  exist  in  air- 
craft NG„  emission  estimates,  because  spectroscopic  techniques 
have  yielded  far  higher  values  for  NO  than  have  standard  probe 
techniques.  Progress  toward  reducing  these  uncertainties  is 
being  made.  Data  have  been  accumulated  which  show  that  under 
some  conditions,  NC^  in  a hot  gas  stream  can  be  destroyed  on 
sampling,  and  that  measured  NO  contents  vary  with  probe  details. 
No  probe  design,  however,  has  yielded  NO  results  in  hot  gas 
streams  as  high  as  those  found  using  ultraviolet  spectroscopic 
techniques  at  the  Arnold  Engineering  Development  Center  (AEDC). 
The  spectroscopic  techniques  are  difficult  and  are  currently 
being  investigated  in  detail.  Final  resolution  of  the  differing 
results  obtained  by  probe  ana  spectroscopic  techniques  may  well 
show  that  both  techniques  have  beer,  in  significant  error  under 
some  circumstances.  A better  understanding  of  proper  measure- 
ment techniques  should  result. 

5 . 2 . 2 Fleet  Projections 

A recently  published  (February  1977)  FAA  document,  Aviation 
Futures  to  the  Year  2000,  which  is  limited  "almost  exclusively 
to  the  United  States,"  shows  minimal  expectation  of  SST  traffic 
increase.  Of  five  scenarios  considered,  only  one,  "Expansive 
Growth,"  includes  any  SST  traffic;  this  implicitly  unlikely 
scenario  includes  some  230  SSTs  by  the  year  2000,  with  about 
100  by  1990.  An  aircraft  manufacturer's  projection  (McDonnell- 
Douglas)  includes  118  SSTs  by  1990.  Continued  growth  in  the 
subsonic  fleet  is  generally  expected. 

5 . 2 . 3 Fleet  Emissions 

A detailed  global  emissions  estimate  provided  by  British 
Aircraft  Corporation  is  included  for  a small  Concorde  fleet 
(nine  aircraft).  Ozone  change  calculations  for  this  fleet  using 
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current  chemistry  have  not  been  attempted,  but  would  be  ex- 
ceedingly small. 

S.2.4  N0x  Emission  Reduction 

Some  progress  in  achieving  lower  NO  emission  indices  in 

X 

advanced  combustors  has  been  reported,  but  the  problems  in 
achieving  substantial  reductions  in  practical  systems  are  evi- 
dent. A significant  aspect  is  that  technological  advances 
leading  to  improved  performance  (lower  specific  fuel  consumption) 
appear  to  be  in  conflict  with  those  leading  to  lower  emissions; 
if  conventional  or  modestly  modified  combustor  technology  is 
used,  i.e.,  without  emissions  reduction  techniques,  a substan- 
tial increase  in  NO.^  emission  index  with  time  could  be  expected. 


S . 3 OZONE  MODELING  DEVELOPMENTS 
S . 3 . 1 1-D  Transport-Genera) 

N^O  concentration  vs.  altitude  measurements  at  various 
seasons  and  latitudes  by  Schmeltekopf  et  al . (1977)  have  resulted 
in  new  estimates  of  average  vertical  transport  as  used  in  1-D 
models.  Results  show  average  vertical  transport  in  the  middle 
stratosphere  (17-27  km)  to  be  approximately  twice  as  rapid  as 
derived  earlier  by  some  modelers  from  CH^  measurements.  The 
powerful  effect  of  including  one  tropical  profile  is  evident. 

The  data,  however,  provide  no  information  for  the  region 

near  the  t.ropopause  critical  to  computation  of  aircraft  effects. 
This  region  involves  widely  differing  estimates  by  various 
modelers.  A 1-D  analysis,  (Appendix  B),  of  the  transport  of 
Zr-95  from  a source  at  40°N  across  the  tropopause  indicates 
stratosphere-to-troposphere  transport  rates  to  vary  seasonally 
but,  on  the  average,  to  be  more  rapid  than  implied  by  some  K„ 
profiles  (e.g.,  Kunten,  1975;  Turco  and  Whitten,  1977)  and 
slightly  slower  than  In  others  (e.g.,  Chang/1976);  it  is  recog- 
nized, however,  that  the  analysis  is  preliminary  in  view  of 
certain  simplifications  introduced.  A brief  examination  of 


available  3-D  modeling  results  from  NOAA-GFDL#  indicates 
(qualitatively,  because  the  computations  are  not  fully  reduced) 
support  for  1-D  models  (such  as  Cnang/1976)  which  show  rather 
gradual  transition  across  the  "globally-averaged"  tropopause; 
the  3-D  modeling  work  also  has  shown  earlier  that  feedback 
( troposphere-to-stratosphere ) of  excess  C-19 is  significant 
after  one  year  or  so,  suggesting  that  this  tracer  is  applicable 
primarily  at  relatively  early  times  following  injection.  It 
can,  of  course,  be  debated  whether  tropospheric-stratospheric 
transport  can  be  adequately  represented  in  such  models  for 
sources  near  the  tropopause. 

5.3.2  1-D  Transport  of  Water  Vapor 

Injections  of  water  vapor  into  the  stratosphere  by  aircraft 
have  become  of  increasing  interest  in  terms  of  effects  on  ozone 
and  climate,  as  noted  above.  Appropriate  procedures  for  model- 
ing this  tracer  are,  however,  in  some  doubt:  sources  and  sinks 
are  not  well  established.  The  tropical  tropopause  cold  trap 
concept  does  not  operate  satisfactorily  in  1-  or  2-D  models. 

More  study  on  this  problem  is  needed,  particularly  if  thermal 
(radiative)  feedback  effects,  which  are  evidently  important, 
or  both  chemistry  effects  and  water  content  itself  (through 
changes  in  the  tropical  tropopause  temperature)  are  to  be  esti- 
mated. Only  3-D  modeling,  perhaps  using  models  more  advanced 
than  those  presently  available,  would  seem  to  be  adequate  for 
the  feedback  estimates;  as  a practical  necessity  however,  tech- 
niques to  decouple  this  aspect  from  the  detailed  chemical  cal- 
culations must  be  considered. 

5.3.3  Trace  Species  Measurements 

Of  the  many  trace  species  measurement  programs,  two  merit 
particular  mention  because  of  the  questions  raised  by  the  results 
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obtained.  The  first  of  these  is  the  various  Cl  and  CIO  mea- 
surements of  Anderson  et  al.  (1977);  one  set  of  these,  in 
particular,  shows  an  inexplicably  large  stratospheric  CIO  con- 
tent (~  8 ppb ) . Chlorine  content  and  chlorine  chemistry,  both 
of  which  (along  with  the  measurements  themselves)  are  placed  in 
question  by  these  measurements,  affect  the  response  of  strato- 
spheric ozone  to  NO  injections.  The  second  relates  to  the  N02 
measurements  by  Noxon  (1978),  who  finds  NO-  contents  in  the 
troposphere  to  be  an  order-of-magnitude  smaller  than  previously 
believed.  N0o  content  in  the  troposphere  affects  tropospheric 
chemistry  (OH,  0-,  contents)  and  thereby  the  response  of  tropo- 
spheric (and,  in  fact,  stratospheric)  ozone  to  NC^  and  numerous 
other  injectants. 

S.3.4  Chemistry 

The  single  most  significant  recent  development  in  ozone 
chemistry  and  its  modeling  has  been  the  direct  measurement  of 
the  rate  of  the  reaction  of  H0o  with  NO  (Howard  and  Evenson, 
1977),  indicating  a 40-foid  increase  in  rate  at  room  temperature 
over  a previous  (CIAP-era)  evaluation.  The  temperature  coef- 
ficient has  been  reported  to  be  moderately  positive,  the  rate 
increasing  slightly  with  decreased  temperature  (Howard,  19)7")- 
Use  of  this  rapid  rate  leads  to  substantially  increased  ozone 
content  in  model  "natural"  atmospheres,  an  increase  which  is 
particularly  evident  in  2-D  models  and  is  of  sufficient  magni- 
tude to  be  of  concern.  However,  direct  measurement  of  the 
h’O^  + 0^  reaction  rate  has  also  been  reported  very  recently 
(Zahniser  and  Howard,  1978);  the  new  rate  should  ameliorate 
the  excess  ozone  problem,  but  evaluation  runs  are  not  yet  avail- 
able. K02  chemistry  is  apparently  becoming  moderately  well 
established,  although  important  uncertainties  still  exist  (as, 
e.g.,  in  the  HO.,  + OH  and  methane  oxidation  reactions)  and 
possible  pressure  dependencies,  temperature  effects,  etc., 
still  need  investigation. 
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A number  of  species  have  been  noted  as  being  of  new  or 
altered  significance;  among  these  are  HCuNO., , H^O^,  and  H0C1. 
Certain  reactions  and  trace  species  contents  in  the  troposphere 
are  also  of  indirect  but  potentially  great  significance  to 
stratospheric  ozone;  reactions  affecting  troposphere  OH  content 
affect  the  lifetime  of  methane  and  hydrogen-containing  chloro- 
carbons,  and  thus  the  flux  of  such  materials  into  the  strato- 
sphere . 

S.3.5  SST  Effects  on  Ozone 

Three  recent  1-D  studies  are  reviewed:  as  indicated  above, 
none  of  these  specifically  include  the  latest  HO^  + 0,  reaction 
rate,  but  available  parametric  evaluations  indicate  clearly 
that  the  new  rate  will  increase  the  probability  that  N0x  in- 
jections at  SST  altitudes  will  lead  to  an  increase  in  the  ozone 
column.  In  general,  as  computed,  the  three  studies  reviewed 
agree  that  recent  chemistry  changes  (primarily  the  H0„  + NO 
change)  lead  to  local  enhancement  of  ozone  from  NO  injected 

X. 

at  SST  altitudes  (17-20  km).  The  net  column  effect  for  "Moder- 
ate" injections  (of  the  order  of,  say,  the  expansive  growth 
230  SST,  year  2000,  fleet  noted  above,  is  probably  positive  as 
indicated,  but  depends  nevertheless  on  the  exact  choice  of  re- 
action rates,  emission  indices  (NO  /Ho0  ratio),  and  a number  of 
modeling  details  (K^  profile,  feedback  effects).  For  "very 
large"  injections,  however  (the  magnitude  depending  on  the  K 
profile  used),  column  ozone  depletion  is  expected  due  to  "satura- 
tion" of  the  ozone-enhancement  effect  in  mid-altitudes  and  an 
increase  in  destruction  at  higher  altitudes.  Water  vapor  addi- 
tions cause  slight  ozone  decreases  in  the  two  models  where  its 
effects  have  been  studied  under  fixed  temperature  assumptions; 
if  thermal  feedback  (stratospheric  cooling  effects  of  the  added 
water)  effects  are  included,  the  effect  of  water  additions  on 
‘he  ozone  column  become  either  negligible  or  enhancing.  As  in- 
dicated above,  higher  order  modeling  of  such  feedback  effects  is 
desirable,  as,  obviously,  would  be  better  knowledge  of  thermal  an 


possibly  pressure  effects  on  rate  coefficients.  The  water/!: 

X 

ratio  is  clearly  a variable  of  interest;  there  may  be,  depei.ii 
on  further  evaluation  of  feedback  effects,  an  optimum  cruise 
value  leading  to  minimum  risk  to  the  ozone  column,  probably 
varying  with  fleet  size  and  altitude. 

5 . 3 . 6 Subsonic  E f f ec ts  on  Ozone 

No  new  results  are  available  for  the  effects  of  subsonics 
on  the  ozone  column.  The  new  rate  for  HO,  + NO,  however,  can 
be  expected  to  increase  ozone  enhancements  fcr  specified  NO 
inj ections . 

5 . 3 . 7 Model  V a 1 i d a t i o n 

As  noted  above,  use  of  "current"  chemistry  (i.e.,  wi 
the  late  HO^  + 0^  rate  result)  to  model  the  "natural"  at 
involves  some  new  questions,  as  there  is  a tendency  fcr  the 
models  to  generate  too  much  ozone;  the  problem  is  :v  • evident 
with  2-D  models.  The  newly  reported  HO.  + 0-,  rate  will  allev'- 

<-  J 

ate  this  problem,  but  by  how  much  is  not  yet  established.  For 
the  perturbed  atmosphere,  some  partial  success  in  matching 
short-term,  observed  ozone  depletions  to  model 

been  achieved  in  both  1-D  and  2-D  models  for  solar  \ r ton  < vents, 
but  unexplained  aspects  exist.  Moreover,  significant  apparent 
anomalies  still  exist  over  the  longer  term  in  comparing  , i 
destruction  predictions  to  ozone  records  in  the  csso:;  ;f  nucl 

weapons  tests  and  chlorofluorocarbon  releases.  A i t : i 

tor  appears  to  be  variation  in  solar  ultravi  let  output  with 
time.  The  various  modeling  efforts  all  ne<  11  be  r<  xai  Ln<  3 
in  the  light  of  up-to-date  chemistry.  Studies  witl  - r 3-1 
models,  recognizing  the  array  of  known  ongoing  perturbinj 
ences,  will  be  of  particular  Interest  to  the  questions  involved. 

A time-dependent  2-D  "catalog"  of  perturbing  influences  vei  t:. 
period  1955-1975  has  been  prepared  (Bauer,  1973). 


■HP 
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S . 4 CLIMATIC  EFFECTS 


Little  new  information  has  been  developed  specific  to  the 
climatic  af t f i rcraft  effluents.  Climatic  research,  of 

course,  continue.-:,  and  all  such  work  has  some  relevance.  A 
brief  review  '■  included  of  current  problem  areas  and  probable 

•ff--  • ' . It  is  noted  that  the  probable  climatic  effect 

- m ■ ■ ;-live<  nthr  p genic  pollutants  (C0_,  N?0,  halo- 

f rreni  concern  is  to  cause  warming;  these  effecl  s 
wi'l  add  co  any  effects  (also  probably,  a small  warming)  due  to 
aircraft  effluent o. 

S . 5 RECOMMENDATIONS 

5 . 5 . 1 Ozone  Model  Development 

axis ring  models  need  examination  to  determine  the  cause 
a n d ‘ : re:.,  v , a tendency  for  such  models  to  yield  too  much 
n<  Ln  natural  tmosphere;  the  new  H0„  + 0,  rate  constanl 
("ahniser  and  Howard,  117?)  should  help. 

5 . 5 . 2 Fleet  Ef f ec  ts 


lew  i Lng  runs  j re  needed  to  develop  est Lmates  vs.  1 
of  the  eff  n the  ozone  column  of  future  (particularly 

subsonic)  ts . Th<  s<  calculati  >ns  should  be  time-dependent 
and  shot  i ii  lude  expected  time-dependent  variations  in.  the 
atr.osph  -r!  • hal  ••arbon  concentrations.  Consideration  should 
be  giver;  to  whether  (and  when)  the  expected  ozone  column 
effect.'  dv  ' . ur. sonic  aircraft  may  become  a matter  of  concern. 

• . : b Ln  both  1-D  and  2-D,  and  should  be 
' : a*  .i.iv  when  (and  if)  modeling  pecularities  introduced 

by  [•  • :n  ■ Ln  chemistry  are  clarified. 


S.5.3  Ozone  Model  Validation 


c 


iumn  in  recent  iecade 


factor  causing  changes  in  the  ozone 
appears  to  be  solar  ultraviolet 
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v ar iab i L ty , but  the  variations  r < j 3 1 al  1 juanti- 

catively . It  is  essential  for  model  valid- '.ion  t . ot hair  dr.- 
on  solar  ultraviolet  variability  with  time  if  It  r.  • - rr  n:r  :.ro- 
p genic  influences  are  to  be  verified. 

5.5.4.  Chemistry 

Continued  study  of  H0-,  chemistry  is  called  f r to  remove 
doubts  about  pressure  effects  and  competing  r ■.  c:  enisr.u ; ! r. 
some  cases  (HO,,  + OH,  e.g.),  additional  rate  >r  emperature 
sensitivity  data  are  needed.  Some  spe  ies  , . . , 

KCHO ) , and  their  formation  and  breakdown  path:  an.  pr  ducts 
■ 1 mori  Li  n stigation.  Some  reactions,  considers  ui  Lmj  *tan 
in  initial  screening  (such  as  OH  + 0C1)  merit  ex- n ir 
view  of  ar.  evident  (if  perhaps  unlikely)  r t e:  t ’ • ' ' •' '.  • ? • 

Methane  oxidation  chemistry  1 sc  needs  ii  ir  ased  tt<  Lon. 

5 . 5 . 5 Transport  and  Feedback  Effects 

Wide  discrepancies  still  exist  among  tin  v ia  - - i-l 

on  the  rate  of  transport  across  the  tropopausc.  These  aiscre- 
p ancles  need  resolution  through  examination  of  the  behavior 
f rious  tracers,  particularly  from  thos<  lei  sited  inti 
stratosphere  at  known  altitudes  and  latitudes,  fro:,  studies 
- ant ial  vorticity , and  from  three-dimensi  ■ • ff  rl 

Water  vapor  transport  is  of  particular  interest;  .audios  f H7C 
from  nuclear  weapons  tests  may  be  of  value.  Thermal  feedback 
effects,  as  on  stratospheric  arid  tropical  trep^pause  tern;  era- 
tures,  need  investigation,  alm<  si  ainly  r<  iuii  '.-  us 

. .■  i : y ' r v 1 ; 3-D  me  le  Ls  . Pol  ■.  T phenomenological 
studies,  parameterized  feedback  procedures  in  lov.f  diner.c;  u- 
ality  models  may  be  feasible  for  stud’es  of  the  r . ne  cheni  c-try 
involved . 

5.5.6  N0x  Emission  Index  Uncertainties 

Possible  errors  in  probe  techniques  hav-  I • - • in 

w rk.  Errors  may  also  exist  in  th<  s.  :tros  • tech- 
niques. Work  underway  should  be  continued. 

S-ll 
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S . 5 . 7 Climate 


thermal  feedback  effects  have  been  shown  to  be  significant 
- 1 ■ ■ ( better  quantify  these  complex  effects , 

however’  and  particularly  to  estimate  changes  in  stratospheric 
w'  -r  content  resulting  from  changes  in  tropical  tropopause 
' ; ‘re,  iditional,  pi  rat  ly  3-1  , rr  delir  st  idle  ar< 
r.ee  je:i . Studies  in  1-D  models  of  the  coupled  climatic-chemical 
: neou;  injection;  f NO,  , H?0,  SO  . .•  . -■ 

: ' s > w ' th  r<  sultant  char  g(  . in  he  alti- 
■'  1 *'  ne,  md  t her  trace  . : eci<  . . nd, 

however,  be  f interest. 


1.  INTRODUCTION  AND  BRIEF  REVIEW 

This  report  presents  continuing  developments  in  the  problem 
of  predicting  effects  that  future  fleets  of  high-flying  aircraft 
may  have  on  the  ozone  column  and  on  surface  climate.  The  issues 
arose  in  the  1970-71  time  frame  during  debates  on  the  advisability 
of  continuing  the  United  States  supersonic  transport  (SST)  pro- 
gram. Initial  concern  was  with  water  effluents,  with  a postu- 
lated decrease  in  ozone  resulting  and  with  some  possible  impact 
on  climate.  The  concerns  shifted  laoer  (1971)  as  rate  constants 
were  reevaluated,  to  emphasize  NO  which,  in  a catalytic  ozone 
destruction  cycle  (Crutzen,  1970;  Johnston,  1971),  was  postulated 
to  lead  to  substantial  ozone  destruction  [3  to  50  percent, 

Johnston  (1971),  depending  on  assumptions],  should  large  fleets 
of  SSTs  come  into  operation  (Crutzen,  1971,  1972;  Johnston, 

1971).  Concern  with  the  climatic  effects  of  SC,,  arising  from 
sulfur  in  the  fuel  and  the  sulfuric  acid  derived  therefrom  were 
also  expressed  in  this  time  frame. 

Although  the  United  States  SST  program  was  canceled,  the 
issues  and  the  scientific  problems  remained,  as  other  SSTs  were 
in  development  elsewhere  and  the  possibility  remained  that  new 
SSTs  might  be  developed  in  the  future.  To  pursue  these  matters, 
a large  multimillion  dollar,  multidisciplinary  program,  the 
Climatic  Impact  Assessment  Program  (CIAP),  was  sponsored  by  the 
Department  of  Transportation,  beginning  in  1971;  the  Report  of 
Findings  from  this  program  was  released  in  early  1975  (Grobecker 
et  al.,  197A).  This  study,  as  well  as  concurrent  studies  by 
the  National  Academy  of  Sciences  (NAS,  1975)  and  by  the  United 
Kingdom  (COMESA,  1975),  largely  affirmed  that  IJO  emissions  from 


aircraft  would  have  a deleterious  effect  on  the  ozone  column,  al- 
though the  estimated  magnitudes  of  the  effect  differed.  Compari- 
son to  pre-CIAP  estimates  was  difficult  due  to  the  ’wide  ranges 
ranges  quoted,  the  differing  averages  involved  (global,  hemi- 
spheric, or  corridor),  differing  estimates  concerning  the  average 

flight  hours  per  dav  r.er  SST  and  thus  NO  emissions  to  be  ex- 

x 

pected  per’  year,  etc.  A nominal  Northern  Hemisphere  figure  for 
- rge  . ivanced  mach  2,7  SSTs  operating  at  19.5  km  (not  500  cf 
the  canceled  U.S.  SSTs)  according  to  the  CIAP  Report  of  Findings 
would  be  9 percent  depletion  in  the  Northern  Hemisphere;  the  NAS 
prc  lure  whic)  . • ft  flow  estimates,  gave  about 

15  percent.  Effects  per  unit  of  (continuous)  NO  emission  were 
noted  to  decrease  with  decreasing  altitude;  however,  even  sub- 
sonics,  and  particularly  advanced  subsonics  operating  to  nearly 
km,  were  implied  to  represent  a significant  threat  to  the 
ne  column  in  view  of  expected  growth  rates.  Water  emissions 
m ■ nsidered  tc  have  little  effect  on  ozone.  Th<  JIA.E  report 
gave  considerable  discussion  to  the  cooling  climatic  effect  cf 
. w h w is  reportedly  larger  (at  0.05  percent  sulfur  in  the 
fuel),  than  the  warning  effect  of  water  vapor  and  the  small 
her  exhaust  ingredients.  Climatic  effects  were 
.•.til  i Implied  to  be  rather  small,  and  of  more  remote  concern  in 
i.  time  sense,  than  were  the  ozone  effects. 

:r.  1975  j in  order  to  implement  the  findings  of  and  to  extend 
the  investigations  carried  out  in  CIAP,  the  High  Altitude  Pollu- 
tion Program  (HAPP)  of  the  Federal  Aviation  Administration  was 

\bou1  tl  is  time  also,  the  effects  of  halocarbons  (chloro- 
fi uzrome thanes  and  others)  on  the  ozone  layer  became  a matter 

i ;■  i . i at5  Dual  scientific  interest,  and  much  research 
was  undertaker:.  This  research  made  evident  the  greater  complex- 
ity of  stratosr  'neric  chemistry  than  had  been  recognized  in  CIAP; 
m re  importance  to  the  aircraft  problem,  however,  certain 
1 n rates,  previously  evaluated  somewhat  Indirectly,  were 
begun  z : be  measured  directly. 
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All  of  these  matters  were  reviewed  in  a previous  report, 
prepared  under  FAA/  HAPP  sponsorship  (Oliver  et  al.,  1977),  based 
on  information  available  as  of  December  1976.  The  report  empha- 
sized developments  in  the  1975-76  time  period  and  noted  changes 
which  had  occurred  relative  to  findings  and  assumptions  reported 
in  CIAP  (Grobecker  et  al.,  197*0  and  in  NAS,  1975-  A brief 
review  of  that  work,  conducted  in  the  1975-76  period  and  of  the 
developments  which  resulted,  follows. 

1.1  FLEET  PROJECTIONS  AND  ALTITUDE  DISTRIBUTION  OF  EMISSIONS 

Changes  in  fuel  prices  and  general  economic  conditions 
between  the  initiation  times  of  CIA?  ana  of  HAPP,  and  recognition 
of  the  need  for  better  altitude  resolution  for  emissions  near 
the  tropopause,  made  new  projections  necessary.  These  were  made 
and  resulted  in  substantially  altered  and  more  detailed  projec- 
tions. These  projections  did  not  include  a large,  mach  2.7 
advanced  SST  by  1990,  as  had  CIAP.  In  addition,  future  fuel 
consumption  by  SSTs , expected  fleet  growth  rates,  and  fuel  con- 
sumption by  advanced  subsonics  at  peak  cruise  altitudes  (near 
1A  km)  were  all  found  to  be  much  lower  than  in  CIAP.  However, 
a caveat  persisted,  in  that  NO  emission  estimates  were  in  rues- 

X 

tion  because  of  a large  measurement  discrepancy  between  accepted 
probe-sampling  techniques  and  a new  in  situ  ultraviolet  spec- 
troscopic technique.  The  spectroscopic  technique  indicated  much 
large  values  for  N0x. 

1.2  N0x  EFFECTS  ON  THE  OZONE  COLUMN 

A review  of  the  chemistry  (reaction  rates)  controlling  the 
production  and  destruction  of  ozone  in  an  N0x-HCx-0  atmosphere 
indicated  that  uncertainties  were  much  larger  than  recognized 
in  CIAP;  it  was  shown  that  N0x  emissions  from  SSTs  might  actually 
increase  the  ozone  column  slightly  rather  than  decrease  it,  or 
the  effects  could  be  substantially  more  deleterious  than  reported 
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in  CIAP,  all  of  these  potential  results  being  within  reaction 
rate  uncertainties.  The  uncertainties  were  largely  in  reactions 
involving  the  HO,,  radical.  Uncertainties  in  parameterized  trans 
port  also  affected  the  results.  Overall,  for  a specified  NO 
injection  rate  corresponding  to  a very  large  fleet  (2000  molecui 
cm^-sec  over  1 km  at  17  km,  globally  equivalent  to  about  1,00 G 
Concordes),  results  varied  from  -17.0  percent  to  + 2.2  percent 
global  ozone  change. 

Water  effects  from  aircraft  exhaust  were  also  studied 
briefly  in  these  uncertainty  studies.  It  was  found  that  when 
NO  negative  effects  on  ozone  are  maximized  by  a choice  of  re- 
action  rates,  water  vapor,  which  is  also  present  in  the  exhaust, 
increases  ozone;  where  MO  effects  on  ozone  are  minimized  or 
positive,  the  water  emissions  reduce  ozone. 

During  the  course  of  this  and  related  work,  recommended 
reaction  rates  were  being  changed  and  interactions  of  HGy  and 
NO  with  CIO  chemistry  were  being  recognized.  One  particular 
reaction  rate  (for  HC , + HO  -*■  Ho0  + 0^)  was  revised  downward 
(somewhat  excessively  by  current  estimates)  ten-fold;  this  de- 
creased computed  C5T  effects  on  ozone  some  10  percent,  but  in- 
creased effects  of  halocarbon  injections.  Interactions  with  CIO 
cycles,  which  yield  reduced  (cr  positive)  SST  NO  effects  on 
ozone  with  increased  chlorine  content,  and  other  changes  also 
combined  to  show  reduced  deleterious  effects  of  SST  NOx  on  ozone 
In  addition,  ozone-forming  reactions  involving  methane  oxidation 
and  nitrogen  oxides,  recognized  but  not  included  in  CIAP,  were 
put  into  some  models.  Inclusion  of  methane  oxidation  reactions 
in  2-D  model  studies  showed  that  subsonics,  in  fact,  appear  to 
increase  rather  than  decrease  the  ozone  column;  in  addition, 
while  supersonics  tended,  in  these  calculations,  to  decrease  the 
ozone  column,  the  combined  effect  of  the  projected  1990  "high" 
traffic  fleet,  including  both  supersonics  and  subsonics,  was  to 
increase  ozone  slightly.  The  increase  varied  with  season  and 


latitude,  with  a maximum  enhancement  of  ab  ut  1.4  r-reent.  Un- 
certainties were  recognized  to  be  large  ir.  1.  th'-se  cal  .'Ulatl  r.s . 

1.3  CLIMATIC  EFFECTS 

Little  new  climate  modeling  work  war  ur.dert ekor.  ir.  the  HA:  F 
program.  However,  as  part  of  the  review  ir,  . ' • - : ■ .,  1 97", 

the  postulated  effects  of  the  various  exhau  * -dients  were 

examined  in  detail,  using  available  models,  estimates  of  resi- 
dence times , etc . The  work  was  aimed  ‘ determining,  ' 
degree  possible,  the  effect  of  SST  exhaust  or  climate,  ra- her 
than  focusing  on  the  effects  of  SO  , which  results  from  re- 
movable (at  some  cost)  impurity  ir.  the  fuel.  These  cal.-ul.rl  s 
indicated  that,  according  to  available  me  1<  Is,  h<  fe  :t 

(warming)  is  greater  than  the  sulfur  effect  (al  .05  I 


sulfur) , 

or  other  effects 

( N , 

ozone  change.: ) ; 

the  net  resul 

appeared 

to  be  a warming 

due  to 

SST  aircraft  ex 

Vq  *"•  ’ ' +■  v*  *n  f*  Vo  & V» 

than  a c 

ooling.  The  same 

c a 1 c u 

lat  1 ns  ..  r n<  t 

applied  tc  th 

subsonic 

fleet,  it  being 

argued 

that  the  models 

, which  dc  not 

include 

contrail  effects 

(which 

may  \ mi nan 1 

a re  Li  ippli- 

cable.  Climatic  effects  for  the  SST  port 

HAFP  fleet  were  found  to  be  small  (Q.‘M  1 . 0.  h ; i.  :• , : 

large  fleets  of  advanced  SSTs , climatic  effects  rr  v : 

might  well  limit  fleet  size,  . sumin  me  specif!  rj 

value  (0.1  K,  for  example)  fron  Ircr  fi  ild  i be  ex  eded. 


It  is  evident  from  the  above  that  many  ciu--  res  h • he  ier- 
standing  of  aircraft  effects  occurred  in  the  per:  r:  r • •*  w n 
December  1974  (the  date  of  the  CIAF  Report  T Fi:  lings;  Irsbecker 
et  al. , 1974)  and  December  1976,  the  off  lat  • th<  Iver 
et  al.,  1977,  report.  The  rate  of  change  was  sucl  that  furtl 
shanges  ;ould  b<  inticipated.  rhes<  indeed  icurred , r- 

ticularly  in  the  ozone  modeling  are  . Relatively  . ' • • . ri  ent 
attention  has  been  given  t o olimatic  ff<  •'  f rcraf  ithaust, 
and  equivalently  dramatic  changes  are  lacking.  This  rep  r:  sum- 
marizes these  recent  developments,  as  of  Turr  1,  197'  ; wh-  r«- 


feasible  and  critical  however,  certain  inf< rmatj  i presented 
the  World  Meteorological  Organization  in  Toronto  (June  2(1-30, 
1978)  has  been  included.  No  attempt  to  include  other  develoj - 
ments  which  may  have  been  reported  during  June  1978  has  been 
made.  An  attempt  is  made,  however,  t giv<  enougl  addi 
background  material  so  that  coherer.ce  and  perspective  are  r.ei:  - 
tained . 

ihe  material  that  follows  is  organized  in  the  same  general 
fashion  as  have  been  previous  reports  (Grobecker  et  al . , 19~-; 

- liver  et  al.,  1977);  thus,  aircraft  emissions  are  first  Jis- 
cussed,  followed  by  a treatment  of  ozone  effects  and  (minimal iy) , 
c 1 imatic  effects . Conclus ions  and  rec  ommendati  ns  fron  thes< 
efforts  have  been  brought  forward  to  the  summary  section,  al rear- 
presented  . 
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2.  AIRCRAFT  EMISSIONS 


2.1  NO  MEASUREMENT  UNCERTAINTIES 

2.1.1  Introduction 

As  discussed  in  the  earlier  report  (Oliver  et  al . , 1977), 
optical  measurements  of  NO  in  aircraft  exhaust  have  disagreed 
sharply  with  measurements  made  according  to  standard  probe  tech- 
niques (SAE  ARP  1256).  The  disagreements  were  first  noticed  in 
work  begun  at  AEDC  in  1972  (Davidson  and  Domal,  1973)  and  are 
not  yet  resolved,  but  some  new  information  is  available.  A 
brief  history  of  the  problem,  some  reinterpretation  of  old  data, 
and  these  recent  developments  follow. 

2.1.2  Background  Data 

Representative  data  obtained  by  Davidson  and  Domal  (1973) 
are  given  in  Table  2-1.  The  tests  were  made  with  a J93  engine 
mounted  in  an  altitude  test  cell.  The  optical  procedure  (see 
McGregor  et  al.,  1973)  utilized  the  NO  (0,0)  y-band,  measuring 
absorption  in  individual  spectral  lines  in  the  ultraviolet  (225 
to  227  nm) . The  probe  procedure  used  standard  techniques,  with 
the  gas  samples  being  brought  through  an  orifice  and  a sampling 
line  to  analytical  apparatus. 

Note  that  large  disagreements  occurred,  particularly  at 
higher  mach  numbers  and  altitudes  and  under  afterburning  con- 
ditions. No  data  were  obtained  at  subsonic  conditions. 


TABLE  2-1.  AVERAGE  NO  CONCENTRATION  IN  THE  EXHAUST  OF  A YJ93-GE-3  ENGINE 
COMPARISON  OF  VALUES  MEASURED  IN-SITU  BY  UV  ABSORPTION  AND 
MEASURED  BY  SAMPLING.  (DAVIDSON  AND  DOMAL,  1973) 


Simulated  Flight 

Condi tions 

Concentration 

Measurements,  ppm 

Mach 

Al ti tude, 
thousands  of  ft 

Power  Setting 

Absorption 

Sampl ing 

1.4 

35 

Mil itary 

75 

56 

2.0 

55 

Mi  1 itary 

165 

70 

Minimum  Afterburning 

175 

80 

Maximum  Afterburning 

278 

110 

2.6 

66 

Mi  1 i tary 

323 

100 

Maximum  Afterburning 

617 

130 

It  is  of  interest  that  analytical  data  obtained  on  sample 
gases  which  had  passed  through  a probe  and  sample  line  were  in 
reasonable  agreement  using  several  techniques,  including  chemi- 
luminescence, dispersive  infrared,  and  wet  chemistry  techniques 
(and,  as  will  be  shown  shortly,  by  UV  spectrometry).  Dispersive 
infrared  data  (on  the  probe-sampled  gases)  wer-  reported  by 
Forney  (1973)  and  were  within  10  or  15  percent  of  the  chemilumi- 
nescence techniques  under  all  conditions.  Comparison  of  the  wet 
chemistry  data  to  that  based  on  the  chemiluminescence  data  were 
not  made  by  Davidson  and  Domal,  but  are  reported  in  Table  2.2 
using  their  data. 

These  large  disagreements  prompted  further  work.  In  one 
study,  at  Wright -!  ittersoi  AFB,  T-56  jet  engine  combustor  was 
used,  and  two  different  t ! :al  in-situ  techniques  (IF.  and  UV) 
were  used  along  with  conventional  sampling  procedures.  Two 
groups  of  Investigai  rs  w<  re  inv<  lv<  ne  gi  ui  fr  m AEI  ' (Few 
et  al . , 1976)  and  one  from  Aeronutronic  (Gryvnak  and  Burch,  1976) . 
The  inlet  air  temperature  to  the  combustor  was  756  K (900°F) , 
pressure  was  120  psig,*  and  a fuel-' • -air  rati  of  0.013  was  used. 

W.  Roquemore  , Wright -Tat ter.: on  AFB,  private  communication, 

12  April  1977. 
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It  is  assumed  that  analyses  are  on  a dry  basis.  CO  is  neglible  except  in  the  fourth  case,  which  is 
in  question  otherwise. 

Calculated  0P  content  in  original  air.  Dry  basis.  Should  equal  21.07. 
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1700°F)  passed  through  a choker. 
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air  was  entrained) , 

and  through 
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out 

18  ft  of  a 

prior 

to  sampling  [at  about  933 
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he 

~i  p 

in.  duct. 

3 in. 

downstream  was 

about  -47  K 

. 

The 

se  .: onditi 

ently 

intended  to  be 

indi cat ive 

0 i* 

the 

emiss ions 

atmo  s 

phere  behind  an 

aircraft . 

Result 

from  the 

from  the  Aeronutronic  grou{  are  given  ir.  Table  2-3.  Note  that 
both  I'Y  and  IR  techniques  implied  a destruction  of  N(  In  sampling 
(even  after  passage  through  a choke  and  ft  of  on  which 

surface  reacti >ns  could  have  t a/  place).  Becaus<  nditions 
differed,  it  is  difficult  to  compar<  re  ilts  those  ob- 

tained on  the  J93  engine  by  AEDC . 


TABLE  2-3.  COMPARISON  OF  OPTICAL  AND  PRCBE  MEASUREMENTS 
FOR  NO  USING  A T-56  COMBUSTOR 


al . , 1976) 

Optical 

178 

114 


3 in.  downstream 
18  in.  downstream 


Optical/Probe 

1.63 

1.63 


Aeronutronic  (IR)  (Gryvnak  and  Burch,  1976) 

Optical/Probe 

Flat  profile  1.24 

Bell  profile  1.37 

Asymmetric  1.51 

A related  study  was  carried  out  by  Ly  : et  1 1 . , 1975;  see 
also  Lyon,  L975»  In  this  work,  a detailed  study  was  made  of 
emissions,  both  nonafterburning  and  from  afterburning  jet  engines 
using  J8!  : : J79  engines.  I Llutant  concentrations  were  ob- 

tained at  the  engine  exit  and  at  a variety  of  stations  to  60 
feet  downstream.  The  effort  used  prob<  techni  jues  nly , 1 ut 


these  were  used  under  the  widely  varying  conditions  between  the 
engine  exhaust  plane  and  points  downstream.  While  some  anomalies 
were  noted,  no  obvious  change  occurred  in  the  N0x  emission  index 
with  distance  downstream,  with  or  without  afterburning,  leading 
the  investigator  to  conclude  that  probe  sampling  gives  accurate 
results,  even  under  afterburning  conditions.  It  is  possible  to 
interpret  the  same  data  in  a somewhat  different  fashion,  as  is 
.one  in  this  report  in  Appendix  A,  with  results  that  suggest, 
but  do  not  prove,  that  reactions  were  taking  place  in  the  probe 
system,  at  least  in  some  of  the  runs. 

In  another  investigation  carried  out  in  197^-1975  (Few  et 
al.,  1977),  the  validity  of  the  UV  analytical  technique  was 
studied  by  comparing  NO  measurements  using  UV  and  conventional 
methods  on  gases  drawn  through  a probe;  results  were  also  ob- 
tained in-situ,  on  the  hot  gas  stream. 

The  results,  as  indicated  earlier  and  as  shown  in  Fig.  2-1, 
showed  excellent  agreement  downstream  of  the  probe  by  conventional 
and  in-situ  spectroscopic  techniques,  but  the  values  sc  obtained 
were  about  one-sixth  those  found  in  the  free  stream.  The  impli- 
cation, if  the  UV  technique  is  correct,  is  that  some  five-sixths 
of  the  free  stream  NO  was  destroyed  in  the  sampling  probe.  The 
conditions  used  (about  600  K combustor  inlet  temperature,  55  psia, 
and  fuel/air  ratios  of  0.01  to  0.05)  correspond  roughly  to  full 
power  at  about  mach  0.9  and  3 6,000  ft,  and  were  not  exceptional; 
i.e.,  afterburning  or  high  altitude,  high  mach  number  flight  was 
not  simulated.  Furthermore,  since  a 6-fold  ratio  was  found  at 
the  combustor  exist,  but  a much  smaller  ratio  implied  at  engine 
exit  by  the  J93  data,  it  could  be  speculated  that  much  NO  cre- 
ated in  the  combustor  is  being  destroyed  in  the  turbine. 

2.1.3  Related  Investigations 

The  possibility  that  probe  sampling  of  hot  reactive  gases 
can  give  erroneous  results  has  been  recognised  for  many  years, 
and  was  almost  certainly  considered  when  developing  ARP1256. 

'V  Cl 

<-  J 


J 
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Nevertheless,  recently  obtained  results — in  addition  to  those 
cited  above — Imply  that  much  remains  to  be  learned.*  A brief 
review  follows. 


Samuelson  and  Harman  (1977)  and  Benson  and  Samuelson  (1977) 
review  past  work  and  the  possible  homogeneous  and  heterogeneous 
reactions  that  can  take  place  in  sampling  combustion  gases  con- 
taining nitrogen  oxides  (NO  and  N09),  and  report  on  experiments 
involving  mixtures  of  NO  (500  ppm)  and  NO,,  (75  ppm)  in  an  Np 
carrier  gas  containing  various  other  gases  (Op,  C0n,  Hp,  CO, 

CH^,  CpH^,  C^Hg),  with  two  probe  materials  (316  stainless  steel 
and  silica).  A temperature  range  from  25  to  400°  C was  investi- 
gated. The  following  results  seem  to  be  of  particular  interest 
here : 


3. 


In  general,  silica  probes  show  less  reactivity  than 
stainless  steel  probes. 

In  many  of  the  tests,  NOp  was  converted  to  NO,  but  the 
sum  of  the  two  remained  unchanged. 

In  0,--free  mixtures  containing  Hp  or  C^Hg  (but  generally 


, v/u,  or  CpH^),  NO  was  removed,  in  some 


not  with  CH^,  CO. 
cases  quantitatively,  depending  on  reducing  gas  concen- 
tration, at  temperatures  of  400°C  or  lower,  particularly 
over  stainless  steel.  The  concentrations  used  were  500 

and  1000  ppm  C0Hr,  and  0.5  to  3 percent  H~. 
jo  • 

These  results  are  of  considerable  interest,  indicating  as 
they  do  the  possibility  of  reactions  taking  place  on  the  probe 
tips  or  in  the  sample  lines.  Of  course,  the  mixtures  used  were 
not  representative  of  typical  combustor  gases,  lacking  excess 
oxygen  and  combustion  moisture.  (Results  with  CO  might  have 


* 

Thus,  while  MO  at  room  temperatures  is  known  to  be  thermo- 
dynamically unstable,  its  kinetic  stability  is  very  great. 
Nevertheless,  it  can  decompose  slowly  heterogeneously . Re- 
cently published  data  (Robertson  et  al.,  1977),  for  example, 
show  that  NO  in  N2  (250  to  1000  ppm)  decomposes  by  3.3  percent 
in  6-10  months  in  steel  cylinders  and  by  3.8  percent  in  17 
months  in  commercial  cylinders. 
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been  different  had  H„0  been  present,  as  water  gas  shift  reactions 
could  have  resulted  in  Hp  reactions  with  NO.  Whether  these 
would  take  place  in  oxygen-containing  streams  is  unknown.) 


Effects  of  probe  materials  and  probe  designs  have  also  been 
studied  recently  at  AEDC  (R.  Bryson).*  The  probe  designs  used 
are  shown  in  Fig.  2-2.  An  AVCO-Lycoming  AGT-1500  combustor  was 
used.  This  combustor,  5-5  inches  in  diameter  and  12  inches 
long,  was  designed  for  use  in  a tank  gas  turbine,  is  symmetric 
in  configuration,  has  a well-defined  flow  field  and  was  thus 
considered  particularly  suitable  for  these  tests.  Results  are 
shown  in  Table  2-4.  Their  results  showed  the  following: 

1.  The  quick-quench  probe  as  recommended  by  GE  (Lyon  et  al. 

1975)  gives  lower  NO  values  than  the  AEDC  baseline 

probe,  which  latter  involves  a large  sample  tube  and 

a noncritical  pressure  drop. 

2.  Baseline  probe  results  were  independent  of  tip  material. 

3-  Strange  results  were  evident  for  CO.  High  CO  values 

were  associated  with  low  NO  . CO  values  determined  by 

x 

probe,  as  well  as  MO  values,  thus  appear  to  be  in  doubt. 


roughly  from  Few  et  al.  (1977),  as  being  about  175  ppmv  at 
0.01,  605  ppmv  at  0.02,  and  860  ppmv  at  0.03.  The  highest  values 
found  by  probe  at  corresponding  conditions  (39  ppmv  at  0.01,  133 
ppmv  at  0.02,  228  ppmv  at  0.03)  were  always  much  lower  than  those 
estimated  from  UV  spectrometry.  It  follows  that  probe  design 
is  important,  but,  particularly  in  light  of  independence  of 
probe  material,  does  not  appear  to  explain  the  differences  be- 
tween the  two  methods. 

In  another  AEDC  study  (Few,  1977)  of  particular  interest 
to  the  subsonic  aircraft  case,  optical  NO  measurements  were  made 
of  the  exhaust  from  an  F-101-GE-100  engine  at  simulated  altitudes, 


Private  communication  October  11,  1977-  See  Bryson  and  Few,  1978 


TABLE  2-4. 


COMPARISON  OF  EMISSIONS  MEASURED 
WITH  THE  VARIOUS  PROBES 


i. oio 

j. ulO 
O.Ol  I 


0.030 

3.030 


0.024 

O.G2*> 


BLC 


NO 

N°x 

CO 

CO- 

Air  H '.'rt . 

Con>UuStor 

San  pie 

Ui 1 uent 

Content . 

1 ontent , 

Content, 

Coritt  nt , 

■ : 

nessure,  kP 

a Pr«oe 

Pressure  mho' 

• ■ : t 

PPIIW 

ppmy 

Perc<  ’ 

J.dlb 

265 

BLSS 

1.04 

30 

39 

165 

266 

BLQ 

1 . Oo 

- 

33 

155 

2. 79 

268 

BLC 

1 .04 

- 

32 

43 

165 

2.77 

265 

00 

l . -33 

26 

33 

1 16 

2.07 

23 

• 

263 

4 . 96 

* 

25 

32 

126 

2.07 

16 

0.816 

262 

OIL 

4.87 

29 

34 

131 

OIL 

2 .57 

1.43, 

26 

32 

138 

r if. 

6 

OIL 

2.  30 

2.32,  V 

13 

13 

70 

1 .03 

47 

l.  *16 

OIL 

1.50,  A‘ 

•4 

97 

? 03 

• 

?6b 

OIL 

2.37 

1 . '6 , A 

11 

13 

62 

.83 

15 

0.816 

263 

OIL 

3.00 

1 .7f  , He 

22 

37 

113 

2. 32 

16 

263 

OIL 

2.23 

3.36,  He 

31 

38 

131 

2.14 

4 

• 

311 

BLSS 

1.06 

- 

114 

120 

70 

5.25 

310 

BLQ 

1.12 

- 

114 

129 

pc. 

6 ?7 

0.RI6 

311 

BLC 

1.06 

- 

118 

133 

70 

5.33 

3.82) 

311 

QQ 

2.29 

- 

73 

92 

77? 

4JT2 

1 1 

309 

5.M 

- 

71 

or 

337 

4 . 02 

11 

. I( 

OIL 

73 

96 

a.  ic 

c 

. ■ 1 

309 

OIL 

3.01 

1.36, 

63 

83 

41  7 

4.11 

4 

0.8*6 

DIL 

2.44 

2.24,  V 

44 

60 

315 

3.21 

4 

J.616 

308 

OIL 

2.76 

1 £ 

59 

75 

46C 

3 

0.516 

303 

DIL 

2.01 

i . <4  , - 

39 

51 

T 35 

J.dlf 

308 

DIL 

3.52 

1 . 5t , He 

91 

574 

3 

0.816 

308 

DIL 

2.60 

3.09,  tie 

63 

4 29 

I:  ! 

34  7 

ELSS 

1.07 

- 

192 

205 

124 

6.6} 

0. dl 6 

346 

BLQ 

1.15 

. 

1 

221 

165 

7 n< 

0.816 

34  7 

6LC 

1.03 

. 

6.  )3 

0.cl6 

346 

QQ 

2.00 

_ 

128 

16.3 

1878 

5.91 

0.816 

345 

QQ 

5.21 

. 

134 

163 

5.72 

U 

0.816 

345 

QQ 

6.61 

• 

136 

164 

324' 

5.  72 

14 

3.816 

343 

OIL 

6.72 

1 

1 7t 

6.36 

0.J16 

341 

DIL 

i.  54 

1.38,  N, 

124 

147 

3975 

32 

0.816 

342 

DIL 

2.74 

- -.c  •/ 

110 

13;< 

6759 

5 26 

74 

..  21 

34  3 

DIL 

3.21 

1.5  , A r 

114 

130 

4M0 

5 25 

25 

*.6lh 

342 

Ml 

2.76 

?.02,  A 

101 

116 

39P5 

4.33 

49 

0.8)6 

343 

OIL 

4.11 

' . ' ' , He 

141 

174 

J.816 

341 

ML 

3.06 

3.23,  He 

128 

5ioi 

5 75 

O.dc’l 

3 27 

BLSS 

1.07 

. 

169 

177 

6 '7 

. • 

RLO 

1.15 

1 75 

191 

85 

c 44 

Q 

. . i 

328 

BLC 

- 

174 

194 

75 

6 38 

(1.816 

1.92 

. 

106 

136 

856 

6.20 

M 

3. 31 

106 

134 

1052 

5.16 

0.816 

6.23 

- 

102 

132 

1148 

U HO 

0.816 

329 

OIL 

6 .29 

126 

147 

2400 

5. 32 

10 

' 0. « 1 

329 

3.  .’4 

1.42,  U2 

102 

13? 

2124 

6. 23 

14 

0.816 

OIL 

2. 59 

2.29.  *< 

112 

1848 

4.74 

26 

• 8 1 6 

OIL 

2 . 99 

1.48,  A‘ 

94 

11; 

1891 

4 . 79 

0.816 

3?9 

OIL 

1.94,  A 

33 

09 

1696 

4.26 

14 

DIL 

3.33 

1.66,  He 

115 

14: 

2333 

5.43 

sir 

OIL 

2. 93 

3.03.  He 

102 

134 

2078 

5. 1 5 

16 

Basel tne 

with  sta'nic-ss 

steei  surface 

59  . 

luicr  quench 

r.asel  !nr 

with  quart/  Su* 

•face 

• Uiluttor 

Basel ine 

wi  th  c >pper  sur 

■face 

,'res  St.re 

ratio  is  r » t * o 

of  combustor 

n . • • • - 

> 1 *ne  pressure 

the  v ; ji-P 


tne  volur«? 


2-10 


r 


under  subsonic  conditions,  and  as  a function  of  distance  from 
the  centerline.  Unfortunately,  probe  results  were  not  taken  in 
the  same  test,  but  extrapolations  were  made  from  probe  results 
(an  orifice-type  probe)  taken  on  a similar  engine  earlier.  The 
results  are  shown  in  Fig.  2-3.  Note  that  a rapid  drop-off  in 
NO  occurs  with  distance  from  centerline.  The  apparent  agree- 
ment between  probe  values  and  optical  values  at  larger  radii  is 
perhaps  misleading  due  to  the  extrapolation  procedure  used; 
in  addition,  mass  flux  as  a function  of  radius  is  unknown,  so 
that  it  is  impossible  to  obtain  a reliable  estimate  of  total 
NO^  flux  by  the  two  methods,  although  this  would  be  of  interest. 
If,  however,  uniform  mass  flux  per  unit  area  is  assumed,  and 
the  plot  accepted  as  drawn,  the  mass  weighted  emissions  by 
optical  techniques  would  be  about  35-^0  percent  above  these  by 
probe . 

An  optical  procedure  was  also  used  in  an  attempt  to  measure 
NO ^ • No  discernible  absorption  was  found,  indicating,  for  the 
setup  used,  less  than  20  ppmv  of  NO2*  It  was  concluded  that 
most  of  the  NGx  is  as  NO. 

The  above  results  suggest  significant  difficulties  in  using 
probe-sampling  techniques,  although  there  is  no  proof  that  the 
spectroscopic  values  are  correct;  there  are  other  additional, 
somewhat  fragmentary,  data  which  tend  to  confirm  these  diffi- 
culties and  indicate  the  need  for  further  work.  Thus,  Seery 
and  Zabielski  (1977),  in  reporting  on  work  carried  out  to  study 
chemistry  in  flames  (rather  than  on  measurement  methods  per  se), 
give  a comparison  of  NO  content  in  a flame  formed  by  methane 
and  air  ( stoichiometric  at  0.1  atm),  as  measured  by  molecular 
beam  samples  (with  mass  spectroscopy)  and  by  a quartz  micro- 
probe. The  microprobe  samples  were  found  to  show  far  lower 
(factors  of  5-10)  NO  values  than  those  found  by  the  molecular 
beam  system. 
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FIGURE  2-3.  NO  concentrations  and  static  temperature  prefiles.  Intermediate 
power.  The  NO  concentration  values  shown  are  at  10.2  cm  down- 
stream. The  probe  values  shown  have  been  extrapolated  (with  un- 
certainties introduced)  from  measurements  at  other  conditions. 
The  temperature  profile  shown  is  at  the  nozzle  exit. 

Source:  Few  (1977) 
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A related  investigation  is  underway  in  which  measurement 
methods  are  at  issue.  N conclusive  results  are  yet  available, 
but  soma  ; ns  il  th  ic  curacy  f the  ultraviolet  spec- 

troscopic technique  as  developed  at  AEDC  have  arisen 
(M.F.  Z bielski , private  communi cation,  1973 ) ; it  does  not 
appear,  howev<  r,  th  unc<  1 ainl  t issue  would  account 

ix-f  . as  found  at  AEDC  in  some  investiga- 

• i Lg.  -1.  Phi!  preliminary  finding  may  be  impor- 
tant and  may,  I n c dth  j r be  s id±«  s,  indicate  that 

: tl  sc  pic  results  have  been  in  some  error 

...  luestions,  there  is  alsc  recent  evidence 

t hat  het er  • s may  enter , wh Lch  could  be  involved 

on  un  Led  . lpling  afl  rt  xri  ing  mixtures  or  could 

take  place  t h s : art Lcles  carried  in  the  hot  stream.  Thus, 
Myersi  n (1978)  found  that  : ' : ■ • ■■  f carbonaceous 

deposit,  highly  re  ictiv  NC  decomposition,  could  be  formed 

on  tnc  walls  fa  eramj  r icl  r be  ing  used  in  the  study  of 
homogeneous  re  duel  . of  NO  by  hydrocarbons  (as  in  auto  exhaust 
mixtures)  at  i ' : • ruer  " 1 '-00-lo00oC . The 

oxygen  in  the  M appe  red  a:  !0,  the  nitrogen  as  Nn.  The  input 


gases  in  the 

■ xperiments  typica 

lly  involved  mixtur 

■es  in  N2 

con- 

taining  0.9  r 

lercent  i-C  H,  , , 0 . 

M 10 

4 percent  0^ » and  M 

JO  at  the 

level 

o 

o 

o 

o 

pi 

These  conditions, 

of  course,  do  not 

exist  in 

usual 

jet  engine  stro  ms,  although  soot  is  always  present  to  some 
l<  ree.  Myerson  also  refers  to  earlier  work  which  showed  that 
raphite  r sisl  mt  1 attack  by  x;  ;en  w uld  react  with  NO  (at 
1233  K) . The  homogeneous  reactions  are,  of  course,  also  of 
int ere st,  with  Samuelson's  work  already  cited.  Presumably, 
at  al  >g  processes  were  involved  in  work  of  Sugiyama  et  al. 
(1973)  which  showed  that  fuel-bound  nitrogen  is  converted  to  N2 
in  thf  ; resence  of  small  amounts  of  02  during  pyrolysis  at  600- 
800°C . 

Finally,  in  still  another  investigation,  two  sources  of 
error  Ln  heird  imii  metl  ds  have  been  studied  by  Mathews 
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et  al.  (1977).  These  have  to  do  with  relative  quenching 
clencies  of  different  third  bodies — water  vapor  showing  a : : 
quenching  efficiency,  for  example.  It  appears  that  err^  rn  f 
20  percent  or  so  may  be  involved  with  different  mixtures 
combustion  gases.  In  addition,  the  possibility  exists  '‘.ha- 
other  nitrogen-containing  gases  may  be  converted  to  I JO  in  th 
analyzer.  Obviously,  the  questions  involved  merit  attention 
in  any  measurement  program  using  these  techniques. 

2.1.4  Significance:  A Probable  Upper-bound  Estimate  of  MO 
Emission  Index  of  Concorde 

The  significance  of  the  NO  measurement  uncertainty  to 
obviously  relates  to  the  question  of  how  many  aircraft  woul i 
cause  a specified  effect.  What  would  be  needed  in  order  to 
bound  such  effects  are  measurements  using  a validated  "V  s:  - - 

trometric  technique  (if  probes  are  shown  to  be  unreliable- 
existing  engines  under  appropriate  operating  conditions  rather 
than  on  combustors  or  other  research-type  setups).  In  the  ace 
of  such  information,  the  significance  must  be  determined  in- 
directly. In  doing  so,  it  must  be  considered  that,  at  present 
both  probe  and  UV  spectrometer  data  are  suspect,  and  ratios  of 
the  two  techniques  cannot  necessarily  be  applied  from  engine 
to  engine. 

The  first  comparison  of  interest  is  that  of  the  emission 
index  data  derived  from  the  T-56  combustor  test  (Few  et  al. , 
1977)  to  that  obtained  from  engine  tests  (Davidson  and  Dor:  ., 
1973)  both  at  AEDC,  as  shown  in  Table  2-5-  The  nonafterburr.':, 
case  at  military  power  is  used.  Here,  the  engine  tests,  ■ . 
when  conducted  under  considerably  more  severe  conditions,*  >-.:v 
¥ 

According  to  Davidson  and  Domal  (1973)  » the  N0X  emission  hhx 
increases  by  the  2.5  power  of  the  combustor  inlet  teraperaturf  . 
The  emission  index  for  the  T-56  combustor  at  800  K inlet  vn  . : - 
by  this  relationship,  be  twice  that  at  600  K,  or  about  l'"  gm 
NOo/kg  fuel  by  UV  techniques,  using  the  data  in  Table  . . !. 

value  (100)  is  about  half  the  equilibrium  value  for  dry  air, 
this  inlet  temperature,  at  an  equivalence  ratio  of  0.9 
(Blazowski  et  al.,  1975). 
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lower  emission  indices  by  spectroscopic  techniques  than  were 
obtained  at  the  output  of  the  combustor  for  the  combustor  tests 


r 


using  the  same  technique.  This  would  seem  to  imply  that  con- 
siderable NO^  reduction  takes  place  in  the  turbine  (which  seems 
doubtful ) . 

TABLE  2-5.  NO  EMISSION  INDEX  DATA  COMPARISONS 
x 


T-56  Combustor  Test  (Few  et  al.,  1977) 


Temperature 

Pressure 

o 

X 

m 

I.  (as  NO^)9 

Inlet,  K 

Inlet,  psia 

f/a 

(Probe) 

(~ln-Situ) 

598 

56.4 

0.010 

4.8 

28.0 

598 

46.2 

0.020 

8.3 

48.9 

600 

54.4 

0.030 

10.4 

46.8 

604 

52.1 

0.040 

8.9 

50.3 

605 

57.0 

0.051 

3.6 

52.3 

J 93 

Engine  Test 

(Davidson 

and  Domal  , 

1973) 

NO  E.I.  (as  NO,) 
x v 2' 


i enper  d tur  te 
Inlet, 

r r e b b u r e ^ 

Inlet,  psiau 

f/a 

(Probe) 

( In-Si tu )c 

Comment 

590 

82 

0.0211 

4.60 

6.04 

Mach  1.4, 
35,000  ft 

672 

53 

0.0196 

6.70 

14.90 

Mach  2.0, 
55,000  ft 

800 

59 

0.0150 

12.03 

35.9 

Mach  2.6 

66,000  ft 

aBased  on  NO  alone.  NO,  data  are  not  reported  but  should  be  negligible. 

K 

Trom  Fig.  12,  p.  44  (graphical)  Davidson  and  Domal  (1973). 

r 

These  values  are  derived  from  the  reported  NO  ratios  and  the  N0X  emission 
indices  from  probe  measurements,  assuming  the  NOo  component  is  constant. 

As  another  comparison  of  interest,  the  J93  data  can  be  com- 
pared to  the  Concorde  data  (CIAP  Monograph  2,  pp.  4-34),  as  in 
Table  2-6.  Note  that  the  J93  probe  data,  if  converted  to  con- 
ditions used  in  the  Concorde,  agree  well  with  Concorde  probe 
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data  (16.7  vs.  16.8).  Also,  if  the  spectroscopic  data  are  cor- 
rect, and  emission  index  (E.I)  of  about  50  gm  NO^/kg  fuel  might 
be  assigned  to  the  Concorde  at  cruise.  The  agreement  in  probe 
results  is  somewhat  surprising  in  that  the  probe  data  obtained 
by  AEDC  were  apparently  somewhat  (~  10  percent)  low  compared  to 
other  turbojet  data.  (See  Davidson  and  Domal,  1973,  p.  52.) 


TABLE  2-6.  COMPARISON  OF  CONCORDE  (PROBE)  AND  J93  N0x  RESULTS 

N0*  E'L 

1 f 1 1 D lArtr  f*  1 1 T m c +•  I 


Temperature 
Inlet,  K 

Pressure  In- 
let, psia 

X 

(Test  No.); 

Condi tions 

f/a 

(Probe) 

(uv) 

Concorde 

801 

102 

0.0147 

16.8 

(4)  Mach  1.99 
58,879  ft,  ISA-9  K 

824 

95 

0.0141 

17.7 

(7)  Mach  2.00, 
58,000  ft,  ISA+5  K 

J93  (measured) 

800 

59 

0.0150 

12.0 

35.9 

Mach  2,  55,000  ft 

J93  (converted  to  Concorde  conditions)* 

800 

102 

0.0150 

16.7 

49.9 

Estimated  values,  assuming  P^  dependence  based  on  results  presented  by 
Blazowski  et  a 1 . , 1975;  prorating  the  quoted  results  at  59  psia. 

2.2  FLEET  FORECASTS  OR  PROJECTIONS 

Fleet  forecasting,  the  problems  inherent  in  it,  and  in  the 
fleet  figures  generated  thereby,  have  been  discussed  in  earlier 
reports  (see  Oliver  et  al.,  1 9 7 7 » Pozdena  1976  , Athens  et  al . , 
1976,  CIAP  Monograph  2,  1975)  and  will  not  be  discussed  in  de- 
tail here.  Nevertheless,  it  is  of  some  interest  to  consider 
two  papers  which  have  become  available  since  our  previous  report. 
One  of  these,  Aviation  Futures  to  the  Year  2000  (published  by 
the  Department  of  Transportation,  FAA,  February  1977),  was 
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sponsored  by  the  FAA;  the  other,  Geddes  (1977),  is  a projection 
by  McDonnell-Douglas , updating  similar  earlier  publications. 

In  Aviation  Futures , five  scenarios  are  assumed,  considera- 
tion being  given  "almost  exclusively"  (not  further  defined)  to 
the  U.S.  fleet.  The  scenarios  are  based  on  sets  of  assumptions 
as  to  population,  business  productivity,  cost  of  domestic  crude, 
etc.  Certain  data  extracted  from  the  report  are  given  in  Table 
2-7.  Note  that  SST  traffic  is  included  only  in  the  most  opti- 
mistic "expansive  growth"  scenario.  Note  also  the  assumed 
prices  of  oil:  in  the  "muddling  through”  scenario,  with  $14 
oil,  the  jet  fuel  consumption  is  less  in  the  year  2000  than  at 
present . 


TABLE  2-7.  UNITED  STATES;  YEAR  2000  ESTIMATES* 


Key  Variables 

1974 

Limited 

Growth 

Muddling 

Through 

Resource 

Allocation 

Individual 

Allowance 

Expansive 

Growth 

Population  (millions) 

212 

250 

297 

250 

250 

297 

Gross  National  Product 
101J  1973  dollars 

1.3 

1.9 

2.1 

2.9 

4.1 

4.3 

Cost  of  Oomestic  Crude 
1973  dol lars/barrel 

11.0 

8.5 

14.0 

8.0 

6.0 

7.5 

Enplaned  Passengers 

Mi  1 1 ions/yr. 

207 

406 

272 

471 

788 

1113 

Total  Revenue 

Passenger  Miles 

Bill  ions/yr. 

131 

259 

167 

304 

485 

597 

Jet  Fuel  Consumption 
Millions  of 
barrels/yr. 

190 

317 

158 

317 

517 

850 

Carrier  Fleet 

Total  Aircraft 

Mix  (percentage) 

2420 

3790 

2430* 

3540 

5660 

7670 

747-type 

5 

4 

8 

9 

10 

707- type 

0 

4 

0 

0 

0 

DC -10  type 

18 

12 

20 

17 

20 

727- type 

62 

64 

59 

55 

41 

Small  (50 
passengers) 

15 

16 

10 

11 

9 

SST 

0 

0 

0 

0 

3 

Large  (1000 
passengers) 

0 

0 

3 

3 

4 

Jet  ST0L  (150 
passengers) 

0 

0 

4 

5 

13 

Source:  FAA,  1977 


*2670  in  198S 
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Direct  comparison  of  the  Aviation  Futures  results  with  those 
of  McDonnell-Douglas  (Decides)  is  difficult,  as  time  periods 
differ,  some  questions  arise  as  to  what  is  included  in  the  two 
forecasts,  and  there  is  a difference  in  focus,  one  (Geddes) 
emphasizing  passenger  miles  and  the  other  passenger  enplanement 
and  operations.  Nevertheless,  some  figures  have  been  developed 


for  comparison  of  the  twro  sources, 

as  shown 

in  Table 

2-8. 

TABLE  2-8.  FLEET 

COMPARISONS 

1985 

1990 

2000 

Total  Aircraft 

Aviation  Futures  ("Expansive  U.S.") 

3890 

-- 

7670 

McDonnell-Douglas  (world) 

7217 

7908 

-- 

McDonnell-  Douglas  (U.S.)* 

3115 

3270 

-- 

SSTs 

Aviation  Futures  ("Expansive  U.S.")  - 

- 20 

( 100)2 

230 

McDonnell-Douglas  (world) 

16 

118 

-- 

U.S.  Passenger  Miles/yr.  (Billions) 

Aviation  Futures  ("Expansive  U.S.")  ~ 260J 

* 3503 

5973 

McDonnell-Douglas 

302 

- 380 

— 

Prorating  U.S.  numbers  by  U.S.  share  of  world  schedule/passenger-miles. 

2 

Interpolated 

3 

In  this  scenario,  a rapid  growth  in  jet  STOL  takes  place  between  1985 
and  2000,  reducing  average  passenger-miles  per  enplanement.  The  approxi- 
mate figures  shown  were  estimated  by  this  author,  assuming  630,  585,  and 
540  (their  numbers)  miles  per  trip  average  in  the  respective  years. 

The  McDonnell-Douglas  projections  show  a new  SST  introduced 
in  1986.  The  Avi  ation  Futures  document  does  not  indicate  the 
SST  types.  Aviation  Futures  shows  a 1000-passenger  plane  to  be 
in  operation  (~  40)  in  the  expansive-growth  scenario  by  1985; 
McDonnell-Douglas  does  not  include  such  a craft. 
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2.3  CONCORDE  EMISSIONS  ESTIMATES 


In  a previous  report  (Oliver  et  al.,  1977),  NO^  emission 
estimates  were  provided  for  the  Concorde  aircraft  as  well  as 
for  a "1990-High"  estimate  of  a world  fleet.  No  new  world  fleet 
results  are  available,  but  new  estimates  have  been  provided  for 
the  Concorde  by  P.C.  Beadle  of  the  British  Aircraft  Corporation 
(Beadle,  1977).  The  new  estimates  are  for  only  a nine-aircraft 
fleet,  but  cover  a wide  network,  and  are  the  only  estimates 
available  for  the  full  altitude-latitude  range.  The  data  are 
still  subject  to  adjustment,  as  certain  approximations  have  been 
made  in  treating  the  various  atmospheres  at  different  latitudes. 
Nine  aircraft,  of  course,  are  of  little  environmental  signifi- 
cance, but  the  distribution  may  be  of  interest  to  traffic  mod- 
elers . 

The  network  Beadle  studied  is  that  planned  "for  likely 
operation"  in  1978-1979-  The  route  network  shown  in  Table  2-9 
is  served  by  nine  aircraft,  totaling  57-5  hours  daily  usage, 
from  6. A hours  per  aircraft  per  day  (or  2330  hours  each  per 
year) . 

The  NO  (as  N0„)  emission  index  data  used  by  Beadle  is 
x 2 

shown  in  Table  2-10.  These  were  generated  by  him  and  involve 
a much  larger  correction  to  zero  humidity  than  used  by  Scott 
(197*0  or  recommended  by  Blazowski  et  al.  (1975),  so  that  the 
emission  index  values  at  cruise  are  20-40  percent  higher  than 
given  previously  by  Scott.  Beadle  argues  that  his  correction 
gives  a better  fit  to  available  data  for  supersonic  flight  con- 
ditions.* Other  emission  indices  recommended  by  Beadle  are 
given  in  Table  2-11;  the  data  are  from  COMESA  (1975) - 


Beadle's  correction  seems  excessive  in  view  of  the  data  re- 
ported by  Scott  where  an  80  percent  reduction  in  inlet 
humidity  resulted  in  an  increase  from  17.7  to  18.4  to  19.3 
kg/N02/gm.  Scott  corrected  for  the  remaining  20  percent 
water  with  a 3 percent  further  increase.  Beadle's  (continued) 


TABLE  2-9.  CONCORDE  ROUTE  NETWORK  SUMMARY 


Route 


..  ij  r ^ e . 


1 'J  7 7 


Or  i y i r. 

Destination 

Distance, 
nm  i 

Flight  Ti me , 
minutes 

Fuel, 
k n 

Wee 
F ren 

London 

Washington 

3,306 

222 

73,591 

7 

New  York 

3,060 

203 

69,385 

14 

Bahrain 

3.056 

222 

75,115 

4 

Washington 

London 

3,423 

209 

77,945 

7 

Paris 

3,546 

215 

79,853 

7 

New  York 

London 

3,115 

1 89 

68  ,059 

1 4 

Paris 

3,266 

201 

69,964 

14 

Bahrain 

London 

3,070 

235 

73,709 

4 

S i ngapore 

3,624 

235 

77,785 

3 

P a r i s 

Washington 

3,464 

226 

80,971 

7 

New  Y o r k 

3,221 

208 

76,774 

1 4 

Daka  r 

2,485 

171 

58,520 

2 

Santa  Maria 

1 ,474 

1 1 9 

35 ,760 

1 

5 i n na  po  re 

Bahrain 

3,040 

220 

79,866 

3 

Mel  bourne 

3,387 

221 

77  ,310 

3 

Ju  ka  r 

Paris 

2,512 

166 

54  ,260 

2 

Rio  de  Janeiro 

2,789 

130 

63,192 

p 

Santa  Maria 

Paris 

1 ,454 

1 1 1 

32,315 

1 

Caracas 

2,785 

1 75 

61  ,026 

1 

Melbourne 

Singapore 

3,353 

223 

78,994 

3 

i o d e Janeiro 

Dakar 

2,779 

1 73 

60  ,849 

2 

Buenos  Aires 

1,150 

88 

33,013 

2 

( iracas 

Santa  Maria 

2,793 

1 72 

60,195 

1 

iuenos  Aires 

Rio  de  Janeiro 

1,150 

87 

29,481 

2 

(1^70) 

(1975) 


correction  for  the  remaining  20  percent  of  water  is  25  per- 
cent or  more,  which  is  difficult  to  understand  in  view  of 
Blazowski's  theoretical  arguments  (Blazowskl,  1975). 

Elazowski  suggests  that  the  difference  in  NO  emission  index 
between  air  at  300  ppm  and  air  at  3 ppm  at  stratosphere  condi' 
tions  should  be  about  1/2  percent. 
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TABLE  2-10.  EMISSION  INDICES  FOR  CONCORDE  (g  NO„/kg  Fuel) 


Source:  Beadle,  1977 


TABLE  2-11.  CONCORDE  EMISSION  INDICES  (OTHER  THAN  NO  , SO?) 

Source:  Beadle,  1977 


Emissions  other  than  NO 

Cl 

Water  Vapor 
Carbon  dioxide 
Carbon  monoxide 
Unburnt  hydrocarbons 
Particulates 


Emission  Index  (g/kg  fuel) 

1270 

3140 

3.5 

0.2 

0.1 


Beadle's  analysis  (like  Scott's)  was  based  on  use  of  an 
International  Standard  Atmosphere  (which  is  the  same  as  the  U.S. 
Std.  Atm.  1976  for  purposes  here)  plus-or-minus  varying  incre- 
ments in  temperature,  with  the  standard  atmosphere  assumed  for 
the  troposphere  (to  13  km)  and  a uniformly  cooler  or  warmer 
atmosphere  above,  rather  than  standard  atmospheres  for  other 
latitudes.  Some  uncertainty  is  recognized  in  so  doing,  but  the 
resulting  estimates  should  still  be  suitable  for  modeling  pur- 
poses, in  view  of  the  other  uncertainties  involved. 

The  Concorde  is  limited  to  a pressure  altitude  of  1 8 . 3 km 
or  60,000  ft  (Scott,  197*0,  or  about  72  mbar.  The  geometric 
altitude  resulting  varies  with  season  and  latitude  from  about 
17*6  km  (arctic  winter)  to  about  18.6  km  (30°N  in  summer). 

Beadle  gives  fuel  flows  and  NO  (as  N0o)  emissions  from  0 
o x 2 

to  19  km  by  1-km  bands  and  from  U0°S  to  55°N  by  5°  latitude 
bands.  His  results  are  given  in  Tables  2-12  and  2-13- 

Beadle  also  analyzes  normalized  emissions  by  altitude  bands 
in  the  stratosphere  (15-16,  16-17,  17-18  km)  as  done  earlier  by 
Oliver  et  al.  (1977,  pp.  2-32)  in  adjusting  A.D.  Little  (Athens 
et  al.,  1976)  estimates.  Beadle  finds  Oliver's  adjusted  values 
agree  reasonably  well  with  his  estimates,  if  taken  to  apply  to 
the  Northern  Hemisphere  wherein  most  of  the  traffic  lies.  How- 
ever, surprisingly  large  relative  amounts  of  N0x  enter  the  18-19  km 


TABLE  2-12.  FUEL  FLOW  FOR  A NINE-CONCORDE  FLEET.  kg/WEEK.  Source:  Beadle,  1977 


Because  such  distribu- 


band,  differing  from  earlier  estimates, 
tions  are  sensitive  to  the  assumptions  made  in  treating  the 
atmospheric  temperature  profile,  such  discrepancies  may  not  be 
surprising  (or,  probably,  too  important). 

The  NO  emissions  shown  by  Beadle  are  almost  certainly  too 
small  to  be  noticeable  in  a computer  perturbation  study.  A 50- 
fold  increase  (corresponding  to  450  aircraft)  would  perhaps  be 
appropriate  for  modeling  purposes  for  accuracy;  however,  to 
test  possible  nonlinearities  (see  Section  3-5.5),  several  fleet 
sizes  would  need  to  be  examined. 

The  Beadle  results  just  shown  are  nominally  for  geometric 
altitude  models.  Resolution  is  probably  such,  however,  that 
little  additional  uncertainty  will  be  introduced  if  these  dis- 
tributions are  used  in  pressure  coordinate  models.  It  is  impor- 
tant to  recognize  that  a substantial  part  of  the  total  emissions 
is  released  in  climbing  to  altitude,  so  that  models  assuming 
aircraft  at  a single  altitude  may  give  misleading  results. 

2.4  N0x  EMISSION  REDUCTION  EFFORTS 

The  CIAP  (Grobecker  et  al.,  1974)  and  NAS  (1975)  studies 
indicated  that  both  subsonic  and  supersonic  aircraft  would  re- 
duce the  ozone  layer  if  fleet  growth  were  to  continue  as  then 
expected  unless  N0x  emission  indices  are  substantially  reduced 
(6-fold  or  more).  This  goal  thus  became  an  add-on  to  efforts 
which  had  been  underway  for  some  time  to  reduce  emissions  at 
takeoff  and  landing,  based  on  EPA  proposed  standards,  and  pro- 
grams were  initiated  under  NASA  sponsorship.  All  these  efforts 
were  reviewed  in  a conference  held  at  NASA  Lewis  in  May  1977 
(Aircraft  Engine  Emissions,  1977).  The  Stratospheric  Cruise 
Emission  Reduction  Program  (SCERP)  (Diehl  et  al.,  1977,  op.cit.) 
is  of  particular  interest.  Much  of  the  "clean  combustor"  work 
directed  to  EPA  goals  is  of  general  technical  interest,  but  the 
modest  reductions  achievable  by  the  combustor  modifications 


investigated  (2-fold  perhaps)  make  them  po.r  candidates 
the  hold  reductions  called  for  in  CIAP  and  by  the  NAS  (which 
! ’ o-  course,  were  set  by  the  stratospheric  chemistry  then 
ei  -t  • will  not  be  reviewed  here).  This  is  not  to  say 
lodest  goals , wit  their  perhaj  more  practi 
j will  not  eventually  find  use  ;ruise  : : iiti  ns. 
x values  rep  rted  in  the  conference  were  based  on  stan- 
probe-sampling  techniques. 

^ r ' • ii-s  a.  plot-  (tig.  2— 1|)  which  summarizes  the 

for  t tl  bs  nic  and  supersonic  cruise.  Several 
should  be  noted: 


'he  inapplicability  of  clean  combustor  technology  to 
n~  expressed  goals  is  clear. 


>STs,  t he  great  it  st  i n N ..  emiss  i n 


m ex  using  conventional  combustors  is  due  to  expected 
i:!  r'  combustor  pressures  and  temperatures,  presumably 
aljed  for  by  efforts  10  achieve  lower  specific  fuel 
msumr  *.  1 on  at  cruise.  Note  that  the  mach  2.32  air- 
• -0  f e 1 (It  km)  design  condition  flies 

an  lone  : . Jlean  combustor  tech- 
• ' ' i r ; 1 maintain  future  NO  emission 

X 

I • ' i •?«  • urrent  lev-  1 • . 


• • • re  irnnended  levels  ' f Me  emissions  are  approached 


rc  ti  n , " remix— 


• atalytic  :hniques.  rhe  prob- 

pr<  ti  engine  using  any  of  these 
niques  se vere ; t hose  us irig  a t a ly s u probabl 

- ' ■ u*  hybrid  m,  isin( 

’ t : . r ly  t promising.  he  L sys- 
h ■ • • 1 nve:? tig:  ted  to  the  greatest  degree. 


: ’ i witl  th  t - Lqu  i ide  jues- 

xi  • ashba  c,  variat  - met  y , auto- 


"'w  iiM-  limits  in  term: 


if  <5 


st 


1 metry ) , 


SUBSONIC  CRUISE  NOx  EMISSION  OUTLOOK 


M - 0.  85 
ALT  - 35  000  ft 


EMISSION 
INDEX,  10- 
gN02/kg  FUEL 


■APPROXIMATE 
RECOMMENDED  LEVEL 


CS-77-537 


CONV'L  CLEAN  FORCED  LPP  CATALYTIC 
COMB.  CIRC 


SUPERSONIC  CRUISE  NOx  EMISSION  OUTLOOK 


EMISSION 

INDEX, 

gNOjj/kg  FUEL 


CS-77-513 


CORE  EMISSION  LEVEL 
ENGINE  EMISSION  LEVEL 

M - 2.  32 
ALT  - 52  000  ft 


CURRENT  SST 


APPROXIMATE 

RECOMMENDED 

LEVEL 


CONV'L  CLEAN  FORCED  LPP  CATALYTIC 
COMB.  CIRC 


FIGURE  2-4.  NOx  emission  outlooks  according 
to  Diehl  et  al.,  1977 


and  effects  of  engine  transients  (relight  requirements,  compres- 


sor stall,  etc.)-  All  these  problems  (and  the  SCERP  program) 
are  in  a relatively  early  state  of  development.  Nevertheless, 
some  interesting  results  have  been  achieved.  Thus  (Fig.  2-5), 
very  low  N0V  emissions  can  be  achieved,  such  as  0.3  gm  NCu/kg 
fuel  at  the  inlet  conditions  shown,  at  an  equivalence  ratio  of 
0.1  (which  is  uncomfortably  close  to  blowout)  with  greater 
than  99  percent  combustion  efficiency.  It  is  also  of  interest 
that  pressure  has  little  theoretical  effect  on  NO  formation 
at  lean  conditions,  as  shown  in  Fig.  2-6,  in  contrast  to  the 
usual  0.5  or'  0 . 6 oo  --  d - ncy  found  with  conventional  com- 
bustors. The  fragmentary  experimental  data  are  not  very  illu- 
minating. 


The  operating  margin  question  and  fuel  air  preparation  prob- 
lem are  illustrated  by  Figs.  2-7  and  2-8.  Note  that  higher  in- 
l-.-t  air  temperatures  increase  the  operating  margin;  also  note, 
however  (Fig.  2-8)  that  the  probability  of  preignition  also 
increases.  A number  of  other  problems  are  also  discussed 
ana  illustrated  by  Diehl  et  al. 

Forced  circulation  techniques  are  discussed  by  Roberts  and 
But:. ; (1977)  in  the  same  document.  These  involve  use  of  lean 
premixed  and  prevapcrized  fuel-air  mixtures  with  a recircula- 
tion. none  for  flame  stabilization.  Studies  on  these  concepts 
were  begun  in  March  197^  and  were  in  the  final  stages  at  the 
time  of  the  Aircraft  Engine  Emissions  Conference  (May  1977) - 
The  two  concepts  studied  were  jet-induced  circulation  (JIC) 
and  vortex  air  blast  (VAB);  these  are  illustrated  in  Fig.  2-9- 
The  goals,  in  general,  were  low  NC'x  emissions  (1.0  gm/kg)  along 
with  low  CO  (0.5  gm/kg)  and  total  hydrocarbon  (0.5  gm/kg) 
emis.L  n indices  at  a combustor  inlet  temperature  and  pressure 

°R  and  5 atmospheres,  with  a combustor  outlet  temperature 
f 3<L  ~0°R,  using  Jet  A-l  fuel.  The  techniques  are,  in  effect, 
spec:,  f ' : mechanical  approaches  to  carrying  out  the  lean  premixed 
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FLAME  TUBE  EMISSIONS 
5.5  atm.  980 0 F 


NO, 

EMISSION 

INDEX, 

gNO?/kg  FUEL 


CS-77-549 


FIGURE  2-5. 


NOx  and  CO  emission  indices  under  lean  conditions. 
The  auto  ignition  and  blowout  limits  are  shown  for 
the  conditions  indicates.  "WSR"  refers  to  a "well 
stirred  reactor"  theoretical  result. 

Source:  Diehl  et  al.,  1977 


PRESSURE  EFFECT  ON  NO*  EMISSIONS  FOR  JET  A 


FIGURE  2-6.  Theoretical  and  experimental  pressure 
effects  in  N0X  formation  at  various 
equivalence  ratios  as  a function  of 
pressure. 

Source:  Diehl  et  al . , 1977 
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PRESSURE  - 5 atm 
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'SCERP  GOAL 
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FIGURE  2-7.  N0X  emission  index  at  various  equivalence 
ratios  and  inlet  temperatures,  showing 
stable  flame,  blowout,  and  operating 
margins . 

Source:  Diehl  et  al . , 1977 


JET  A AUTOIGNITION  DELAY  AND 
DROPLET  EVAPORATION  TIMES 


FIGURE 


COMPRESSOR  EXIT  PRESSURE,  atm  cs-77-533 


2-8.  Auto  ignition  regions  and  droplet 
evaporating  times  as  a function  of 
combustor  exit  pressure  and  droplet 
diameter. 

Source:  Diehl  et  al.,  1977 
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JIC  COMBUSTOR  DETAILS 
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FIGURE  2-9.  The  jet-induced  combustor  (JIC)  and  vortex  air  blast 
(VAB)  combustor  concepts. 

Source:  Roberts  and  Butze,  1977 
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prevaporized  combustion  processes  over  a range  of  : • 
ditions  (idle,  takeoff,  cruise).  In  the  JIC  system,  the  : r 
and  fuel  are  premixed  and  fed  (at  30°  relative  to  the  com:  . 
axis)  into  a combustion  chamber  in  a direction  counter  to  the 
eventual  flow.  The  VAB  system  involves  concurrent  flow,  wit! 
variable  angle  swirl  vanes;  reaction  stabilization  occurs  in 
a resulting  vortex  recirculation. 

Early  efforts  showed  that  low  NO  . emissions  at  cruise  con- 

X 

ditions  could  be  obtained,  with  the  VAB  concept  the  better  one 
from  an  emissions  standpoint  at  cruise,  and  essentially  meeting 
the  goal  of  1.0  gm/kg.  However,  as  was  noted  in  earlier  d’ s- 
cussion,  the  lean  stability  deteriorated  (and  CO  emissions  in- 
creased) rapidly  with  decreasing  inlet  temperatures,  as  shown 
in  Fig.  2-10,  so  that  the  simple,  lean,  premixed,  constant- 
geometry  combustor  cannot  be  considered  for  use  across  the  range 
of  engine  conditions. 

More  recent  work  with  these  concepts  has  been  concerned 
with  seeking  ways  to  avoid  the  problems  indicated  above,  using 
variable  dilution,  variable  geometry,  staging,  etc. 

The  foregoing  discussion  has  been  intended  to  highlight 
briefly  some  of  the  possibilities  and  problems  in  achieving 
ultra-low  N0x  combustors.  It  would  appear  obvious  that  con- 
siderable time  and  money  will  be  required  before  such  concepts 
can  be  brought  to  the  level  of  safety  and  sophistication  re- 
quired for  aircraft  use.  Stationary  gas  turbine  systems  or 
automotive  applications  using  these  approaches  might  well 
come  before  those  for  aircraft. 

2.5  SOME  COMMENTS 

The  foregoing  discussion  permits  only  a few  general  obser- 
vations. Thus,  uncertainties  in  emission  Indices  still  exist 
as  they  have  since  1973-  Aircraft  fleet  projections,  particu- 
larly of  S3Ts , have  dropped  dramatically.  Progress  is  being 


VAB  REACTION  ZONE  CO  EMISSIONS 


FIGURE  2-10.  CO  emissions  with  a vortex  air  blast 

(VAB)  combustor  simulated  idle  pressure 
with  various  inlet  temperatures,  illus- 
trating rapid  rise  in  CO  at  low  inlet 
temperatures . 

Source:  Roberts  and  Butze,  1977 


.xde  on  NO^  cruise  emission  reduction  programs  (although  the 
riginai  basis  for  these  efforts,  as  will  be  described  more 
'ully  in  the  following  sections)  has  been  shown  to  be  in  sub- 
rr.iai  doubt)  but  the  attendant  problems  in  achieving  large 
•eductions  in  NO„  emission  index  values  are  all  too  evident. 


3.  OZONE  MODELING  DEVELOPMENTS 


3.1  INTRODUCTION 

As  noted  in  Section  1.0,  the  cut-off  date  for  input  to  our 
previous  summary  report  under  the  HAPP  effort  (Oliver  et  al . , 
1977)  was  December  1976.  Since  that  time,  there  have  been  a 
number  of  developments  affecting  the  computed  effects  on  ozone 
of  aircraft  emissions,  chiefly  in  the  kinetics  area.  The  prin- 
cipal effort,  in  a national  sense,  in  this  period  has  been 
uevoted  to  the  halogen  problem;  the  fertilizer  question  has 
also  been  subjected  to  some  study.  The  sheer  volume  of  this 
other  work,  and  its  somewhat  peripheral  relevance  to  the  air- 
craft question,  preclude  its  detailed  review  here.  An  attempt 
will  be  made  to  draw  from  these  work  developments  which  seem  to 
be  of  particular  interest  to  aircraft  effects  and  to  the  uncer- 
tainties therein. 

As  will  be  evident  in  the  following  material,  changes  in 
chemistry  which  have  generally  reduced  computed  effects  of  NO 

A 

additions  have  brought  what  have  been  secondary  issues  to  the 
fore,  such  as  the  effects  of  thermal  feedback,  of  simultaneous 
water-vapor  injection,  and  of  detailed  choices  in  reaction 
rates;  the  sign  of  the  effect  on  the  o zone  column  change  can 
vary  with  choices  in  these  matters,  as  well  as  with  other  model- 
ing details.  A recently  passed  law  is  of  interest  in  this  con- 
nection and  merits  comment. 

3.2  SOME  COMMENTS  ON  PL95-95 

In  connection  with  modeling  "details,"  the  apparent  impact 
of  Public  Law  95-95,  passed  on  August  7,  1977,  should  be  noted. 
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This  act.  Clean  Air  Act  Amendments  of  1977,  covers  many  aspects 
of  air  quality;  however,  amendments  to  Section  126  are  of  prime 
interest  here.  This  amendment  adds  a Section  B to  this  section 
of  prior  legislation  entitled  "Ozone  Protection;”  the  amendment 
refers  only  to  ozone  in  the  "stratosphere".  The  stratosphere 
by  definition  (Sec.  152)  is  "that  part  of  the  atmosphere  above 
the  tropopause,"  and  thus  includes  the  mesosphere,  etc.  The 
act  states  that  the  "Administrator  [of  the  EPA]  shall  conduct 
a study  of  the  cumulative  effect  of  all  substances,  practices, 
processes,  and  activities  which  may  affect  the  stratosphere, 
especially  ozone  in  the  stratosphere"  and  goes  on  to  cite  halo- 
carbons,  other  sources  of  chlorine,  bromine  compounds,  and 
emissions  from  aircraft.  EPA  thus  becomes  the  lead  agency  for 
protection  of  stratospheric  ozone,  and  through  a Coordinating 
Committee,  reviews  and  comments  on  research  by,  and  research 
plans  of,  other  agencies.  An  interim  report  by  the  NAS  has, 
in  fact,  already  been  prepared  under  this  act  (NAS,  1977). 


While  the  full  implications  of  this  act  are  still  developing, 
an  apparent  (if  perhaps  temporary)  result  has  transpired  due  to 
emphasis  in  the  act  on  changes  in  stratospheric  ozone,  rather 
than  on  changes  in  the  ozone  column,  including  changes  in  the 
upper  troposphere,  which  total  effect,  of  course,  controls 
changes  at  the  earth's  surface.  This  emphasis  can  lead  to  awk- 
ward reporting  problems  for  investigators.  Thus,  in  connection 
with  climatic  effects,  release  of  halocarbons  may  cause  depletion 
of  ozone  in  the  stratosphere,  which  would  involve  a cooling  at 
the  earth's  surface;  however,  on  inclusion  of  the  tropospheric 
greenhouse  effect  of  the  halocarbons  themselves,  the  net  effect 
appears  to  be  a surface  warming.  Similarly,  1-D  atmospheric 
ozone  modelers,  if  the  act  is  taken  literally,  may  decide  not 
to  include  the  troposphere  in  their  models,  and  may  obtain 
thereby  misleading  results,*  particularly  if  aircraft  effects 
If 

It  should  be  noted  that  modelers  who  exclude  the  troposphere 
have  argued  the  opposite;  that  is,  tropospheric  chemistry,  in 
their  view,  is  so  poorly  understood  that  inclusion  of  the 
troposphere  can  lead  to  misleading  results. 
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are  being  computed.  With  aircraft  (see  e.g.,  Hidalgo  and  Crutzen 
1977),  particularly  realistic  fleets  with  principal  traffic  in 
the  upper  troposphere,  significant  changes  can  take  place  which 
may  be  more  important  at  ground  level  than  the  very  small 
change  (possibly  of  the  opposite  sign)  which  may  result  in  the 
stratosphere.  Any  model  excluding  the  troposphere  necessarily 
fixes  conditions  more  or  less  arbitrarily  at  the  lower  model 
boundary  (the  tropopause),  thereby  affecting  or  suppressing  the 
changes  which  would  take  place  at  this  altitude  in  a model  with 
a troposphere  included.  Erroneous  or  misleading  estimates  of 
ozone  effects,  like  those  for  climate,  could  result  using  such 
procedures . 

- > 

3.3  TRANSPORT  PROCESSES:  1-D  PARAMETERIZATION 
3.3.1  Background 

The  fundamental  problems  in  representing  atmospheric  trans- 
port processes  or  "dynamics”  in  a one-dimensional  parameteriza- 
tion (the  "Kz"  profile)  have  been  described  at  length  elsewhere 
(see,  e.g.,  Mahlman,  1975;  COMESA,  1975;  NAS,  1976,  Appendix  B; 

Oliver  et  al . , 1977).  Here  we  note  only  that  (1)  In  principle, 

* 

each  tracer  has  its  own  profile,  depending  on  its  sources  and 
sinks;  water  vapor,  e.g.,  very  likely  behaves  differently  than 
does  either  MO  or  ozone;  (2)  fundamental  ambiguities  as  to 

X 

interpretation  of  1-D  results  are  evident,  with  some  modelers 
arguing  that  their  models  represent  "global  (time  and  spa^e) 
averages"  and  others,  mean  mid-latitude  conditions,  etc.;  and 
(3)  for  sources  near  the  tropopause,  no  satisfactory  means  of 
taking  into  account  time  and  latitude-varying  tropopause  heights, 
or  regions  of  upwelling  or  subsidence,  has  been  developed.  In 
spite  of  these  and  other  problems,  however,  the  great  computa- 
tional simplicity  introduced  by  carrying  out  computations  in 
one-dimension  has  insured  their  continued  use.  Attempts  to  re- 
fine such  representations  are  of  continuing  interest,  as  are 
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intercomparisons  of  the  estimates  used  by  different  authors. 
Some  recent  developments,  with  relevant  background  information 
follow. 

3.3.2  N £ 0 Tracer  Studies  and  Derived  K z Profiles 

Undoubtedly,  the  most  extensive  recent  work  in  characteris- 
ing 1-D  transport  is  that  of  Schmeltekopf  et  al.  (1977,  1977a). 
These  authors  reported  on  an  extensive  series  of  measurements 
of  N as  a function  of  altitude  (to  about  35  km),  latitude 
(7 B°S  to  63°N)  and  (to  s ' ■ with  interpn  tati  n 

thereafter  as  to  K profiles.  The  w rk  ' : ■ r nt  ly  t y<  t 

fully  reported;  the  fol  i:  prii  ipal  ly  fr<  m Schmeltekopf 

et  al.  (1977). 

Schmeltekopf  et  al . (1977  ),  argue  that,  while  both  CH,.  and 
NpO  are  suitable  tracers,  the  loss  processes  involved  have 
greater  uncertainties  with  CH^ , so  that  N,;0  is  a more  suitable 
tracer.  Thus,  CkL  is  destroyed  principally  by  OH,  the  profile 
of  which  is  not  well  established,  as  well  as  by  CK^D)  (ar.d  by 
Cl,  although  not  noted  by  Schmeltekopf),  whereas  the  principal 
loss  process  (~  80  percent)  for  N^O  is  photolysis 

::£0  + h v ( A < 230  n m)  ->  N?  + 0(1D) 

with  secondary  destruction  mechanisms  (~  20  percent  from  0(XD) 
formed  by  photolysis  of  C., , giving 

•J 

N^O  + 0(^D)  -*■  2N0  (50  percent) 

or 

NpO  + O(^D)  N ^ + 0n  (50  percent) 

Schmeltekopf  et  al.  (1977)  state  that  the  overall  uncertainty 
in  the  photodissociation  rate  is  ± 50  percent.  (This  is  evi- 
dently less  than  the  overall  uncertainty  in  [OH]  concentration, 
on  which  only  limited  data  exist — see  Anderson,  1976). 
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Several  aspects  of  this  work  are  of  particular  interest: 

1.  Local  K„  profiles  vary  tremendously  with  latitude,  with, 
e.g.,  Alaska  values  being  up  to  two  orders  of  magnitude 
smaller  than  Panama  values.  In  the  authors'  words, 
"Individual  profiles  probably  have  little  physical 
reality . " 

2.  Inclusion  of  a Panama  profile,  of  which  only  one  was  aval 
able  at  the  time  of  the  1977  paper,  roughly  tripled  the 
global  average  eddy  diffusivity  otherwise  obtained,  in  th 
17-27  km  altitude  range  principally  discussed. 

3.  The  area-averaged  eddy  diffusivity  coefficients  derived 
were  found  to  be  a factor  of  1.5  to  3 times  those  found 
by  Hunter.  (1975)  using  CH,.  data  taken  at  Palestine, 

Texas.  This  result  relied  heavily  on  the  Panama  data. 

In  this  author's  opinion,  the  profiles  derived  by  pro- 
cedures, such  as  those  used  by  Schmeltekopf  et  al.,  seem  most 
appropriate  to  a material,  such  as  N„0  or  one  of  the  long-lived 

C. 

halocarbons,  which  involve  a surface  source,  with  long  life  and 
good  mixing  in  the  troposphere,  and  a stratospheric  sink.  To 
what  extent  the  mean  K profiles  so  derived  are  applicable  to 
stratospheric  sources,  concentrated  in  terms  of  latitude,  and 
with  tropospheric  sinks,  is  a matter  of  some  debate. 

The  procedures  used  by  Schmeltekopf  et  ai.  (1977)  give  K„ 
profile  data  for  roughly  the  17-27  km  region,  and  as  such  are 
not  alone  sufficient  for  use  in  1-D  modeling.  However,  Crutzen 
et  al . (1978),  from  the  same  group,  in  a study  of  chlorocarbon 
effects  on  ozone,  utilized  these  and  other  data  to  form  a K„- 
altitude  profile  which  they  characterized  as  follows  (in  cm^/sec) 
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^ Zj 

= 7 x 10  0 s;  z s 8 km 

Y.  r,  = -4  x 10'  + (4  x 10“  - 7 x 10~)(z  -15)/  3 ^ z ^ 15  km 

K„  = 4 x 10 J exp  [(z  - 15  )/5 . 54  ] 15  <.  z ^ 29  km 

4 

K,  = ; x 10  29  ^ z s.  50  kj. 

X . = 5 x 10  exp  [(z  - 50)/10.9J  z > 50  km 

The  values  from  15  to  30  km  are  based  on  the  N-0  data.  At 

a 

higher  altitudes,  the  values  are  based  on  methane  data  (Ehhalt 
et  al . , 1974),  adjusted  by  a factor  of  2 to  account  for  day- 

night  averaging.  As  noted  earlier,  the  Schmeltkopf  1577  K. 

z 

values  are  considerably  greater  in  the  15-25  km  region  than  are 


the  Hunten  1975  values. 

This  profile  was  included  in  comparison  studies  herein  cf 
effective  values  (see  Section  3-3-6  below)  and  is  also  of 
general  interest.  For  aircraft  (and  other)  effects  studies, 
in  a later  paper,  using  a 1-D  steady  state  model,  with  a 1C -km 
lower  boundary,  Crutzen  and  Howard  (1978)  used  a modified 
profile,  still  based  on  the  N?0  data,  as  follows: 

Above  15  km 

= 5 x 10"  exp  [(Z-15)/8j,  but  < 10 

From  10-15  km,  two  profiles  were  used 

K,,  (high)  = 5 x 10“  exp  [(15-Z)/2] 

K (low)  = 5 x 103 
*-* 

with  Z in  km  and  in  cnL/sec. 

The  statement  is  made  that  the  K (high)  relationship  is  more 
appropriate  to  global  conditions,  representing  a continuous 
transition  from  troposphere  to  stratosphere,  whereas  K„  (low) 
is  more  appropriate  at  latitudes  poleward  rf  40°.  These  rela- 
tionships are  sketched  on  Fig.  3-1,  along  with  certain  ether 
profiles  of  interest. 
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3.3.3  The  NASA-Ames  Model  K?  Profile 


Foppoff  et  al . (1973)  and  Turco  and  Whitten  (1977)  describe 
the  NASA-Ames  1-D  model.  The  K profile  used  by  them  (see  Sec- 
tion 3-5-5)  in  a recent  assessment  and  referred  to  as  a "WT"  >r 
"Wofsy-type"  profile  is  plotted  as  Poppoff  et  al.,  1973,  on 
Fig.  3-1 . This  profile  includes  a very  low  eddy  diffusivity 
region  near  the  model  tropopause,  followed  by  a rapid  increase 
in  K„  with  altitude.  The  plot  shows  a value  for  the  tropo- 
sphere taken  from  Turco  and  Whitten  (1977);  their  recent  model- 
ing exercises,  however,  have  not  included  the  troposphere 
(Poppoff  et  al.,  1978).  The  small  eddy  diffusivity  values  in 
the  13-15  km  region  lead  to  slow  transport  of  material  to  the 
troposphere  and  large  accumulations  (long  residence  time,  or 
burden/flux)  from  steady  stratospheric  sources.  The  profile  in 
the  region  was  developed  to  make  the  model  consistent  with 
transport  of  excess  C-lA,  following  Johnston  et  al  (1976).* 

The  Hunten/1975  profile  is  also  shown;  it  also  was  noted  by- 
Johnston  et  al.  (1976)  as  leading  to  a reasonable  representation 
of  the  excess  C-lA  data.  There  has,  of  course,  been  contro- 
versy surrounding  the  use  of  excess  C-14  data  in  this  manner. 

See  Chang,  1975;  Mahlman,  1975;  Oliver  et  al.,  1977  (pp.  3-35 
to  3-89);  also  Mahlman  and  Moxim,  1978  as  discussed  in  Section 
3.3.5  below. 

Turco  and  Whitten  (1977),  in  agreement  with  Johnston  et  al. 

(1976)  note  that  it  would  be  difficult  to  reproduce  excess  C-lA 
data  with  a tropopause  set  higher  than  13  km.  This  height  is 

< 

thus  properly  a global  model  tropopause  which  may  not  (in  our 
opinion)  have  any  physical  reality  relative  to  aircraft  flight 
in,  for  example,  mid-latitudes. 

1? 

The  WT  Kz  profile  also  leads  to  good  agreement  with  the 
measured  NpO  and  CFM  distributions  (R.C.  Whitten,  private 
communication,  1978). 
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3.3.4  Latitude-Tropopause  Adjustments 

Ir.  the  above  vein  and  as  noted  earlier,  the  proper  treatment 
of  aircraft  altitudes  relative  to  tropopause  altitudes,  and  pro- 
per mean  values  thereof  for  various  problems  of  interest,  has 
not  been  given  adequate  attention.  Hunten  (1975)  recommended 
a a -km  upward  adjustment  of  aircraft  altitudes  for  traffic  at 
about  A0°N,  to  compensate  for  tropopause  differences  between 
40°H  and  about  30°N,  for  which  his  model  K„  profile  was  devel- 
oped. Mahlman  (1975)  pointed  out  some  of  the  difficulties  with 
simple  approaches  of  this  type.  In  NAS  (1976)  Appendix  E,  some 
alternative  adjustment  procedures  were  recommended  for  tropo- 
pause height  variations;  if  adopted,  these  approaches  would  at 
least  keep  aircraft  positions  correctly  above  or  below  local 
tropopause  h<  ights.  (See  Oliver  e4-  al . , 1977,  pp . 3-94.)  The 
procedures  used  by  Crutzen  and  Howard  (1978)  were  described 
rlier  in  Section  1.2 . Popj  if 1 et  al.  (1978)  however,  based 
on  1-D  and  2-D  modeling  efforts,  argue  for  direct  use  of  the 
"WT"  profile,  with  a 13-km  tropopause,  without  aircraft  altitude 
adj ustment . 

3.3.5  Some  Relevant  Aspects  of  a 3-D  Tracer  Study 

Mahlman  and  Moxim  (1973)  have  carried  out  a model  experiment 
which  is  relevant  tc  the  questions  just  discussed.  In  this 
experiment,  using  a 3-D  global  general  circulation  model,  an 
instantaneous  pulse  of  Inert  water-soluble  tracer  was  injected 
on  1 January  at  65  mbar  (~  18.8  km)  with  its  center  of  3o°N  and 
l80°E.  The  pulse  was  spread  in  such  a way  that  most  of  this  I] 

tracer  was  contained  in  a "single  stratospheric  disturbance." 

One  aspect  of  their  results  is  of  particular  interest,  repre- 
sented by  Fig.  3-2.  In  this  figure,  time  and  zonally  averaged 
mixing  ratio  profiles  in  October  of  the  fourth  model  year  are 
shown.  While  It  is  not  possible  from  the  plot  alone  to  derive 
K profiles  (a  steady  state  distribution  does  not  prevail),  it 
is  of  interest  to  note  the  smoothness  of  the  "global  average" 
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FIGURE  3-2.  Tracer  study  results 

Source:  Mahlman  and  Moxim,  1978 
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curve  up  to  38  mbar  relative  to  curves  for  other  latitudes, 
general,  a break  in  the  curve  is  caused  by  a change  in  the 
rapidity  of  vertical  transport,  which  in  1-D  modeling  would  be 
represented  by  an  "eddy  diffusivity . "*  A decrease  in  slope  is 
analogous  to  a decrease  in  eddy  diffusivity.  Qualitatively, 
these  3-D  results  (in  our  opinion  at  least)  tend  to  support  the 
use  of  K profiles  with  relatively  smooth  transition  from  tro- 
pospheric  to  stratospheric  values,  at  least  for  "global  average" 
i-D  models. 

As  a separate  matter,  Mahlman  and  Moxim  (1978)  note  (refer- 
ring to  Mahlman,  1975,  for  more  detail)  that  feedback  of  C— 1 h 
from  the  troposphere  to  the  stratosphere  would  affect  apparent 
removal  rates  of  excess  C-14  during  the  second  year  following 
cessation  of  testing.  Thus,  the  relevance  of  excess  C-14  at 
longer  times  as  a tracer  is  brought  into  question  (excluding 
other  considerations);  however,  Mahlman  and  Moxim  (1978)  also 
find  that  the  water-soluble  gaseous  tracer  is  removed  somewhat 
more  slowly  than  were  particulates  from  weapons  tests,  indica- 
ting that  sedimentation  does  play  an  important  role  over  longer 
times  (perhaps  more  than  one  year). 

3.3.6  Studies  Using  Zr-95  Data 

In  a previous  report,  Oliver  et  al . , 1977  (Appendix  C)  work 
by  Bauer,  Oliver,  and  Wasylkiwskyj  was  reported  in  which  the 
behavior  of  Zr-95  injected  by  Chinese  nuclear  tests  (40°N,  90°E, 
approximately  18  km.  altitude)  was  examined  in  a 1-D  model  as  a 
function  of  time  after  several  events  at  different  seasons. 

This  work  has  since  been  substantially  revised  (Bauer  et  al., 
1978)  and,  for  completeness,  is  included  as  Appendix  B of  this 
•eport . 

if 

3-D  models  do  not  use  eddy  dif fusivities , except  in  subgrid 
parameterisations  in  which  very  low  values  are  used.  In  the 
particular  experiment  described,  zero  vertical  subgrid  eddy 
di ffusivities  were  used. 


The  reason  advanced  by  Bauer  et  al.,  1978,  for  examination 
of  the  Chinese  data  is  that  these  data  more  closely  resemble  a 
mid-latitude  stratospheric  aircraft  source  with  a tropospheric 
sink  than  do  any  of  the  other  data  used,  such  as  data  on  CH,,  or 
NpO  which  involve  a globally-dispersed  surface  source  and  a 
stratospheric  (primarily  high-altitude  tropical)  sink.  In 
Bauer's  approach,  the  available  data  are  treated  in  terms  of 
a two-step  model  (a  constant  tropospheric  and  a constant 
stratospheric  K ) which  is  justified,  at  least  for  an  initial 
investigation,  on  tne  basis  of  simplicity;  i.e.,  a minimum  ir. 
the  number  of  variables  is  needed  to  describe  the  data  involved. 
The  appropriate  form  for  mid- latitude  stratospheric  sources 
is  not  well  established  in  any  event.  A correction  is  included 
for  sedimentation.  Two  mean  tropopause  values  were  considered 
(11  and  14  km)  to  determine  the  sensitivity  of  results  to  this 
variable.  The  data  show  that  removal  rates  are  much  more  rapid 
in  the  winter  than  in  the  summer,  so  that  difficulty  Is  encoun- 
tered in  attempting  to  develop  an  appropriate  mean  K for  use 
in  1-D  time-  and  space-averaged  modeling.  By  making  what  were 
considered  to  be  suitable  assumptions,  however,  the  authors 
were  able  to  develop  approximate  mean  stratospheric  eddy  aiffu- 
sivities  (tropopause  to  18  km)  and  to  calculate  therefrom  mean 
residence  times  and  injection  coefficients  (a)  in  the  manner  of 
McElroy  et  al.  (197*0,  or  Hunten  (1975  or  1975a).  Results  are 
shown  in  Table  3“1.  Mean  K values  (K)  found  for  the  various 
tests  are  shown  in  the  tabulation.  The  recommended  values  are 
those  shown  with  an  assumed  11-km  tropopause. 

The  results  in  Table  3“1  suggest  that  residence  times  or 
injection  coefficients  for  18-km  injections  by  either  the  Chang/ 
1974  profile  or  the  Chang/1976  profile  are  somewhat  low,  but 
that  some  of  the  other  model  values,  particularly  using  the 
2-km  adjustment  on  the  injection  altitude  recommended  by  Hunten 
(1975),  are  apparently  high.  As  noted,  a strong  caveat  is 
implied  by  the  great  seasonal  variability  found  and  the  approxi- 
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TABLE  3-1.  COMPARISON  OF  RESULTS  FROM  DIFFERENT  MODELS  FOR  18-km  INJECTION* 


Ot , t-n  , 

- 17  2 K 

Source  10  cm  -sec  yr 

This  analysis  (Two-Step  Case) 


11-km  tropopause 


with  sedimentation  (K  = 7.25  x 103  cnr/sec) 

2.4 

2.0 

— 3 2 

without  sedimentation  (K  = 9.75  x 10  cm  /sec) 

1.8 

1.5 

14-km  tropopause 

with  sedimentation  (K  = 2.44  x 103  cm3/sec) 

5.0 

3.5 

— 3 2 

without  sedimentation  (K  = 4.55  x 10  cm  /sec) 

2.5 

1.9 

Other  K-Profiles* 

Chang/1974 

1.7 

1.4 

Chang/1976 

1.8 

1.6 

Hunten/1975 

4.6 

3.6 

Hunten/1974  (+2) (18  km,  latitude 

7.1 

4.6 

adjusted  to  20  km) 

Crutzen/1974 

1.8 

1.7 

Crutzen-Isaksen/1975 

1.9 

1.7 

Wofsy/1975 

3.2 

2.6 

Crutzen  et  al . , 1978 

2.3 

2.0 

Poppoff  et  al . , 1978** 

4.4 

3.6 

Crutzen  and  Howard  1977/1978 

High 

2.4 

2.1 

Low 

3.3 

3.2 

★ 

See  Appendix  B (or  Bauer  et  al . , 1978)  for  details. 

**  4 

A Kz  of  7 x 10  in  the  0-10  km  region  was  assigned  for  these  calculations. 
The  troposphere  is  not  included  in  the  Poppoff  et  al.,  1978  work. 
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mations  needed  to  get  annual  means.  It  is  suggested,  however, 
that  since  the  large  Soviet  tests  in  1961-1962  took  place  in  the 
late  fall,  overestimates  of  pollutant  burdens  may  result  if  a 
year-round  average  value  is  applied  in  computing  effects. 

Further  analysis  of  these  and  other  data  using  a more  com- 
plex model  (K  increasing  with  altitude  above  the  tropopause)  is 
anticipated  in  follow-on  work. 

3.3.7  The  Water  Problem 

a.  Stratospheric  Water.  The  original  focus  of  potential 
SST  effects  on  ozone  was  on  effects  due  to  added  water  vapor 
(see,  e.g.,  H.  Harrison,  1970;  Hearings,  1971,  p.  451  ff).  A 
change  of  0.2  ppm,  as  was  suggested  might  result  from  an  SST 
fleet,  was  estimated  to  cause  a global  average  reduction  in 
ozone  of  about  0.8  percent  based  on  calculations  by  London  and 
Park  (Hearings,  1971).  The  reaction  rates  used  were  subsequently 
revised,  with  water  vapor  becoming  of  much  less  importance.* 

The  N0x  issue  was  raised  in  1971  (Crutzen,  1971;  Johnston,  1971) 
and  became  of  dominant  interest.  Lately,  however,  with  further 
revisions  in  reaction  rates  and  with  the  N0x  problem  becoming 
of  lesser  significance,  the  water  problem  has  arisen  again,  as 
will  be  discussed  more  fully  in  a later  section.  This  leads  to 
the  question  of  how  water  vapor  should  be  handled  in  1-D  and 
2-D  models. 

* 

These  early  calculations  were  based  on  an  assumed  value  for 
the  rate  of  the  HO2  + Oj  reaction  of  10"^  cm3  molec-Isec~3 , 
as  used  by  Hunt  (1966)  in  modeling  the  ozone  column  using  an 
HOx  - 0X  atmosphere.  It  was  later  pointed  out  by  Nicolet 
(1970)  that  the  rate  was  probably  much  too  high;  modelers 
accepted  his  arguments  and  (initially)  set  the  rate  to  zero 
in  investigating  the  N0X  catalytic  cycle  and  Its  effects  on 
ozone  (Crutzen,  1971,  197*0.  It  is  of  interest  that  the 
early  water  vapor  efforts,  as  well  as  later  efforts  involving 
N0X,  illustrate  the  fact  that  a model  which  reasonably  dupli- 
cates a "natural"  atmosphere  may  not  correctly  predict  the 
effects  of  a perturbation  to  that  atmosphere.  Duplication 
of  the  "natural"  atmosphere  in  a model  is  generally  considered 
to  be  a necessary,  but  clearly  not  a sufficient  condition,  to 
predict  effects  of  perturbations. 
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There  are  several  related  aspects  to  the  water-vapor  question 
First,  perhaps,  is  the  limited  data  base  in  terms  of  reliable 
measurements  at  various  times  and  places.  (See  Harries  1976, 
Pendorff  1977 > Ellsaesser  197^- ) According  to  Ellsaesser  (137*0* 
the  data  that  do  exist  suggest  that  the  conventional  tropical 
tropopause  cold-trap  argument  (Erewer,  19*0)  is  quantitatively 
unsatisfactory  in  explaining  stratospheric  water  vapor;  a sink, 
specific  for  Ho0  vapor,  is  needed,  which  Ellsaesser  suggested 
may  exist  in  the  Antarctic  night.  More  data  will  be  required 
(temperatures,  upwelling  velocities,  relative  humidities)  to 
fully  test  these  suggestions.  It  seems  possible  that  the  appar- 
ent problem  may,  in  fact,  be  an  artifact  of  averaging  processes.* 
The  vapor  pressure  of  water  (ice)  varies  by  a factor  of  about  2 
with  a A°C  change  in  temperature  near  -80°C,  and  regions  in  the 
tropics  exist  (as  pointed  out  by  Ellsaesser,  citing  Angell  and 
Korshover  197*0  wherein  12-month  running  mean  temperatures  are 
-8a.5°C  to  -81. 5°C , significantly  colder  than  the  -80°C  nominal 
tropical  tropopause  used  by  Ellsaesser  in  his  analysis.  It  may 
be  precisely  these  cold  regions  through  which  air  is  entering 
the  stratosphere.  Similarly,  tropospheric  air  enters  the  strat- 
osphere on  the  anticyclonic  side  of  jet  streams  in  mid-latitudes; 
in  these  regions,  temperatures  considerably  below  zonal  mean 
values  are  not  uncommon  (J.D.  Mahlman,  NOAA  GFDL,  private  com- 
munication, 1978).  These  zonal  inhomogeneities  can  be  (and  are) 
modeled  in  existing  3-D  models  (as  at  GFDL) , resulting  in  a 
satisfactorily  dry  stratosphere.**  (A  CHj,  contribution  is  not 
included  in  the  GFDL  model.)  None  of  this,  of  course,  proves 
¥ 

I am  indebted  to  Dr.  J.D.  Mahlman  of  NOAA/GFDL  for  pointing 
out  this  and  other  "3-D"  features  of  the  water  vapor  question. 

It  is  of  interest  that  the  GFDL  model  referred  to  has  an 
excessively  cold  polar  night,  resulting  in  an  excessive 
estimate  of  the  Antarctic  sink  proposed  by  Ellsaesser. 

Even  so,  the  resulting  "sink"  is  of  almost  no  signifi- 
cance. (J.D.  Mahlman,  private  communication,  1978.) 


- 
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that  the  node Is  duplicate  the  real  atmosphere,  as  the  data 
against  which  to  test  the  models  is  inadequate. 


The  above  discussion  suggests  that  water  vapor  transport 
into  (and  out  of)  the  stratosphere  is  dependent  on  processes 
which  can  be  represented  in  3-D  modeling,  but  which  lead  to 
difficulties  in  1-D  or  2 -D  models.  In  these  simpler  models, 
the  cold-trap  concept  cannot  be  applied  if  mean  temperatures 
are  employed;  rather,  some  further  prescription  ' s necessary. 

This  is  conventionally  done  in  1-D  models  by  specifying  a water 
vapor  content  above  (or  at)  some  model  tropopause,  and  in  some- 
what more  complex  ways  in  2-D  models.  It  is  not  overly  diffi- 
cult to  develop  some  prescription  in  these  models  which  leads 
to  a reasonable  "natural"  stratosphere,  but  it  is  difficult  to 
establish  that  this  prescription  will  permit  an  estimate  of  the 
perturbations  to  be  expected  from  a stratospheric  water  source. 
This  is  certainly  true  if  thermal  feedback  effects  of  water 
vapor,  or  effects  on  the  tropical  tropopause  temperature  (with 
changes  in  water  vapor)  are  in  question. 

b . Water  Vapor  in  the  Upper  Troposphere--Imp1 i c a t i o n s . 

That  the  stratosphere  is  dry  (very  low  relative  humidity,  except, 
presumably,  near  the  tropical  tropopause)  in  certain  regions  is 
well  established.*  There  seems  to  have  been  less  attention  given 
to  upper  tropospheric  water  contents  in  mid-latitudes,  but  the 
information  that  is  available  may  have  some  implications  for 
1-D  and  2-D  modeling  problems.  Thus,  data  exist  which  show  that 
upper  tropospheric  air,  in  general,  apparently  being  mixed  with 

dry  stratospheric  air  (see  Danielsen  and  Louis,  1977),  is  also 

— — - 

Manabe  and  Wetherald  (1967),  after  reviewing  rather  limited 
data,  modeled  relative  humidity  h as  a function  of  Q(e  ambient 

pressure/surface  pressure)  according  to  h = h*  y~-r'4v  , where 

h*  = 0.77,  with  adjustment  to  prevent  the  mass  mixing  ratio 
from  dropping  below  3 ppmm.  Manabe  and  Wetherald  (1975)  show 
model  (2-D)  relative  humidity  values  (zonal  means)  of  140-50 
percent  in  large  areas,  including  the  tropical  tropopause, 
which  figures  presumably  reflect  the  averaging  processes 
described  in  Section  3-3-7a  above. 
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fairly  dry  (i.e.,  low  relative  humidity).*  The  Amarillo  data 
plotted  on  Fig.  3.3  at  about  30oN,  show  a relative  humidity 
about  19  percent  at  the  tropopause.*  It  is  of  interest  alsc 
that  the  Amarillo  data,  as  well  as  data  at  Washington,  D.C., 
and  Trinidad  (CIAP  Monograph  1,  p.  3-52)  show  a gradual  champ 
in  mixing  ratio  across  the  tropopause;  however,  other  data  (sane 
reference)  over  Greenland  show  abrupt  increases  below  about  250 
near  (which  might  represent  the  tropopause).**  With  some  excep- 
tions, as  in  the  cold  regions  near  jet  streams  referred  to 
earlier,  and  data  such  as  that  for  Greenland,  the  air  near  the 
tropopause  in  mid-latitudes  is  not  saturated.  The  existence, 
on  average,  of  rather  dry  upper  tropospheric  air  suggests  that 
vertical  mixing  from  below  in  this  region  is  not  (on  the  average) 
extremely  rapid,  as  implied  by  models  using  large  values  of 
tropospheric  eddy  dif f usivity : further,  the  lack  of  gradient 
across  the  tropopause  suggests  that  extreme  step  changes  in 
eddy  diffusivity  at  the  tropopause  are  not  representative  of  the 
actual  atmosphere.  There  is  some  additional  support  for  this 
lack  of  step  change  in  the  data  on  HTO  distributions  with  alti- 
tude across  the  tropopause  (see  Mason  and  Ostluna,  1976)  wherein 

abrupt  changes  in  concentration  or  in  gradient  were  not  found  in 

+ 

crossing  the  tropopause. 

c.  Mater  Vapor  - Genera).  The  somewhat  murky  picture  (in 
a quantitative  sense)  of  water  vapor  transport  in  the  upper 

if 

Saturation  mixing  ratios  in  Fig.  3"3  have  been  computed  from 
the  1966  U.S.  Standard  Atmosphere  and  from  the  Smithsonian 
Meteorological  Tables  (List,  1958). 

The  data  of  McClatchey  et  al.  (1972)  as  plotted  by  Ramanathar. 
(1977)  alsc  indicate  low  relative  humidities  in  the  upper 
troposphere  in  the  Arctic  (winter  and  summer),  mid-latitudes 
(winter  and  summer;,  and  the  tropics  (~  15  percent  at  10  km). 

The  plotted  data  appear  to  be  suspect,  however,  in  the  stra- 
tosphere (~  0.6  ppmm) . 

4. 

‘It  is  of  interest  that  Poppoff  et  al.  (1978)  view  this  same 
data  in  a directly  opposite  way.  We  plan  to  examine  these 
data  more  carefully  in  future  work. 


_ 
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troposphere  and  stratosphere  is  (following  Dan' els-  . 

1977  , who  do  not  discuss  water  but  do  discus:-  tr  \ 
the  principal  source  of  stratospheric  air  is  f rcr.  ,<pwe.. 
certain  regions  of  the  tropics  with  little  \ ,-r4  ica i tr-: 
elsewhere  and  general  subsidence  poleward.  Th<-_-  air 
the  stratosphere  carries  with  it  a small  arrmuir 
which,  in  general,  is  then  transported  by  bulk:  *"1  w 
so-called  tropopause  gaps  into  the  tr  ; is:  here,  1 

further  mixed  longitudinally  and  meridionally.  S - 
obviously  also  is  carried  into  the  si  ratosr.  her*.  b 
from  mid-latitudes,  as  evidenced  by  regions  of  r 
water  vapor  mixing  ratios  (6  ppmm  or  10  ppm:;  ...  . 
traverse  (80°N  to  50°S  at  17-19  km)  by  Hilsenratr 
Antarctic  sinks  may  exist  (Ellsaesser,  197--;  dr.anf- 
Kuhn,  1977),  as  may  Arctic  sinks,  but  these 
in  terms  of  the  water  budget  and  transpor- . It  w . 
water  vapor  introduced  above  the  tropopause  In  :-- 
has  an  uncertain  fate;  some  meridional  transj 
the  sloping  tropopause  would  certainly  occur.  Becau: 
general  subsidence,  vertical  transport  would  be  so  . 
r.ot  eliminated. 

For  reasons  such  as  these,  and  because  was  i 

ent  consequence,  CIAP  modelers  generally  made  - 

include  water  vapor  injections  (COMESA  modelers  did  . 
of  interest  that  the  K_  profile  generated  by  COMESA 
vapor  ’’behaves"  better  than  that  for  NO  (see  - : n. 
latter  involves  negative  regions  near  18  km.  How- v 
COMESA  model  shows  much  too  high  a water  content  in 
sphere  (~  50  ppmm  at  20  km  in  the  Northern  Hemisphere, 
ing  with  increasing  altitude,  and  varying  with  seaso:  . 
of  the  order  of  5 ppmm  at  33  km  in  the  fall;  see  Fin'.  - 
COMESA,  1975),  so  that  it  is  not  clear  tc  what  ext  :n  ' * 
rived  K„  profile  for  H,0  duplicates  the  "true"  atm  r- 

duced"  to  one  dimension. 
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Kjim2  sec) 

FIGURE  3-4.  Globally  averaged  one-dimensional 
Kz  profile  based  on  water,  ozone 
and  NO2  fluxes  and  gradients;  derived 
from  3-D  results. 

Source:  COMESA  1975 

Ir.  some  cases  in  CIAP  and  in  HAPP,  estimates  of  the  effects 
f water  vapor  have  beer,  made  using  adjusted  water  contents 
based  or.  residence  time  or  other  arguments  (as  at  Lawrence 
Livermore  Laboratory;  see  Oliver,  1977);  in  other  cases,  rain- 
out  c efficients  (essentially  fractional  rates  of  decay  at  some 
model  altitude)  have  been  used.  In  still  other  cases,  modelers 
have  fixed  the  water  content  at  some  level  (as  NASA  Ames, 

3 x 10'-’  mol/cm*'  at  10  km;  Turco  and  Whitten,  1977),  and  assumed 
the  same  L ..  profile  applies  to  water  as  to  other  species.  "Rain- 
out  coefficients"  have  also  been  applied  in  2-D  models  over  cer- 
tain altitude  and  latitude  bands  in  the  troposphere  (G.  Widhopf, 
private  communication,  1978). 

Clearly,  as  water  vapor  effects  on  ozone  and  thermal  feed- 
back effects  become  of  significance,  more  debate  is  called  for 
on  the  "proper"  way  to  model  water  vapor  in  reduced  dimension 
mode  Is . 
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3.4  SOME  IMPORTANT  MEASUREMENTS  OF  TRACE  SPECIES 


A considerable  number  of  field  measurements  of  trace  species 
has  taken  place  or  been  reported  in  the  period  between  this  and 
our  previous  report.  Two  of  these  efforts  seem  particularly 
worthy  of  note  in  view  of  the  questions  raised  by  the  measure- 
ment s . 

The  first  of  these  (the  Cl  and  CIO  measurements  of  Anderson 
et  al.,  1977,  and  Anderson,  1978)  is  of  particular  importance 
to  the  halocarbon  problem,  but  is  significant  to  the  general 
understanding  of  ozone  chemistry  and,  because  background  C1X 
content  is  important  to  the  computations,  to  the  specific  ques- 
tion of  the  response  of  stratospheric  ozone  to  perturbations  by 
aircraft  effluents.  In  brief,  Anderson  et  al.  (1977)  have  re- 
ported three  sets  of  profile  measurements  of  Cl  and  CIO  taken 
28  July  1976,  2 October  1976,  and  8 December  1976  over  Palestine 
Texas,  at  altitudes  from  roughly  25  to  40  km.  A 14  July  1977 
flight  (Anderson,  oral  presentation,  Toronto,  1978)  also  took 
place.  The  Cl  measurements  involved  no  surprises:  the  CIO 
measurements,  however,  showed  wide  variability  among  the  dif- 
ferent flights — a factor  of  about  10.  Furthermore,  three  of 
four  profiles  show  a peak  CIO  content  higher  than  called  for  by 
any  of  the  models;  the  profile  for  8 De  :emb<  1 . , in  ; art  ! • ila 

>h  wed  ~ S ppb  of  CIO  at  30  km,  some  three  t imes  the  highest 
value  predicted.  As  at  least  some  of  the  C1X  is  nec< ssarily  in 
the  form  of  HC1,  the  total  C1X  content  of  the  stratospheric 
air  implied  by  the  flight  was  clearly  far  larger  than  can  be 
accounted  for  by  any  known  sources  of  chlorine.  These  results 
have  not  yet  been  explained;  a discussion  is  included  in  Hudson, 
1977.  The  July  1977  flight  also  showed  excessive  chlorine 

7 ppbv  It  is  of  interest  that  in  this  flight,  the  OH  and 
Q,  contents  looked  normal  (Anderson,  r;  1 : resentation,  Toronto, 
L97  . ; 1 the  excess  chlorine,  if  present,  was  havinj 
vi  lent  Impact  on  ozone.  Wofsy  (1976  ),  in  a detailed  review  of 
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nd  experiment  with  regard  to  stratospheric  trae<  ;ases , 

: r:  :tes  an  excess  of  CIO,  and  notes  difficulties  in  recon- 

H,  '10,  and  0^  observations  in  the  upper  stratosphere. 

finds  the  theory  to  explain  adequately  the  major  features  of 

. • . de  and  seasonal  distributions  of  NO  and  UNO.,.  It  should 

x 3 

:led  that  important  rate-constant  revisions  have  been  made 

ati<  s wer<  :arried  ut , s t h - 1 t he  re s tilts  — 
..it  :h  r all  of  the  work  in  this  field--may  need  reexam- 


noted  above,  C1X  content  of  the  stratosphere  affects 
■ie  results  for  NO  injections  with  increased  chlorine,  leading 
! n ed  or  more  positive  effects.  This  has  been  shown  by 
. . t . in  is  at  Lawrence  Livermore  Laboratory,  with  their  results 
Lvei  el  al.,  1977  ; no  si  mi 1 a v results  with  the 
■ mistry  are  available. 

’ r d measurement  program  relates  to  the  NC  . content  of 

A 

. r-e.  The  measurements  involved  are  those  by  Noxon 
, ...  oer-rts  tropospheric  N0o  (no  NO  ) measurements,  and 

r :.  spherin'  NO  tc  be  perhaps  ar.  ordc r-«. f-magnitude 
ighl  previ  us3  y . Th<  results  art  Lmj  rl  anl 
ions  nd  1 hus  t 'H  , K,  Cl,  ther  hydro- 

*4 

iarbons , i 11 feti mi  s , and  thereby  t strate- 


gist ent 


Rest  j.nce  of  trorosphei 


;u  ! i ma;  als  lei  nd  n the  abs 


the  tropospher 
■;  winici  nan  •. 


ted,  h we n i , 


s i t , i i t : ■ L c a t i oi  a 1 Least, 

■ ' n . , has  been  criticized  ns  giving  erratic  re- 

T.netsov  and  7!  — aiuliina  (1977'-  Tl.e  implicat  i ns 
I 1 ialgc  (1978),  but  no  new  model! 
aircraft  NO  injections  into  the  troposphere  are 
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3.5  CHEMISTRY  CHANGES  AND  MODELING  RESULTS 


3.5.1  Introduction 

In  the  sections  just  preceding,  certain  developments  have 
been  noted  on  the  parameterization  of  transport  processes.  These 
developments  have  not,  in  general,  been  fully  "folded"  into  the 
modeling  results,  partly  because  ambiguities  are  involved,  but 
also  and  primarily,  because  other  developments — in  chemistry-- 
have  occurred  which  have  been  of  more  immediate  significance. 

These  developments  in  chemistry  have  reduced  net  effects  of  NO 

x 

injection  on  the  ozone  column  substantially  while  somewhat  in- 
creasing H20  effects;  whether  the  net  effect  on  the  ozone  column 
is  positive  or  negative  depends  on  a number  of  modeling  variables, 
included  among  which,  however,  is  the  parameterization  used  for 
transport.  Transport  processes  are  thus  still  of  significance. 

3.5.2  Chemistry  Changes  and  Uncertainties 

a.  Previously  Identified  Problem  Reactions.  In  a pio- 
neering study  of  uncertainties,  Duewer  et  al.  (1977  and  1977a) 
pointed  out  the  importance  of  five  key  reactions  in  evaluating 
perturbations  of  stratospheric  ozone  by  aircraft  (in  a chlorine- 
free  atmosphere):  These  are,  using  their  reaction  numbers: 


HO-  + NO  -*■  OH  + NO-  ( R 2 6 ) 

d.  2 

HO-  + OH  ■+  02  + K-0  (R19 ) 

HO-  + HO ^ - H ^ 0 ^ +'o2  (Rl8) 

HO-  +0-+0H+20-  (RlU) 

OH  + N02  “ HNO-  ( R20 ) 


These  reactions  have  since  been  reviewed  by  the  same  group 
(Duewer  et  al.,  in  Luther,  1977)-  Recommended  values  for  these 
(and  numerous  other  reactions)  have  also  been  provided  in  Hudson 
(1977),  along  with  discussion  and  estimated  uncertainties. 
Crutzen  and  Howard  (1978),  and  Poppoff  et  al.  (1978)  also  dis- 
cuss the  various  rate  data  for  these  and  other  reactions  of 
significance . 


The  rate  of  reactions  involving  H0o,  particularly  at  strato- 
spheric temperatures  and  pressures,  are,  in  general,  still  not 
fully  established,  although  important  progress  has  been  made. 
Eased  on  at  least  one  example  (see  below)  the  uncertainties 
quoted,  in  our  opinion,  tend  to  be  underestimated.  In  fact,  the 
most  important  changes  in  computed  aircraft  effects  since  our 
previous  report  have  resulted  from  the  measurement  at  room  tem- 
perature by  Howard  and  Evenson  ( 1 S 7 7 ) of  the  rate  of  Reaction  26: 


H0„  + NO  - OH  + N0„ 

2 


(R26) 


- 12  1 

The  race  found,  (8.1  ± 1.5)  * 10  cm  /sec  at  296  K,  has  since 

_1P 

been  confirmed  (±  0.1  x 10  “ ) by  four  groups  (Howard,  1978). 

This  figure  can  be  compared  to  the  "survey"  value  quoted  in 
CI.'.P  Monograph  1 (p.  5-227)  at  300  K of  2 x 10-°  cirr/sec  (with 
an  uncertainty  factor  quoted  of  3 either  way),  a factor  of 
40  slower. 

Ir  very  recent  work  (reported  June  1978),  Howard  (1978) 
reports  the  following  temperature-dei  indent  rate  for  this  re- 


. ion  as 


= (3. 


) x 10  ^ exp  [(25^  ± 50)/T] 


over  the  range  230-000  K.  This  positive  temperature  dependence, 
and  the  rate  measured,  are  in  marked  contrast  to  earlier  values, 
.'ivonaitis  and  Heickler.  (1977  ) reported,  e.g.,  a rate  constant 
at  room  temperature  of  1 x 10~  ^ cm~Vsec  and  a temperature  de- 
pendence, based  on  certain  assumptions  as  to  other  reactions 
-f  . - 1400  ± 500)/T.  Hudson  (1977)  gives  no  temperature  de- 
reud-nce,  but  gives  an  uncertainty  at  230  K,  corresponding  to 
a •cultiplier  of  1.4  to  0.35. 

In  studies  prior  to  the  Howard  (1978)  work,  authors  assumed 
t : .perature  dependence  as  follows,  with  k values  shown  at  room 
t-vperature  and  at  a nominal  stratospheric  temperature: 


\ 


i 


^ • 


Crutzen  and  Howard,  1978 
k26(HE)  = 3 x I0""1  exp  (-390/T) 


K296  K k230  K 

3.0  x 10"1"  5-5  x 10 


Luther  et  al.,  1977  (p.  47) 

k26(rIE)  = 4,28  x 10-liexp(-500/T)  7-9  x 10-^"  4.9  x 1C-12 


In  later  work  at  Lawrence  Livermore  Laboratory  (LLL)  (Wuebbles, 
1978),  the  above  temperature  dependence  for  this  reaction  was 
dropped  (See  Section  3.5.4c.). 

The  central  value  expression  for  k„g  developed  by  Howard 
(1978)  yields  9-96  x 10~12  at  230  K.  The  much  higher  than 
heretofore  measured  rate  for  this  reaction,  its  positive  tem- 
perature coefficient,  and  the  fact  that  measurements  by  the 
technique  recently  used  by  Howard  (and  others)  require  low 
pressures  (<  10  Torr,  Burrows  et  al.,  1978)  have  raised  the 
possibility  that  a pressure-dependent  three-body  effect  may  be 
involved.  Howard  (1978)  argues  however,  based  on  certain  re- 
lated data,  that  the  measured  rate  as  reported  is  applicable 
at  atmospheric  pressure. 

This  reaction,  with  this  new  rate,  has  a powerful  influence 
on  both  tropospheric  and  stratospheric  ozone  chemistry  in  at- 
mospheres containing  NO  , as  shall  be  shown. 

A second  important  reaction,  which  accounted  for  much  of 
the  change  in  computed  NO  effects  on  ozone  between  the  period 

X 

of  the  CIAP  report  (Grobecker  et  al.,  1974)  and  our  recent  re- 
port (Oliver  et  al.,  1977  ),  is  Reaction  19'. 


H02  + OH  - 02  + Ho0  (R19 ) 

The  CIAP  figure  of  2 x 10  cm8/sec  was  revised  in  the 
previous  Duewer  et  al.  evaluation  (1977  ) to  2 x 10-'*'2  cm'Vsec; 

the  Hudson  report  recommends  3 x 10-'*'2.  Burrows  et  al.  (1977  ) 

■“11 

reported  a direct  measurement  of  5.1  x 10~  , but  this  figure, 
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based  on  other  evidence,  is  considered  to  be  high  both  by  LLL 
(Duewer  et  al.,  1977;  Luther,  1977)  and  by  the  NASA  panel.  (One 
might  argue,  in  general,  for  the  validity  of  direct  measurements 
over  indirect,  in  view  of  experience  with  the  HO^  + NO  rate  data.!' 

The  reaction 

H0„  + H0o  -+  H~0„  + 0o  (Rl8 ) 

d d d c. i d 

is  also  clearly  of  significance  in  view  of  the  importance  of 
HC^;  e.g.,  in  its  reaction  with  NO.  The  importance  of  the  re- 
action also  depends  on  the  rate  of  destruction  of  the  ^oriT!ec*> 

as  discussed  by  Pcppoff  et  al.  (1978).  Several  room  temperature 
measurements  of  the  rate  have  been  made  (Luther,  1977),  which 
appear  to  be  in  good  agreement  (at  2.5  x 10  “ with  ± log  k of 

0.3;  Hudson,  1977).  However,  there  is  concern  (Hudson,  1977) 
in  light  of  the  Howard  and  Evenson  (1977)  measurement,  that 
competing  reactions  may  have  been  inadequately  accounted  for 
(Burrows  et  al . , 1978).  An  argument  is  advanced  by  these  authors 
for  the  formation  of  Ho0^  as  an  intermediate  and  extended  to 
argue  for  possible  pressure  dependence  in  HO^  + NO  and  CIO  + NO 
reaction.  No  temperature  dependence  data  exists  for  Reaction  18; 

+1000 

the  Luther  (1977)  report  suggests  a range  of  e — ^ — which  alone 
represents  a much  larger  uncertainty  than  indicated  in  Hudson 
(1977). 

Crutzen  and  Howard  (1978)  in  parametric  studies  used  two 
values  for  this  reaction: 

k (standard)  = 5 x 10~'lc  e~500/.  (9-3  x lO-^  at  298  K) 

k (high)  = 1.5  x 10-11  e_500/T  (2.8  x 10-1^  at  298  K) 

A fast  rate  constant  for  this  reaction,  according  to  Crutzen  and 
Howard  (1978)  and  Duewer  et  al.  (1977),  increases  ozone  deple- 
tion. 

The  rate  of  the  reaction 

H0~  + 0,  OH  + 2 0„  ( R1 A ) 

c 3 2 
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is  also  (still*)  quite  significant  to  modeling  of  ozone  in  nat- 
ural and  perturbed  atmospheres.  A very  recent  (June  1978)  ex- 
pression has  been  reported  by  Zahniser  and  Howard  (1978)  as 


kl4  = (1.4  ± 0.4)  x 10' 


exp  (-580  ± 100 ) /T 


over  the  range  250  - 370  K.  At  230  K,  the  central  value  using 

-15  3 

this  expression  is  1.1  x 10  cm  /sec . 

The  rate  expression  used  by  LLL  for  Reaction  14  is  1C 

exp(-1250/T) , which  yields  4.4  x 10_J'J  cm^/sec  at  230  K.  The 

3 -14 

rate  quoted  in  Hudson,  1977  is,  in  cm  /sec,  k^4  - 7.3  x 10 

e 127 j/T^  wj_£h  + log  k equal  to  0.3  (a  factor  of  2).  The  ex- 

-16 

pression  yields  2.9  x 10  at  230  K.  Crutzen  and  Howard  (1978) 
use  two  values,  a "standard"  and  a "high"  value  as  follows 

k standard  = 10-1^  exp  (-1525/T),  or  1.3  x 10-1^  at  230  K 

and  _-.fr 

k high  = 2 x 10  3 exp  (-1250/T),  or  8.7  x 10  at  230  K 

The  NASA  (Hudson,  1977)  rate  is  sensitive  to  the  rate  of  HC^  + 

HO.,,  which,  as  pointed  out  earlier,  is  uncertain,  and  may  be 
influenced  in  turn  by  reevaluation  of  the  impact  of  the  HO^  + 

NO  reaction.  Note  that  the  new  rate  expression  given,  at  230  K, 
a value  3.8  times  that  given  by  the  expression  in  Hudson  (1977). 
The  uncertainty  quoted  in  Hudson  (1977)  corresponded  to  a fac- 
tor of  2. 

Crutzen  and  Howard  (1978)  [as  well  as  Poppoff  et  al.  (1978), 
also  Duewer  et  al.  (1977  and  1977a),  Table  3.5,  later  herein] 
show  that  if  the  reaction,  H00  + 0-,  is  "fast,"  the  "natural" 
ozone  column  is  reduced,  but  NO  additions  tend  to  enhance 
ozone.  The  effect  is  a subtle  one,  involving  a decrease  in  a 
reduction  (which  thus  becomes  an  enhancement);  i.e.,  as  ex- 
plained by  Crutzen  (private  communication,  August  1977).  If 
Cycle  1,  for  which  the  H0„  + 0.^  rate  is  controlling, 


See  first  footnote  Section  3.3.7a. 
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OH  + C>3  -*■  H02  + 02 
HO„  + 03  - OH  + 2 02 


is  rapid,  then  this  ozone  destruction  cycle  is  significant.  As 
NO  is  added,  however,  the  cycle  is  slowed  by  the  competitive 
reaction  of  H02  with  NO  as  follows: 

OH  + 0^  H02  + 02 

H02  + NO  ->  N02  + OH 
N02  + hv  -*•  NO  + 0 


No  net  reaction 


which  decreases  the  destructiveness  of  the  HO  cycle. 

The  fifth  reaction  listed  by  Duewer  et  al.— that  forming 
HNO^,--still  involves  a number  of  uncertainties  but  other  reac- 
tions, as  below,  are  of  greater  current  interest. 

Hudson  (1977)  cites  two  additional  reactions  in  the  HO 
series  which  the  NASA  panel  felt  to  be  worthy  of  particular 
attention , 

OH  + 0o  H0„  + 0o 

and  3 2 2 

H02  + 0 -+  OH  + 0? 


The  uncertainties  in  these  rates,  in  addition  to  the  first 
four  cited  above,  contribute  significantly,  according  to  Hudson 
(1977),  to  the  overall  uncertainty  in  sti atospheric  models. 

b . Some  Other  Reactions  and  Species  of  New  or  Possibly 
Increased  Significance.  There  has  been  considerable  effort 
in  the  past  few  years  devoted  to  stratospheric  chemistry, 
which  is  particularly  important  to  the  halocarbon  problem  (see, 


e.g.,  Hudson,  1977;  Watson,  1977)-  As  most  of  this  work  is  of 
secondary  interest  to  the  aircraft  problem,  no  attempt  will  be 
made  to  review  it  here.  Rather,  discussion  will  be  limited  to 
a few  items  which  seem  to  be  of  interest. 

HOgNOg  (peroxynitvic  acid).  HO^NO^  was  not  considered  in 
modeling  done  under  CIAP  nor  has  it  been  included  in  most  reac- 
tion schemes  since.  Attention  to  this  species  was  apparently 
drawn  by  Simonaitis  and  Heicklen  (1976)  with  subsequent  atten- 
tion being  given  it  by  Cox  et  al.  (1977)  and  by  Jesson  et  al. 
(1977)  • 

The  species  is  an  analog  of  the  well-known  peroxyacetyl 
nitrate  (PAN)  important  in  smog  chemistry.  According  to  Cox 
et  al.  (1977),  HC^NC^  is  not  important  in  the  troposphere, 
having  a lifetime  of  only  9 seconds  at  298  K.  At  stratospheric 
temperature  (220  K),  however,  the  lifetime  is  estimated  at  about 
one  month  (Jesson  et  al.,  1977)  to  0.3  years  (Cox  et  al.,  1977). 
It  could  thus  be  an  important  "reservoir"  for  both  H02  and  NO^, 
and  could  enter  into  a variety  of  possibly  significant  reactions 
(Jesson  et  al . , 1977);  its  peak  concentrations  would  be  in  the 
20-30  km  region,  and  could  reach  1-3  ppb,  largely  at  the  expense 
of  HNO^  concentrations  (Jesson  et  al . , 1977)  without  NO  and  NO^ 
concentrations  being  much  affected. 

Absorption  data  and  photodissociation  products  are  appar- 
ently inadequate.  Modeling  runs  showing  the  possible  effect  of 
this  species  are  not  extensive;  however,  Poppoff  et  al.,  1978, 
assuming  H02N02  to  be  photolyzed  at  the  same  rate  as  H20o , 
fjund  that  its  inclusion  reduced  ozone  production  by  NO^,  but 
the  absolute  effect  on  the  ozone  column  was  small. 

HOCl.  Jaffe  and  Langhoff  (1978)  have  carried  out  a theo- 
retical study  of  the  ultraviolet  photodissociation  of  HOCl. 

The  results  do  not  agree  with  earlier  experiments  which  indi- 
cated a double-peaked  absorption  spectrum  with  peaks  at  about 
220  and  320  nm;  the  theoretical  results  showed  significant 
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absorption  only  at  the  shorter  wave  length.  The  significance 
of  this  is,  of  course,  the  H0C1  is  not  photodissoc iated  by  the 
relatively  strong  flux  in  the  300-350  nm  region,  so  that  H0C1 
has  a much  longer  life  than  previously  believed  for  Cl  species, 
reducing  the  effectiveness  of  the  Cl  catalysis  cycle  (see  Prasad 
et  al.,  1978).  Jaffee  and  Langhoff  (1978)  estimate  HOC1  life 
due  to  photodissociation  (to  Cl  + OH)  as  5-3  h at  ^0  km  and  5 8 
h at  30  km.  They  point  out,  however,  that  other  processes  as, 
e.g.,  reaction  with  0(  P)  may  predominate  in  controlling  the 
lifetime . 

H0C1  could  have  some  significance  in  the  aircraft  effects 
problem,  thus  such  processes  as 

H02  + CIO  -*  H0C1  + 02 
or  (heterogenously ) 

C10N02  + H20  -*■  H0C1  + UNO,  . 

(These  formation  mechanisms  were  suggested  by  Jaffee  and  Langhoff, 
1978)  . 

H t,0  g . Poppoff  et  al . (1978)  point  out  the  fact  that  a 
"long"  wave  (>  255  nrn)  absorption  tail  for  H-,0o  has  beer,  ne- 
glected in  both  CIAP  and  in  Hudson  (1977)*  This  neglect  pro- 
duces an  erroneously  large  estimate  of  H„0o  in  the  lower  strato- 

C.  C. 

sphere.  In  their  model,  CH^0o  is  assumed  to  behave  like  Ho0„ . 
Neglect  of  the  long  wavelength  tail  would  put  both  of  these 
species  at  fairly  large  values  ('v  1 ppbv)  which  Poppoff  et  al. 
suggest  would  be  incompatible  with  experimental  data. 

The  reaction  CIO  + OH  -*■  products . This  reaction,  which, 
if  it  proceeds  according  to 

CIO  + OH  -*■  HC1  + 0., 

at  moderate  rate,  would  be  important  in  stratospheric  chemistry, 
for  two  reactive  species  are  destroyed  (R.E.  Roberts,  IDA, 
private  communication,  1977)  in  the  process. 
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Classically,  a reaction  such  as  the  one  cited  is  assumed  by 
kineticists  to  be  a four-center  reaction  and  too  slow  to  be  of 
consequence.  In  this  case,  however,  the  extreme  energy  release 
involved  (~  55.6  kcal/mol)  and  its  possible  chain-terminating 
character,  suggest  that  it  merits  investigation.  Its  importance 
also  relates  to  CH^  content;  this  reaction  would  be  of  greater 
importance  with  the  recent  CH;j  profile  measurements  of  Ackerman 
et  al.  (1977)  than  with  earlier  profiles  (W.S.  Smith,  private 
communication,  1977). 

The  speculative  reaction  09  + ClO  -*  09Cl0.  Jesson  (1978) 

I Ci  u 

in  a speculative  suggestion,  commented  that  if  the  reaction 
02  + CIO  -+  02C10  takes  place, 

followed  by 

o2cio  h-*v  cio2  + 0 

and 

cio2  h+  CIO  + 0 

a mechanism  for  odd  oxygen  production  would  result,  but  no  evidence 
is  known  for  such  a reaction. 

3.5.3  Other  Modeling  Change1, 

Recent  models  have  been  improved  in  terms  of  their  ability 
incorporate  such  effects  as  multiple  scattering,  diurnal 
averaging,  and  thermal  feedback.  Each  of  these  effects  can  be 
shown  to  be  of  significance;  in  1-D  modeling,  however,  the  appro- 
priateness of  the  various  procedures,  as  in  1-D  modeling  in  gen- 
eral, can  be  debated:  each  represents  an  improvement  on  an  al- 
ready "globally-averaged"  process  of  sorts.  The  appropriate 
means  of  further  averaging  is  never  fully  clear.  In  COMESA 
(1975)  efforts,  1-D  models  were  used  with  fully  diurnally  and 
seasonally  varying  sun,  but  with  fixed  eddy-dif f usivity  profile. 
Currently,  modelers  feel  this  process  uses  too  much  computer 
time,  and  may  be  difficult  to  justify;  instead,  in  what  seems 
1 be  a reasonable  approach,  the  appropriate  time-averaged 
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species  product  concentrations  are  determined  by  simpler  diurnal 
averaging  studies.  The  object,  of  course,  is  to  avoid,  in  ki- 
netic calculations,  the  multiplying  together  of  concentrations 
of  species  which  exist,  only  at  different  tines  of  the  24-hour 
period.  As  far  as  is  known,  however,  all  models  "turn  the  sun 
on"  and  off  over  the  entire  model  altitude  simultaneously  for 
all  wavelengths  in  these  diurnal  averaging  procedures,  thus 
ignoring  differences  in  sunrise  and  sunset  with  altitude  ana 
wavelength . 

Thermal  feedback  effects  also  involve  obvious  questions  as 
to  whether  the  "average"  heating  or  cooling  processes  are  appro- 
priately considered,  or  whether  the  heating  prccc-s.-  • , like  the 
trace  species  in  diurnal  modeling,  art  in  fact  separated  in 
time  and/or  space  from  the  chemical  reaction  processes. 

In  spite  of  these  questions,  in  this  author's  opinion, 
these  "improvements"  to  the  ] -D  models  are  justified  in  indi- 
cating the  potential  importance  of  various  physical  phenomena, 
but  higher  (preferably  3-D)  modeling  will  be  necessary  to  pro- 
vide a more  quantitative  indication  of  their  significance. 

3.5.4  Effects  of  Recent  Chemistry  Changes  on  Model  Atmospheres 

Incorporation  of  the  Howard  and  Evenson  (1977)  rate  con- 
stant in  existing  models  increases  the  ozone  in  the  model  "nat- 
ural" atmosphere.  The  effect  is  evident  in  1-D  models,  but 
rather  easily  dismissed  since  total  ozone  in  a i-D  model  varies 
significantly  with  K profile  choice  as  well  as  with  chemistry, 
and  considerable  ambiguity  exists  as  to  whether  a i-D  model  is 
a global  average,  a mid-latitude,  or  a hemispheric  model.  In  a 
2-D  model,  the  effect  becomes  quite  obvious,  as  shown  in  Fig.  3 . , 
provided  by  Widhopf  (private  communica  ion,  1978).  Note  the  large 
excess  of  ozone  in  the  model  atmosphere  over  that  measured.  This 
model  had  previously  provided  a more  reasonable  match  between 
calculated  and  measured  atmospheres,*  (See  Wi  ...  pf  et  al . , 1977 ) . 
? 

Chlorine  species  :re  no*  included  in  these  mode]  results. 
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FIGURE  3-5.  Widhopf  model  results  showing  how  use  of  the  Howard-Evenson  Reaction  rate 
distorts  the  model  ozone  profile  relative  to  the  measured  atmosphere  values 
Widhopf,  private  communication,  1978.  (See  also  Widhopf  et  a!.,  1977.) 


It  can  be  anticipated  that  the  fast  HO^  + 0^  reaction,  as 
recently  reported  (Zahniser  and  Howard,  1978;  Section  3-5.2), 


will  reduce  the  excess  of  ozone,  but  to  what  extent  is  not  yet 

known.  In  1-D  model  results  (Zahniser  and  Howard,  private 

communication,  1978),  use  of  the  new  HO^  + 0^  rate  which  is 

some  four-fold  faster  than  the  old  rate,  at  230  K,  brought  the 

1 8 2 

1-D  ozone  column  from  10.6  to  9-6  (x  10x  , This  approx- 

imate  10  percent  reduction  apparently  would  not  be  sufficient 
to  bring  the  computed  Wldhopf  model  into  line  with  the  measured 
values,  but  would  be  expected  to  improve  the  situation. 

3.5.5  Recent  Calculations  of  Aircraft  Effects  (1-D) 

a.  The  Different  Models  Compared.  Potential  SST  effects 
on  ozone  have  been  reestimated  in  detail  in  light  of  recent 
(post-1976)  changes  in  chemistry  by  workers  at  three  organiza- 
tions— Lawrence  Livermore  Laboratory  (LLL),  National  Center  for 
Atmospheric  Research  (NCAR)  and  the  National  Aeronautics  and 
Space  Administration  at  Ames  (NASA-Ames).  These  studies  have 
all  used  1-D  models,  but  these  have  varied  in  important  details. 
No  new  studies  have  been  done  in  the  subsonic  question.  Only 
draft  copies  of  the  studies  are  available  to  us  at  the  time  of 
writing.*  A review  as  of  July  1977  was  prepared  by  Broderick 
(1978),  which  is  scheduled  for  early  publication.  The  paper  by 
Broderick  includes  results  as  affected  by  the  Howard  and  Evenson 
(1977)  rate  for  the  HO^  + NO  reaction;  however,  some  effects  not 
included  in  that  paper  have  since  become  available,  and  these 
results  are  emphasized  in  the  following. 

The  characteristics  of  the  three  1-D  models  which  seem 
most  pertinent  to  this  discussion  are  shown  in  Table  3~2.  The 

r 

models  themselves  are  described  in  detail  elsewhere  (Chang  et 
al.,  1977;  Poppoff  et  al.,  1978,  and  Turco  and  Whitten,  1977; 


The  final  version  (Poppoff  et  al. , 1978)  became  available 
during  the  review  period. 
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TABLE  3-2.  MODELS  AND  RUNS  COMPARED1 


Aspect 

LLL 

NASA-Ames 

NCAR/NOAA 

Troposphere  Included 

Yes 

2 

No^ 

No 

Altitude  Range,  km 

0-55 

10-120 

10-65 

Steady  State  (SS)  or 
Time-Dependent  (TD) 

TD 

TD 

SS 

Vertical  Resolution 

1 km 

2 km 

1 km 

Multiple  Scattering  Included 

Yes 

No 

— 

Albedo  Effects  Included 

Yes 

No 

-- 

Diurnal  Averaging  Included 

Yes 

Yes 

-- 

Thermal  Feedback  Included 

Yes 

No 

-- 

Cl  X chemistry  Induced 

Yes 

Yes 

Yes 

Smog  chemistry  Included 

Yes3 

Yes 

Yes 

"Basic"  profile(s)4 

Chang/1976 

WT;  Chang/19765 

Crutzen  et  al , 1978 

Water  Vapor  Changes 

Parameteri zed^ 

Computed 

Computed 

NO^  injection  rates 

Varied 

Varied 

Fixed 

NO  /Ho0  ratio 

X c 

Varied 

Varied 

Fixed 

HgO  injection  (alone) 

Yes 

Yes 

No 

^Refers  to  models  in  the  configuration  used  to  generate  the  main  body  of  re- 
sults discussed  here.  Each  group  has  a number  of  models.  The  various  ef- 
fects listed  as  "included"  can,  of  course,  be  excised  if  desired;  similarly, 
effects  shown  as  not  included  can,  in  general,  be  included  if  desired. 

2 

NASA-Ames  has  a second  model  which  includes  the  troposphere;  they  find  that 
tropospheric  ozone  is  little  affected  by  SSTs. 

^Simplified.  See  Oliver  et  al . , 1977. 

4By  this  is  meant  the  implicitly  preferred  profile.  Results  are  reported 
for  various  profiles  by  each  group. 

5 

"Wofsy-type,"  a profile  similar  to  that  of  these  authors  but  developed  by 
NASA-Ames  based  on  excess  C - 1 4 data.  See  Section  3.3.3.  This  is  referred 
to  elsewhere  herein  as  the  Poppoff  et  al.,  1978  profile. 

5Two  versions  are  used  to  correct  approximately  for  latitude  (see  Section 
3.3.2).  These  are  referred  to  as  Crutzen  and  Howard  (1978)  "high"  and 
"low."  The  mid-stratosphere  portion  is  based  on  Schmeltekopf  et  al.  (1977). 

3Very  recent  calculations  have  included  computed  changes  in  H^O. 
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Crutzen,  1977).  In  a sense,  the  differences  result  from  dif- 
ferences in  approach  and  purpose.  Thus,  the  LLL  group  has  used 
and  modified  its  model  for  a number  of  years  to  explore  a wide 
variety  of  perturbations;  SSTs , atmospheric  nuclear  explosions, 
and  halocarbons  have  all  been  studied  and  sensititives  (parti- 
cularly as  to  kinetics)  explored.  Ho  claim  (to  our  knowledge) 
has  been  made  that  an  "upperbound"  or  other  policy-oriented  re- 
sult is  intended  by  these  calculations. 

The  specific  NASA-Ames  modeling  work  described  here  is  de- 
voted entirely  to  the  SST  problem;  in  this  work,  the  troposphere 
has  little  effect  on  the  result  as  tested  by  them  in  studies 
with  a second  model  which  included  the  troposphere.  They  argue 
that  tropospheric • chemistry  is  uncertain  in  any  event. 

The  NASA-Ames  group  also  argues  against  the  inclusion  of 
thermal  feedback  effects  for  somewhat  similar  (uncertainty) 
reasons.  The  eddy  diffusivity  profile  used  by  them  is  consis- 
tent with  their  approach.  Crutzen  and  Howard  (1978)  also  argue 
that  tropospheric  chemistry  is  poorly  known  and  thereby  justify 
its  omission  for  studying  SST  effects.  A major  purpose  of  the 
Crutzen  and  Howard  paper  is  to  examine  stratospheric  ozone  sen- 
sitivity to  various  reaction  rates  of  current  importance  in 
stratospheric  chemistry.  In  a sense,  both  the  NASA  Arnes  and 
Crutzen  and  Howard  papers  are  consistent  with  the  apparent  in- 
tent of  PL95-95  (see  Section  3-2).  Variations  at  and  near  the 
tropopause,  which  would  be  expected  with  aircraft  cannot,  in 
our  opinion,  be  taken  adequately  into  account  with  a model 
which  fixes  conditions  (concentrations  or  fluxes)  at  the  tropo- 
pause. None  of  the  available  modeling  studies  report  effects 
of  simultaneous  injection  of  SO^ . 

The  various  results  will  now  be  reviewed  briefly. 

b.  NASA-Ames  (Poppoff  et  al.,  1978).  Poppoff  et  al. 
emphasize  effects  on  stratospheric  ozone  due  to  simultaneous 
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injections  of  NO^  and  H„0  at  various  altitudes,  and  for  "near- 
term"  and  "far-term”  NO  emissions  technology.  "Near-term"  en- 
gine technology  implies  today's  combustor  technology;  "far- 
term"  technology  includes  two  possible  NO^  emissions  levels, 
roughly  40  percent  and  15  percent  of  today's  values.  The  ac- 
tual values  used  are  quoted  below.  Their  model,  as  indicated 
above,  did  not  include  a troposphere.  Their  results,  which 
follow,  are  of  particular  interest  in  that  they  show  clearly, 
for  small  to  reasonable  SST  fleet  sizes,  the  change  in  signi- 
ficance of  H00  and  N0x  emissions  due  largely  to  the  revision  of 
the  KO^  + NO  reaction  rates. 

Two  aircraft  types  were  considered  by  Poppoff  et  al . , 1978, 
with  varying  N0x  emission  indices,  as  shown  in  Table  3~3. 


TABLE  3-3.  AIRCRAFT  AND  FLEET  CHARACTERISTICS 
( Poppoff  et  al . , 1 978) 


Nominal* 


Type 

Nominal 
Mach  No. 

Nomi nal 
Altitude, 
km 

Nominal  Fuel 
Flow  Rate, 
kg/ hr/a i rcraft 

N0X  Emission 
Index , 

9 N02/kg 

Fleet 

108 

N0X 

Emissions 

kg/yr 

H20 

A 

2.7 

20.0 

37,800 

15.6 

1.47 

125 

6.0 

0.57 

125 

2.0 

0.19 

125 

B 

2.3 

17.5 

35,200 

18.0 

1.65 

117 

7.0 

0.64 

117 

3.0 

0.27 

117 

*100 

Ai rcraft, 

7 hours/day. 

From  Poppoff  et 

al.  (1978) 

The  reported  ozone  depletions  are  given  in  Table  3“ 4.  Re- 
sults are  given  for  the  "Dickinson-Chang"  profile  (or  Chang/ 
1976),  and  for  "WT"  (Wofsy-type)  profile. 
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TABLE  3-4.  OZONE  CHANGES  CALCULATED  FOR  100  AIRCRAFT 
WITH  DIFFERENT  NO  EMISSION  INDICES 

% Ozone  Change* 


Type 

Cruise 

Altitude, 

km 

N0X  E.I., 
g NO2/ kg 

Injections , 

103  kg/yr 

N0X  only 

N0X  plus  H2O 

N0X 

H20 

DC1 

WT2 

DC1  WT2 

A 

20.0 

15.6 

1.47 

125 

0.12 

0.12 

0.07  0.03 

6.0 

0.57 

125 

0.05 

0.05 

0 -0.06 

2.0 

0.19 

125 

0.02 

0.02 

-0.03  -0.08 

B 

17.5 

18.0 

1.65 

117 

0.13 

0.13 

0.10  0.06 

7.0 

0.64 

117 

0.05 

0.05 

0.02  -0.02 

3.0 

0.27 

117 

0.02 

0.02 

-0.01  -0.05 

The  Ames  authors  point  out  that  these  near-zero  numbers,  obtained  by  inter- 
polation, involve  substantial  (fractional)  uncertainty. 

^This  is  referred  to  as  the  Dickinson-Chang  profile;  it  has  been  referred  to 
elsewhere  as  the  "Chang/1976"  profile  (or  at  times  as  the  "new"  Chang  pro- 
file). 

2 

This  is  the  NASA-Ames  "Wofsy-type"  profile,  developed  by  NASA-Ames  (see 
Section  3.3.3). 


The  results  given  in  Table  3“ 4 are  plotted  in  Fig.  3~  6 for 
the  two  K„  profiles.  Note  that  the  ozone  column  increases  with 
increased  NO  injection  at  constant  water  injection  rate;  the 
depletion  appears  to  be  due  to  the  H^O  injections.  The  NO^  re- 
sult is  entirely  opposite  to  that  found  in  CIAP.  The  HnO  effect 

CL 

on  ozone  was  never  clear  in  CIAP;  in  one  paper  however  (Crutzen, 
1974),  it  was  reported  that  water  vapor,  according  to  then- 
current  chemistry,  actually  increased  ozone  slightly.  Prior  to 
CIAP  and  the  NO  issue,  however,  as  noted  in  an  earlier  section, 
water  vapor  depletion  of  ozone  was  a major  question  with  regard 
to  the  SST  (Section  3.3-7). 

Poppoff  et  al.  (1978)  also  give  effects  as  a function  of 
N0x  and  H~.0  injection  rate  at  20  km  in  Fig.  3 — 7 , and  in  Fig.  3“8 
as  a function  of  altitude,  for  a specified  NO.  injection  rate. 

X 

In  Fig.  3“7,  the  water  effects  and  NO  effects  are  computed  sep- 
arately. The  scales  are  sized  however  for  an  NO  emission  index 

X 
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FIGURE  3-6.  NASA-Ames  1-D  model  results  plotted  here  to  show  effects  of 
varying  N0X  injection  rates  at  constant  H2O  injection  rates. 
The  water  rates  correspond  to  100  SSTs  at  17.5  and  20  km  for 
Dickinson-Chang  (DC)  and  Wofsy-type  (WT)  profiles.  The  points 
plotted  involve  approximations  so  the  curves  are  not  as  pre- 
cisely located  as  might  be  inferred. 

Data  Source:  Poppoff  et  al.,  1978  (see  Table  3-4) 


0 


30 


5 10  15  20  25 

GLOBAL  NOx  INJECTION  RATE,  108  kg  N02/yr 

FIGURE  3-7.  Separate  effects  of  water  and  N0X  injections  at 
20  km  at  various  rates  on  the  ozone  column  above 
10  km  as  calculated  at  NASA-Ames  (Poppoff  et  al., 
1978)  for  three  Kz  profiles.  The  water  abscissa 
and  the  N0X  abscissa  are  coincident  for  a N0X 
emission  index  of  6 gm/kg  (as  NO?)  and  water 
emission  index  of  1.3  kg/kg.  For  other  N0X 
emission  indices,  the  water  effect  should  be 
scaled;  e.g.,  the  NASA  "current"  technology  SST, 
the  water  effect  is  6/ 15.6  or  0.38  times  the  value 
read  at  a given  N0X  abscissa. 

Data  Source:  Poppoff  et  al . , 1978 
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' IGURE  3-8. 


Effects  of  N0X  and  N0X  and  H2O  injections  at 
various  altitudes  on  the  ozone  column  above 
10  km  as  calculated  by  NASA-Ames  (Poppoff  et 
al.,  1978).  Results  shown  are  for  a N0X  in- 
jection of  7 x 10°  kg  (as  N0?)/year  and  for 
this  N0X  rate  plus  1.5  x 10"  kg  b^O/year, 
for  two  Kz  profiles.  A N0X  emission  index  of 
6 gm/kg  fuel  (as  NO2)  is  implied. 

Data  Source:  Poppoff  et  al . , 1978 
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Note  also  that  according  to  these  results,  an  optimum  Nt;x 
emission  index  should  exist,  varying  with  selected  K profile 
and  fleet  size,  to  provide  minimum  effect  on  the  ozone  column. 
This  optimum  is,  of  course,  a "1-D"  optimum,  and  should  not  be 
given  too  much  credence.  It  is  not  obvious  from  these  results 
that  very  low  emission  index  values  are  a necessarily  attrac- 
tive goal. 

Poppoff  et  al.  (1978)  give  an  extensive  discussion  of  un- 
certainties and  effects  of  various  rate  constants.  They  point 
out  in  particular  the  sensitivity  of  the  results  to  the  rate 
constant  used  for  the  HO^  + 0^  reaction.  Their  results  are  in- 
cluded as  Table  3~5 • Note  the  always  deleterious  effects  of 
water  vapor.  For  details  other  than  those  shewn  in  the  foot- 
notes, the  reader  is  referred  to  Poppoff  et  al.,  1978.  They 
also  point  out  quite  convincingly  that  estimates  of  uncertain- 
ties are  unreliable  and  apt  to  be  misleading. 

c.  Crutzen  and  Howard  (1978).  Crutzen  and  Howard  studied 
ozone  changes  due  to  injections  of  NO  , N2O,  and  chlorof luoro- 
methanes.  As  noted  in  Table  3~2,  a steady-state  stratospheric 
model  was  used.  A combination  of  concentration  (N„0,  CH ,, , CO, 

d.  n 

Ho0,  H9,  CH3CI)  and  flux  (odd  nitrogen,  Cl.,  ozone)  boundary 

c.  u-  J X 

conditions  was  used  at  10  km,  the  lower  boundary.  A major  pur- 
pose was  to  investigate  effects  resulting  from  changes  in  the 
rate  used  for  the  H0o  + NO  reaction,  and  to  consider  effects 
due  to  other  reaction  rate  uncertainties. 

The  authors  chose  equinox  at  40°  latitude  as  being  repre- 
sentative of  global  conditions.  Two  eddy  aiffusivity  profiles, 
differing  only  below  15  km  were  used  to  represent  global,  averay 
(K„  high)  or  conditions  poleward  of  40°  (K^  low)  as  noted  earli 
(See  Section  3*3*4).  An  injection  rate  (over  1 km)  of  1000  NO 
molecures  cm  sec  and  2 x 10^  H2O  molecules  cm  Jsec  was  use 

in  all  cases;  this  corresponds  to  an  injection  rate  of  1.23 

9 10 

x 10  kg  NO2  per  year  and  9*64  x 10  kg  Ho0  per  year;  a fuel 
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TABLE  3-5.  SENSITIVITY  STUDIES  OF  OZONE 
PERTURBATIONS  DUE  TO  SSTsl 
(Source:  Poppoff  et  al.,  1978) 


Model  Modification 


Stratospheric  ozone  column 
changes  due  to  SST  N0X  injection 
without/with  water  vapor  effect‘d 


Increase  background  N~0  by  50% 

^ 3 

Increase  total  chlorine  to  5 ppbv 

Eliminate  chlorine  nitrate  formation 

Eliminate  scattered  radiation 

Eliminate  diurnal  effects 

Use  Moortgat  CH^O  quantum  yields^ 

Include  NO^NO^  formation 

Decrease  the  H0~  + 0,  rate  by  a factor  of  3^ 
Increase  the  HO^  + 0^  rate  by  a factor  of  3 
Increase  the  OH  + HN03  rate  by  a factor  of  36 
Decrease  ambient  H^O  from  1.0  to  2.5  ppbv  at  14  km 


-0.35%/  -0.88% 
+2.13 %/  +0.42% 
-0.80%/  -0.81% 
+0.55%/  -0.54% 
+0.30%/  -0.50% 
+0.68%/  -0.38% 
+0.17%/  -0.74% 
-0.50%/  -1.10% 
+1.80%/  +0. 13%  * 
-1.74%/  -1.90% 
-0.07%/  -1.10% 


A global  N0x  injection  rate  of  7 x 10  kg  NO^/yr  and  1.5  x 10  kg  H^O/yr 
at  20  km  are  assumed;  the  WT  diffusion  coefficient  has  been  used  for  each 
sensitivity  test.  See  Poppoff  et  al . (1978)  for  detailed  discussion  of 
these  effects. 

2 

For  reference,  the  nominal  ozone  column  changes  are  +0.40%/  -0.62%, 
respectively. 

3 

An  increase  of  about  3 ppbv. 

^Represents  an  increase  of  about  50%  (varying  with  altitude)  in  the  yield  of 
hydrogen  radicals  over  baseline  model  values. 

JThe  baseline  HO2  + O3  rate  is  7.3  x 10  ^ e ^^cm'Vsec . Note  that  the 
corresponding  tripled  rate  at  230  K is  8.6  x 10'^  cm^/sec,  a rate  which  is 
lower  than  the  value  of  1.1  x 10“^  cm^/sec  calculated  from  the  Zahniser 
and  Howard  1978  expression  (Section  3.5.2a). 

^The  baseline  OH  + HN0.  rate  is  8.0  x 10"^  cm^/sec. 
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consumption  of  7.71  x 10  ^ kg/year  and  an  NO  emission  index  of 
16.0  g NO^/kg  fuel  is  also  implied.  The  fuel  flow  corresponds 
to  684  of  the  type  A aircraft  (at  7 hours  per  day)  described  in 
the  previous  section.  Injections  at  13,  16,  and  20  km  were 
studied . 

The  effects  due  to  variations  in  three  reaction  rates  and 
the  two  eddy  diffusivity  profiles  were  computed.  The  three  reac- 
tions and  rate  expressions  used  are  as  follows,  using  nomencla- 
ture as  in  Section  3.5-2. 


H02  + OH 

" °2  + 

h2o 

(R19) 

ho2  + o3 

-*  OH  + 

C\J 

O 

OJ 

(R14) 

H0o  + NO 

-*•  OH  + 

OJ 

o 

(R26) 

ho2  + ho2 

H2°2 

OJ 

o 

+ 

(Rl8 ) 

Reaction  rate  effect  studies  used  the  following  values  in  cm  / 
molecule-sec.  The  value  at  230  K is  shown  as  a representative 
stratospheric  value. 


Expressions 


k 230  K 


k19(ST) 
k19 (hi ) 
kl4(ST) 
kl4 (hi ) 
k26 (HE ) 
k26(C) 
kl8 (ST) 
kl8 (hi ) 


2.1 

x 10 

11 

same 

5.1 

x 10“ 

11 

same 

10“ 

13  exp 

(-1525/K) 

1.3 

X 

lO-10 

2 x 

10“13 

exp 

(-1250/K) 

8.7 

X 

10“16 

3 x 

10-11 

exp 

(-390/K) 

5-5 

X 

10"12 

5 x 

io“12 

exp 

(-450/K) 

7.1 

X 

10“13 

5 x 

10“12 

exp 

(-500/K) 

5.7 

X 

10“13 

1.5 

x 10“ 

11  exp  (-500/K) 

1.7 

X 

io“12 
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Here  "ST"  refers  to  "standard"  and  "hi"  to  a "high"  value; 

HE  refers  to  the  Howard  and  Evenson  measurement  (using  a tempera- 
ture coefficient  consistent  with  an  A factor  of  3 x lO-^  and 
their  value  of  8 x lO-"^  cm^/sec  at  296  K);  and  "C"  refers  to  an 
expression  based  on  the  Cox  (1975)  measured  value. 


The  Crutzen  and 

Howard  results  are 

summarized 

in  Table 

3-6 

TABLE 

3-6 

. STRATOSPHERIC 

OZONE  COLUMN  CHANGES 

, PERCENT 

Case 

1 

2 

3 

4 

5 

6 

k19 

ST 

ST 

hi 

ST 

hi 

ST 

k14 

ST 

ST 

ST 

hi 

hi 

ST 

k26 

C 

HE 

HE 

HE 

HE 

HE 

He 

ST 

ST 

ST 

ST 

hi 

hi 

A1 ti tude 
Injection, 

km 

Kz,  HIGH- 

-POLEWARD  40 

0 

13 

-0.23 

+0.07 

-0.20 

+0.71 

+0.22  +0.08 

16 

-1.15 

-0.37 

-1.03 

+ 1.21 

0 

■0.39 

20 

-3.08 

-1.55 

-2.75 

+ 1.11 

-1.01 

■1.64 

Kz,  L0W--GL0BAL 

13 

-0.81 

-0.26 

-0.91 

+ 1.67 

+0.40 

■0.24 

16 

-2.06 

-1.05 

-2.17 

+2.19 

+0.06 

1.05 

20 

-3.95 

-2.37 

-3.91 

+ 1.82 

-1.04 

■2.41 

The  effect  of  the  Howard  and  Evenson  vs.  Cox  measurement 

of  Reaction  26  is  seen  by  comparing  cases  1 and  2;  the  faster 

rate  results  in  smaller,  or  opposite  sign,  effects  on  the  ozone 

column.  Note  that  the  actual  rate  used  in  these  calculations 

is  less  than  the  recommended  NASA  rate  (Hudson,  1977)  of  8 x 
-12 

10  , which  assumed  no  temperature  correction. 
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The  effects  of  the  rate  of  Reaction  19  can  be  seen  by  com- 
paring Runs  2 and  3.  Note  that  the  standard  value  for  this  reac- 
tion (Case  2)  is  taken  by  the  authors  to  be  2.1  x 10-11;  the  re- 
commended value  by  NASA  (Hudson,  1977)  is  3 x 10-11  and  the  only 
direct  measurement  (Burrows  et  al.,  1977)  is  5.1  x 10-1"  (which 
is  used  in  Case  3).  Note  that  the  rate  for  this  reaction  was 
taken  as  2 x 10  ^ in  CIAP,  resulting  in  large  ozone  depletions. 
Case  3 would  seem  to  be  reasonable  as  a perhaps  conservative 
estimate . 

The  powerful  effect  of  Reaction  14  can  be  seen  by  comparing 
cases  2 and  4,  with  a fast  rate*  resulting  in  ozone  enhancements 
at  all  altitudes.  The  reasons  for  this  were  discussed  in  Sec- 
tion 3 - 5 - 2 . Note  that  the  rate  recommended  in  Hudson  (1977  ) 

_ -i  r 

at  230  K is  2.9  x 10~  , close  to  the  geometric  mean  (3-3)  of 

the  two  values  (230  K)  developed  from  the  Crutzen-Howard  ex- 
pressions. Duewer  et  al.  (previous  section)  use  an  expression 
which  gives  4.4  x 10-^  cm^/sec  at  230  K. 

The  small  effects  of  the  rate  of  Reaction  18  can  be  seen 
by  comparing  cases  2 and  6. 

The  effects  of  using  all  high  values  vs.  all  low  values 
can  be  seen  by  comparing  cases  1 and  5-  Assuming  k^g  has  little 
effect.  Case  5 can  be  compared  to  Case  4 to  again  show  the 
strong  effect  of  Reaction  19. 

The  effect  of  a reduced  effective  transport  rate  at  the 
tropopause  can  be  seen  by  comparing  the  K high  to  the  Kr?  low 
values.  As  most  aircraft  traffic  is  poleward  of  40u,  the 
(high)  values  may  be  most  appropriate,  although  the  interpreta- 
tion of  these  results  is  subject  to  some  debate. 


The  "fast  rate"  used  here  is  still  slower  than  the  1.1  x 10 
cm3/sec  value  found  by  Howard  (1978).  See  Section  3-5.2. 
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Crutzen  and  Howard  also  show  the  vertical  redistribution 
of  ozone  that  takes  place  in  Case  4 for  1 6 km  injection  (1.21 
percent  net  increase).  In  this  case,  ozone  was  enhanced  from 
10  to  about  27  km,  above  which  ozone  was  depleted.  Maximum 
percentage  local  ozone  enhancement  was  8.36  percent  at  10  km, 
the  lower  model  boundary. 

Crutzen  and  Howard  conclude  that  ozone  enhancement  from 
SST  injections,  accepting  the  Howard-Evenson  measurement,  would 
result  only  with  some  combination  of  a fast  rate  for  HC^  + 0^ 
and  a slow  rate  for  HC^  + OH.  The  high  rate  for  HO^  + 05  has, 
in  fact,  since  been  reported  (Zahniser  and  Howard,  1978),  as 
discussed  earlier. 

d . Lawrence  Livermore  Laboratory  (LLL)  Results  (Luther, 
1977).  These  workers  have  continued  their  extensive  studies 
of  atmospheric  perturbations  due  to  aircraft,  halocarbons,  and 
nuclear  weapons  tests,  etc.  Several  aspects  of  their  work  are 
of  particular  interest  here. 

Historical  changes  in  aircraft  ozone  depletion  results 
with  changes  in  chemistry  and  K profile.  The  LLL  workers 
break  their  reporting  of  this  work  into  two  periods:  197^-1976, 
1976-September  1977- 

The  LLL  results  for  the  197^-1976  period,  in  which  chlorine 
chemistry  was  not  included,  are  shown  in  Fig.  3-9  For  17-  and 
20-km  injections.  The  Chang/197^  profile  results,  with  197 ^ 
chemistry,  represent  the  starting  point;  the  second  point  (mov- 
ing down)  shows  changes  resulting  from  the  use  of  the  Chang/ 

1976  Kz  profile  (referred  to  by  NASA-Ames  in  the  previous  sec- 
tion in  the  Dickinson-Chang  profile),  with  the  same  chemistry. 
Other  results,  using  the  Chang/1976  profile,  with  changes  in 
chemistry,  are  shown.  The  biggest  single  effect  in  this  time 
frame  was  clearly  the  result  of  the  revision  of  the  HO  + HO^ 
rate  (R19)  from  2 x lO-"^  crrr/sec  in  CIAP  to  2 x 10  cm^/sec, 
their  then-preferred  value.  Note  again  that  a somewhat  higher 


and 


value  (3  x 10  cm'vsec)  is  recommended  in  Hudson  (1977  ), 

that  both  groups  reject  the  only  direct  measurement  cl  5 i >• 

-11 

10  (see  Section  3.4.2).  Details  of  the  other  rate  changes 
are  given  in  Luther  (1977). 


Change  in  total  ozone  - % 


FIGURE  3-9.  Changes  in  calculated  ozone  depletions 
over  the  1974-1976  period  due  to  various 
changes  in  modeling  parameters.  Global 
injection  rate  of  N0X  (as  NO?)  taken  as 
2.46  x 109  kg/yr.  The  Kz  change  refers 
to  the  Chang/1976  profile.  See  Luther, 

1977,  for  details.  See  also  Broderick, 

1978. 
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Following  this  work,  a number  of  changes  were  made  to  the 
LLL  model,  including  the  addition  of  multiple  scattering,  in- 
clusion of  chlorine  chemistry,  and  modification  of  some  addi- 
tional reaction  rates  and  absorption  cross  sections  in  accord- 
ance with  the  recommendations  in  Hudson  (1977).  (A  rate  of 
-11  3 

2 x 10  cm  /sec  was  retained  in  these  modeling  exercises  as 

1 i 

the  then-preferred  LLL 'value;  later  runs  have  used  3 x 10  x^.) 
These  changes  resulted  in  an  increase  in  ozone  reduction  for 
the  same  injection  rate  (2000  mol/cm  -sec  over  1 km;  2.46  x 

Q 

10  kg/yr  as  NO^)  from  about  1.0  percent  at  17  km  and  4.0  per- 
cent at  20  km  as  in  Fig.  3~9  to  1.3  percent  at  17  km  and  4.8 
percent  at  20  km,  respectively.  At  this  point,  introduction 
of  the  Howard  and  Evenson  rate  constant  had  dramatic  effects, 
as  shown  in  Table  3“7,  which  table  also  includes  LLL  results 
using  the  h’unten/1975  profile.  Note  that  these  results  in- 
clude an  assumed  temperature  coefficient  for  the  Howard  and 
Evenson  measurement,  reducing  the  rate  at  a nominal  strato- 

__  i p 

spheric  temperature  (230  K ) to  4 . 9 x 10  J"  rather  than  8.1  x 
-12 

10  , as  recommended  in  Hudson  (1977),  and  presumably  reducing 

the  magnitude  of  the  change  in  effect;  however,  the  slower  rate 
used  for  the  OH  + HC^  reaction  works  the  other  way. 


TABLE  3-7.  OZONE  COLUMN  CHANGES  SHOWING  EFFECTS  OF  THE  HOWARD 
AND  EVENSON  (1977)  REACTION  RATE  MEASUREMENT.  GLOBAL 
N0X  INJECTION  2.46  x 10^  kg/yr  (Source:  Luther,  1977) 


Profile 

N0X  Injection 
Altitude  (km) 

Change  in 

0j  Column, 

Change 

in  NO  , 

X 

Old  Rate* 

HE  Rate* 

Old  Rate 

New  Rate 

Chang  (1976) 

20 

- 4.79 

0.55 

46.7 

47.9 

17 

- 1.31 

1.96 

27.3 

27.8 

Hunten  (1975) 

20 

-10.8 

-6.90 

96.2 

97.1 

17 

-4.35 

0.83 

61.4 

62.3 

★ 

The  "old  rate" 

used  is  2 x 10" 

1 3 3 

cnT/sec;  the 

HE  rate  is  based 

on  the 

Howard  and  Evenson  1977  measurement,  with  an  assumed  temperature 
coefficient  [4.28  x 10"^  exp  (-500/T)]. 
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Injection  rate  effects  and  modeling  comments.  In  studies 
under  CIAP,  ozone  depletions  calculated  using  1-D  models  were 
generally  found  to  be  roughly  linear  with  injection  rate.  Thus, 
the  Hunten  formula  involved  the  expression 

6 = l.H05x  " 0.0105X2 

where 

6 = percent  ozone  depletion 
X = incremental  mixing  ratio  of  odd 

nitrogen,  ppbv,  which  is  proportional 
to  injection  rate 

In  the  range  of  interest,  the  linear  term  dominated.  Sim- 
ilarity, in  correlating  Chang's  1 97 3 — 7 *4  results,  expressions 
were  used  of  the  type 


where  x^  = annual  injection  rate  and  a and  b are  constants. 
Values  of  b ranged  from  1.0  at  9 km  to  0.89  at  23  Km.  (See 
pp.  3-32,  Oliver  et  al.,  1977). 

With  1976  studies  (see  p.  3-70,  Oliver  et  al.,  1977)  using 
the  Chang/1976  profile,  ozone  column  changes  were  found  to  be 
linear  with  injection  rate  at  17  km;  however,  nonlinearities 
were  evident  with  the  Hunten  profile.  Introduction  of  the 
Howard-Evenson  rate  for  the  HO^-NO  reaction  has  now  made  non- 
linearity even  more  pronounced.  This  is  shown  by  the  LLL  data 
given  in  Fig.  3~10  for  17-  and  20-km  injections  of  NO  (and  by 
the  NASA-Ames  results  given  earlier);  the  changes  since  the 
previous  report  (Oliver  et  al . , 1977)  are  shown  for  convenience. 
Water  injections  are  not  included  in  Fig.  3 “10  nor  are  thermal 
feedback  effects.  Note  that  the  total  enhancement  with  either 
profile  is  less  than  1 percent  for  20-km  injections  and  (appar- 
ently) about  2 percent  for  17-km  injections.  Note  also  the 

scale  of  the  Injections:  1000  molecules/cm^-sec  over  1 km  al- 

9 

titude  corresponds  to  1.23  x 10  kg  NO^  as  NO-,  per  year.  With 
the  yearly  NO  emission  figures  (NAS,  1975)  cited  for  Concorde, 
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OZONE  COtllMN  CHANGE,  percent 


HANG  K 


• CHANG,  K2;  1976 


• HUNTEN,  Kz;  1977 


• HUNTEN,  K2;  1976 

• CHANG,  Kz;  1977 


• HUNTEN,  K2;  197' 


’••  CHANG,  K2;  1976 


HUNTEN,  K2;  1976 
n ppb  CLX' 


II#  HUNTEN,  K2,  1976 
2000  (NO  CLX) 


NO  INSERTION  RATE,  MOL  cm  -SEC  OVER  1 km  BAND 


FIGURE  3-10 


Ozone  column  changes  for  17  and  20  km  N0X  injections  calculated 
at  Lawrence  Livermore  Laboratory,  using  their  1-D  model,  for 
two  eddy  diffusivity  profiles  (Chang,  1976  and  Hunter,  1974) 
with  1976  and  1977  chemistry  sets.  The  1976  values  are  taken 
from  Table  3-20  (p.  3-65)(01iver  et  al . , 1977)  for  an  atmos- 
phere without  chlorine.  A number  of  changes  were  made  in 
going  to  the  1977  results. 


[ 

20  km 

I 

this  corresponds  to  about  2000  aircraft  on  a global  basis,  or 
about  1000  in  the  Northern  Hemisphere,  if  the  ozone  depletion 
is  assumed  to  be  that  in  the  "corridor."  A figure  of  500  mole- 
cuies/cm  -sec  at  20-km  corresponds  to  375  (NAS,  1975)  "large 
SSTs"  on  a global  basis. 

An  interesting  point  now  emerges  which  can  be  studied  only 


in  a 2-D  model. 


?his  is  that  injection  conditions  could  exist 


which  on  a global  average  basis  would  imply  enhancement;  how- 
ever, this  might  consist  of  a "corridor"  in  which  depletion 
occurs,  and  a larger  area  elsewhere,  where  enhancement  occurs. 
Obviously,  higher  dimension  model  studies  are  needed. 

Additional  important  points  result  from  these  LLL  calcula- 
tions. First,  as  has  been  mentioned  earlier,  the  NO  injection 

alters  the  altitude  distribution  of  NO  . This  is  shown  in  Fig. 

x 

3-11,  which  shows  percentage  changes  in  local  ozone  concentration 
for  three  injection  rates  with  the  Chang/1976  K„  profile.  The 
ozone  concentration  is  increased  at  low  altitudes  and  decreased 
at  higher  altitudes.  The  absolute  concentration,  of  course,  varies 
with  altitude  so  that  percentage  changes  in  concentration  cannot 
be  interpreted  directly  in  terms  of  column  changes.  Note  also 
that  even  though  the  net  effect  on  the  ozone  column  might  be  zero 
in  a 1-D  model,  a climatic  effect  could  result  from  the  altitude 
redistribution  of  ozone  as  well  as  from  the  added  NCU , Ho0  and 

(-  C. 

S0o.  A final  important  point,  already  made,  is  illustrated  in 
Fig.  3~11,  in  which  ozone  concentration  changes  near  the  tropo- 
pause  (which  in  this  author's  opinion  are  best  developed  in  a 
model  such  as  this  one,  i.e.,  one  with  a troposphere"',  are  given 
for  the  LLL  1-D  model  with  several  NO^  injection  rates  at  20  km. 

Multiple  scattering  and  surface  albedo  effects.  The  LLL 
authors  studied  the  effects  of  including  or  not  including  mul- 
tiple scattering  and  of  changes  in  surface  albedo.  Results  for 
NOx  injections  were  generally  insensitive  (approximately  1 per- 
cent increase)  to  these  changes,  although  in  one  near-zero  case, 
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the  sign  of  the  effect  reversed  (from  -0.07  percent  with  pure 
absorption  to  +0.39  percent  with  multiple  scattering  and  a sur- 
face albedo  of  0.25.  For  halocarbon  injections,  however,  the 
effect  of  multiple  scattering  is  quite  significant,  increasing 
the  depletion  estimate  by  about  20  percent. 


FIGURE  3-11.  Percentage  changes  in  ozone  vs. 

altitude  for  three  injection 
rates  of  N0X  at  20  km.  The  rates 
shown  are  in  molecules/cm3  sec 
over  a 1 km  band,  and  correspond 
to  global  injection  rates  of  0.62, 

1.23,  and  1 .85  x 109  kg  N0X  (as 
NO2)  per  year.  Chang/1976  Kz 
profile. 

Source:  Lawrence  Livermore  Laboratory,  1978 


Water  vapor  injections  and  thermal  feedback  effects.  As 
was  found  earlier  by  Liu  et  al . , 1976,  increases  in  the  strato- 
spheric water  vapor  were  noted  to  lead  to  a decrease  in  ozone 
in  the  tlASA-Ames  modeling  work  (Poppoff  et  al . , 1978),  as  al- 
ready discussed.  These  effects  have  been  further  investigated 
by  ILL,  including  studies  of  thermal  feedback  effects,  which 
later  were  not  included  in  Poppoff  et  al.  (1978). 
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The  LLL  results,  including  feedback  effects  for  parame- 
terized water  contents  with  the  Chang/1976  profile,  and  re- 
vised chemistry  (from  that  used  in  Figs.  3-10  and  3-11),  are 
shown  in  Fig.  3**12.  Details  of  the  chemistry  used  are  given  by 

Luther  and  Duewer  (1973);  important  changes  from  the  1977  chem- 

_ ^ 2 

istry  used  in  Figs.  3-9  and  3-10  include  use  of  8.0  x 10~x 
3 -1 

cm  -sec  , independent  of  temperature  for  the  HG„  + NO  reaction, 

-5  a 


,-H 


cmJ-sec  ^ for  the  H0„  + OH  reaction.  An  NO  in- 

2 x 


and  3 x 10 

jection  of  1000  mol/cmJ-sec  at  20  km  was  used  for  the  injection 
studies.  The  results  show  that  increases  in  water,  with  or 
without  NO  injections,  reduce  the  ozone  column;  however,  with 
thermal  feedback,  the  effect  on  ozone  of  the  increase  in  water 
content  is  largely  eliminated. 


The  LLL  invest igators  also  report  earlier  computations  of 
water  effects  using  other  chemistry  sets,  for  both  the  Chang/ 

1976  and  Hunten/1975  K„  profiles.  Resulting  effects  were  simi- 
lar for  the  Chang/1976  profile  but  differed  somewhat;  the  ozone 
was  depleted  by  the  degree  cf  depletion  decreased  with  increased 
H.;0  with  thermal  feedback  for  the  Hunten/1975  case.  All  such 
results  are,  of  course,  sensitive  to  the  temperature  coefficients 
used  for  the  reaction  rates  and  many  of  these  are  poorly  estab- 
lished . 


Wuebbles  (1978)  reports  on  later  calculations  carried  out 
at  LLL , in  which  the  water  increase  was  calculated  rather  than 
parameterized.  The  reaction  rates  used  were  those  suggested  by 
Hudson  (1977),  which  are  much  the  same  as  those  used  in  the 
computations  shown  in  Fig.  3-12.  A CIO  content  of  1.7  ppbv 
("current"  levels)  was  assumed,  along  with  the  Chang/1976 • pro- 
file. The  computations  included  multiple  scattering,  but  did 
not  include  thermal  feedback.  Results  are  shown  in  Table  3~8. 
Note  that  water  and  NO,  injection  effects  on  ozone  at  17  km  are 

X 

roughly  additive,  but  somewhat  less  than  additive  at  20  km. 

No  results  with  the  new  model  are  available  for  the  Hunten 
K„  profile. 
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Change  in  — ppmv 


FIGURE  3-12.  Lawrence  Livermore  Laboratory  modeling  results 
showing  water  effects  on  the  ozone  column  with 
and  without  consideration  of  thermal  feedback 
effects.  Chang/1976  Kz  profile.  N0X  injection 
rate  1000  mol/cm^-sec  over  1 km.  For  detailed 
chemistry,  see  Luther  et  al , 1977;  critical  rates 
include  (in  cm3/molecule-sec)  8.0  x 10-12  for 
H0o  + NO,  e x 10'11  for  H02  + OH,  and  7.3  x 10“14 
exp  (-1225K)  for  HO?  + 03  . 


TABLE  3-8.  EFFECTS  OF  COMBINED  NOx  - H?0  INJECTIONS 
ON  THE  OZONE  COLUMN  CHANG/1976  Kz  PROFILE; 
HUDSON,  1977  CHEMISTRY. 

(Source:  Wuebbles,  1978) 


Ozone  Column 

Change,  % 

Injection  Rate 

17- km 

20-km 

mo1/cm9-sec* 

Injection 

Injection 

NO..  1,000 

A 

H90  0 

1.57 

1.47 

NO 

h2o 

0 

177,000 

-0.16 

-0.47 

NO 

h2o 

1,000 

177,000 

+1.42 

+1.11 

NO 

X 

1,000 

+ 1.13 

+0.42 

h2o 

531,000 

* 3 

Over  a 1-km  altitude  band,  1000  mol/cm  -sec  of  N0X  corresponds 
to  1.23  x 109  kg/yr  (as  NO?);  177,000  mol/cm^-sec  of  H?0 
corresponds  to  8.5  x 10 10  kg/yr  of  HjO. 
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3.5.6  2-D  Model  Results 


IIo  2-D  nodel  results  incorporating  botn  the  Howard-Evenson 
rate  and  chlorine  chemistry  are  available.  A 2-D  evaluation 
including  these  is  underway  at  NASA-Ames;  the  evaluation  will 
include  Widhopf  modeling  results  as  part  of  their  study.  Widhopf 
and  Glatt.  (197')  dc  report  briefly  on  runs  using  the  Howard- 
Evenson  rate  data,  along  with  ether  revisions,  but  not  includ- 
ing the  very  recent  HCr  + C-,  measurement  or  chlorine  chemistry. 
The  run  used  a model  fleet  of  both  subsonic  and  supersonic  air- 
craft and  included  both  H„0  and  NO  injections;  the  fleet  was 

(-  X 

that  reported  in  Oliver  et  al.,  1977,  except  that  the  NO  injec- 

X 

Lien  rate  above  15  km  (the  SST  cruise  component)  was  tripled. 

The  crime  differences  from  earlier  results  (Widhopf  et  al.,  1977) 
was  an  increase  in  the  total  ozone  column  over  and  above  that 
computed  earlier.  The  total  ozone  column  increased  at  40°N  by 
about  3-5  percent  of  the  background  in  the  fall  vs.  about  1.2 
percent  in  earlier  runs.  A minimum  in  the  increase  was  found 
in  the  winter  with  about  1 percent  in  the  current  result  vs. 
about  0.4  percent  earlier.  The  model  suffered,  as  noted  earlier, 
from  an  excess  f ozone  in  the  "natural"  atmosphere. 

3.6  VALIDATION  ATTEMPTS 

Studies  of  mrdel  estimates  of  ozone  depletion  for  "known" 
events  for  use  in  comparison  with  the  ozone  record  have  contin- 
ued. These  include  effects  of  t he  solar  proton  event  of  August 
1972,  of  atmospheric  nuclear  weapons  tests  in  the  early  1960s, 

and  estimates  •_  f ozone  i * about  40  km  due  to  increas- 

ing halocarbon  production. 

The  August  1972  solar  proton  event  and  its  effect  on  ozone 
near  4 mbar  was  described  and  compared  to  satellite  observations 
by  Crutzen  et  al.  (1975)  and  by  Heath,  Krueger,  and  Crutzen 
(1976).  The  work  has  been  updated  with  current  chemistry  by 
Borucki  et  al.  (1978)  In  a 2-D  model.  In  the  Crutzen  et  al . 
studies,  reasonable  agreement  (with  some  ambiguities)  between 
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theory  and  observation  could  be  argued,  at  least  over  the  early 
period  following  the  events.  The  Borucki  et  al.  work  showed 
that  the  new  chemistry  has  little  effect  on  the  results  for  MO 

x 

perturbations  at  the  altitudes  involved;  substantial  unexplained 
discrepancies,  especially  at  low  altitudes,  were  noted  (R.C. 
Whitten,  private  communication,  1978). 

Injections  of  NGx  from  nuclear  weapons  tests  in  the  1958- 
62  period  have  also  been  studied  by  several  groups,  following 
a suggestion  by  Foley  and  Ruderman  in  1972  (see  Foley  and 
Ruderman,  1972,  1973;  Bauer  and  Gilmore,  1975;  COMESA,  1975; 

NAS,  1975)-  Only  1-D  model  studies  have  been  performed,  al- 
though there  is  evidence  (see  Appendix  B)  that  the  season  and 
latitude  at  which  the  explosion  takes  place  is  important.  Model 
studies  showed  that  hemispheric  average  depletions  of  the  order 
of  6 to  8 percent  should  have  occurred,  but  whether  there  was 
any  evidence  that  it  did  occur  can  be  debated  (see  Bauer  and 
Gilmore,  1975;  NAS,  1975,  Angell  and  Korshover,  1977)-  The 
measurement  by  Howard  and  Evenson  (1977)  of  a much  faster  rate 
for  the  HO,-,  + NO  reaction  than  thought  earlier,  led  to  an  im- 
mediate possible  explanation  for  the  apparent  discrepancy,  be- 
cause it  implied  a much  lower  sensitivity  of  the  ozone  column 
to  NO  injections.  Calculations  by  LLL  (Chang  et  al.,  1978) 
have  indeed  shown  a lesser  sensitivity,  with  reduction  in  esti- 
mated effect  to  about  2-4  percent,  depending  on  (unknown) 

stabilization  altitudes  for  an  NO  cloud.  The  2 percent  lower 

x 

value  would  not  be  inconsistent  with  observations  (Angell  and 
Korshover,  1977);  within  the  uncertainties,  then,  the  LLL  authors 
argue  that  there  are  not  necessarily  any  incompatibilities  be- 
tween theory  and  observation. 

It  should  be  noted  at  this  point  that  a major  point  made 
in  NAS  (1975)  is  not  discussed  in  Chang  et  al.  (1978),  and  re- 
lates to  solar  cycle  effects  on  ozone;  a missing  "increase"  is 
noted  by  NAS.  Indeed,  substantial  variations  in  total  ozone 
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might  be  expected  if  the  solar  ultraviolet  and  the  ozone  column 
variation  with  time  is  as  large  as  may  be  the  case  (global  ozone 
changes  of  7 percent  or  more,  solar  maximum  to  solar  minimum, 
Callis  and  Nealy,  1978).  This  variation,  which  is  quite  un- 
certain, represents  the  largest  known  potential  source  of  ozone 
variability  during  the  1955-1975  time  period  (see  Bauer,  1978). 
Obviously,  if  an  ozone  record  is  to  be  examined  for  comparison 
with  models,  all  such  variations  with  time  should  be  considered. 

Finally,  a perhaps  more  significant  discrepancy  lies  in 
the  ozone  records  versus  the  model  results  in  terms  of  the  halo- 
carbon  problem.  Thus,  Crutzen  et  al.  (1978)  compute  that  ozone 
decreases  of  0.65  percent  per  year  should  have  been  occurring 
at  about  40  km  due  to  transport  of  Cl  species  from  halocarbons 
released  at  the  surface.  However,  during  this  same  period, 
based  on  Umkehr  records,  Angell  and  Korshover  (1977)  show  a 
general  upward  trend,  as  does  Pittock  (1977  ) j based  on  balloon 
flights  over  Aspendale,  Australia.  Taken  at  face  value,  these 
observations  would  represent  a serious  challenge  to  the  C1X 
theories.  Questions  about  the  Umkehr  data,  however,  and  other 
possible  reasons  for  observed  changes,  preclude  any  such  con- 
c lusion . 

These  various  model  validation  questions  are  all  evidently 
confounded  by  a variety  of  time-varying  phenomena  (see  Bauer, 
1978).  A critical  requirement  would  seem  to  be  the  difficult 
one  of  establishing  solar  ultraviolet  variability  with  time, 
preferably  as  a function  of  sunspot  number  or  some  other  easily 
identifiable  variable. 
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4.  CLIMATIC  EFFECTS 


4.1  INTRODUCTION 

Little  specific  research  attention  has  been  given  to  the 
potential  climatic  effects  of  aircraft  exhaust  in  the  period 
since  the  cut-off  date  (December  1976)  for  input  to  our  previous 
report  (Oliver  et  al.,  1977).  The  reasons  for  this,  of  course, 
relate  to  the  complexity  of  the  climate  problem  in  general,  to 
the  lack  of  an  immediate  or  evident  threat  from  aircraft,  and 
to  the  finite  nature  of  resources  with  which  to  study  aircraft 
effects;  in  addition,  from  a national  standpoint,  interest  has 
moved  into  problems  of  more  immediate  policy  relevance,  princi- 
pally the  potential  effects  of  halocarbons  and  of  increases  in 
atmospheric  carbon  dioxide.  Even  for  these  however,  the  actual 
level  of  climatic  modeling  work  accomplished  has  been  modest 
(although  further  work  is  certainly  in  the  planning  stage)  and 
quite  cautiously  reported.  (See  e.g.,  NAS,  1976;  NAS,  1977.) 

In  view  of  this  general  lack  of  new  work  specific  to  air- 
craft climatic  effects,  there  would  seem  to  be  little  to  discuss. 
This  is  not  quite  the  case,  as  climate  modeling  and  climate 
research,  wherever  done,  are  applicable  in  general  to  a variety 
of  problems.  Furthermore,  it  seems  appropriate  to  devote  a few 
paragraphs  to  the  ensemble  of  current  climatic  threats,  and  to 
note  the  common  and  uncommon  features  of  the  aircraft  threat 
vis-a-vis  the  others.  In  addition,  there  have  %been  some  studies 
of  rather  direct  (if  not  necessarily  quantitative)  relevance, 
and  some  apparent  trends  in  the  treatment  of  the  general  problem 
seem  worthy  of  note.  These  follow  in  the  order  just  mentioned. 
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4.2  THE  VARIOUS  (AND  CHANGING)  "THREATS"  TO  OZONE  AND  CLIMATE 


A number  of  the  various  "threats"  postulated  (ir.  recent 
years)  to  climate  and  ozone  are  given  in  Table  4-1,  along  with 
some  characteristics  as  now  perceived.  Note  that  any  material 
which  affects  ozone,  either  in  total  column  cr  ir;  altitude 
latitudinal  distribution,  also  affects  climate;  similarly,  any 
factor  causing  a change  in  climate--or  a change  in  stratospheric 
temperatures — will  affect  this  ozone  column. 

The  aircraft  threat  (NO  ) is  shown  in  Table  4-1  as  leading 
to  an  increase  in  ozone,  a result  shown  by  at  least  some  models 
(see  Section  3),  even  for  aircraft  operating  at  20  km.  For 
from  fertilizer  or  power  plants  (directly  from  combustion,  or 
indirectly  from  fixed  nitrogen) , the  table  shows  a decrease  in 
ozone,  but  the  chemistry  as  used  in  the  calculations  on  which 
this  result  is  based  is  in  need  of  revision;  some  models  do  shew 
an  increase.  From  a climate  (mean  temperature)  standpoint,  both 
aircraft  effluents  and  N~0  increases  may  well  result  in  warming, 
although  in  both  cases  a "netting  out"  of  cooling  and  warming 
effects  involved:  for  aircraft,  see  Oliver  et  al.,  1977;  for 
N^O,  the  tropospheric  greenhouse  effect  of  added  No0,  which 
absorbs  at  7.78  pm  in  a "partial  window"  and  also  at  17  pm  and 
4.5  pm,  all  in  the  infrared,  probably  more  than  compensates  ror 
the  cooling  which  would  be  expected  if  stratospheric  ozone  is 
depleted.  No  model  calculations  including  both  these  effects 
are  known  to  this  author;  however,  Yung  et  al.,  1976,  indicate 
0.5  to  0.75  K warming  with  a doubling  of  N^O,  whereas  Crutzer. 
and  Howard  ( 1978),  indicate  at  most  a 10  percent  reduction  ir; 
stratospheric  ozone  for  a doubling  of  N,0.  A 10  percent  re- 
duction of  stratospheric  ozone  may  involve  a cooling  of  the 
order  of  0.1  K (see  Oliver  et  al.,  1977,  pp.  4-10,  4-11) 
ignoring  ozone  altitude  redistribution  effects,  NO,  changes. 
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TABLE  4-1.  A COMPARISON  OF  ANTHROPOGENIC  EMISSIONS  AS  THREATS  TO  OZONE  AND  CLIMATE 


Release  depends  on  reprocessing  procedures. 

Multiple  effects.  Cooling  due  to  stratospheric  aerosols  and  if  ozone  is  depleted;  warning  due  to  H?0,  contrails,  NG~,  and  0. 
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the  C0o  problem;  however,  a CO^  budget  question  (Kerr,  197?; 
Woodwell,  1978)  exists  which  must  (or  should)  first  be  resolved. 
One  power  plant  exhaust  ingredient,  CCS,  has  no  known  tropo- 
spheric sink,  is  not  removed  by  lime  scrubbers,  and  will  pre- 
sumably increase  the  stratospheric  aerosol  layer  with  some 
compensating  (cooling)  effect.  The  fly  ash  and  SO^-derived 
aerosols  would  also  have  some  effect — whether  warming  or  cooling 
is  perhaps  controversial — but  these  will  likely  be  removed  from 
the  stack  gas  by  scrubbers  and  precipitators.  It  is  of  interest, 
as  pointed  out  by  Robinson  (1977),  that  the  net  loss  in  solar 
radiation  due  to  sulfuric  acid  aerosols  (assuming  2.5  percent 
sulfur  in  the  coal  and  no  scrubbers)  could  exceed  the  heat 
released  by  combustion,  arguing  that  these  effects  should  not 
be  considered  in  isolation.  An  increase  in  CO^  will  cause  a 
coding  of  the  stratosphere,  with  a reduction  in  thermal  reaction 
rates  of  ozone-destruction  cycles,  resulting  in  a slight  increase 
in  ozone  (ceteris  paribus!). 

Nuclear  power  plants,  whatever  their  other  aspects  might 
be,  are  not  normally  considered  to  be  a threat  to  ozone  or  to 
climate,  except  for  waste  heat  effects.  However,  as  pointed  out 
by  Boeck  (1976)  such  plants  do  generate,  and  waste  fuel  repro- 
cessing could  release,  Kr-85  which  is  a g emitter  with  a 10.8 
year  half  life.  Kr-85  would  affect  atmospheric  electricity, 
with  unknown,  but  possibly  important,  effects  (see  Marksor.,  1973) 
also,  some  NO  would  be  formed  as  the  principal  8 (O.2U9  MEV)  is 
slowed  in  the  air.  Both  climate  and  ozone  could  be  affected  by 
the  various  processes  involved.  With  care,  of  course,  Kr-85 
could  be  largely  contained  rather  than  released  to  the  atmos- 
phere . 

These  various  threats  have  common  and  unique  features. 

All,  of  course,  involve  much  the  same  chemistry  and  physics, 
and  all  could  impact  similarly  on  the  biosphere  if,  e.g.,  there 
is  a general  warming,  or  a destruction  of  stratospheric  ozone, 


etc.  But  there  are  important  differences  both  in  a geophysical 
and  policy  sense.  Thus,  aircraft  sources  are  concentrated  in 
mid-latitudes  with  effluents  that  have  a short  life  in  the  trop 
sphere;  halogen  and  N^O  sources  involve  surface  releases  of 
materials  with  long  tropospheric  lives  and  with  stratospheric 
destruction.  Carbon  dioxide  is  inert  in  the  atmosphere  but  is 
absorbed  by  (growing)  plants  and  the  oceans.  These  differences 
mean  that  effects  will  differ  relatively  with  latitude  and 
response  and  recovery  times  will  also  differ.  Prom  a policy 
point  of  view,  it  is  important  to  know  whether  the  atmosphere 
would  recover  in  5 years  or  100  years  if  an  error  were  made. 
Aircraft,  of  course,  have  long  life  times  (approximately  20 
years)  and  policy  decisions  are  difficult  to  reverse;  should  it 
prove  necessary,  the  atmosphere  would  recover  if  the  polluting 
source  were  discontinued,  much  more  rapidly  in  the  case  of  air- 
craft than  in  the  case,  e.g.,  of  halocarbons. 

4.3  CLIMATE  AND  CLIMATE  MODELING  - SOME  DEVELOPMENTS  OF 
INTEREST 

4.3.1  Scope 

Climate  modeling  and  the  study  of  climates  are  fields  in 
their  own  right,  and  no  attempt  will  be  made  here  to  assess  the 
real  or  potential  significance  in  evaluating  aircraft  effects 
of  the  various  publications  in  this  field.  A few  papers  have 
touched  on  matters  of  evident  interest,  and  some  comments  on 
these  matters  follow. 

4.3.2  Predictability  of  Climate  Change 

In  a particularly  sobering  paper,  Robinson  (1978)  examines 
the  prediction  equations  for  complicated  hydrodynamic  flows  "in 
a critical  and  didactic  manner."  His  conclusions  are  obviously 
pertinent  to  the  entire  modeling  problem,  and  seem  to  be  best 
presented  by  direct  quotation: 

"The  only  logically  valid  prediction  with  the  types  of 
equations  now  used  in  meteorological  practice  is  one  of  no 
change . 


"Predictive  equations  ostensibly  based  on  dynamics 
are  empirical  and  must  be  validated  empirically. 

"...Other  predictive  schemes  merit  a priori  the  same 
careful  consideration  as  those  ostensibly  based  on  dy- 
namics . 

"Empiricism,  in  itself,  does  not  limit  the  value  of  a 
predictive  scheme  but  it  does  impose  stringent  require- 
ments for  validation.  The  implications  of  this  are 
particularly  disturbing  in  the  context  of  understanding 
and  predicting  climate  change,  where  the  historical 
record  is  insufficient  to  provide  useful  validating 
statistics,  unless,  as  seems  unlikely,  it  can  be  greatly 
improved  in  detail.  Validation  of  an  empirical  method 
of  climate  prediction  must  proceed  simultaneously  with 
its  application.  Several  prediction  intervals  must 
elapse  before  the  performance  of  the  method  can  be 
assessed;  for  example,  skill  in  predicting  decadal 
averages  might  begin  to  be  quantified  after  a century  of 
application.  There  is  no  reason  to  expect  a more  rapid 
resolution;  even  if  the  relevant  atmospheric  statistics 
are  stationary,  the  integral  time  scale  is  certainly 
more  than  a century." 

Robinson's  critiques  must  obviously  be  kept  in  mind  when 
climate  change  estimates  are  offered.  As  Robinson  indicates, 
however,  the  third  statement  above  represents  the  positive 
implication  of  his  arguments.  In  effect  (and  Robinson  does  not 
make  the  statement),  it  might  be  argued  that  relatively  simple 
formulations  compared  to  large  general  circulation  models  (such 
as  1-D  or  2-D  models)  may  provide  useful,  if  not  necessarily 
quantitatively  reliable,  information;  in  a sense,  all  such  pre- 
dictive techniques  must  await  validation  and  any  prediction  at 
the  present  state  of  understanding  (in  this  author's  opinion) 
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should  be  considered  primarily  as  a measure  of  risk  rather  than 
as  a quantitative  prediction.  Other,  perhaps  more  pragmatic, 
aspects  and  results,  follow. 

4.3.3  Stratospheri c Aerosol  Effects 

The  principal  climatic  effort  in  CIAP  was  in  studies  of  the 
effects  of  aerosols  resulting  from  fuel  sulfur  in  SSTs , (see 
Oliver  et  al . , 1977).  The  conclusion  reached  in  CIAP  was  that 
such  aerosols  will  cause  a cooling  effect  at  the  surface.  The 
water-vapor-induced  warming  was  argued  in  CIAP  to  be  smaller 
"han  the  cooling  effect  of  the  aerosols;  a reexamination  ir.  our 
previous  report  (Oliver  et  al . , 1977)  however,  recognizing  sedi- 
mentation effects  and  attempting  internal  consistency  in  the 
treatments,  suggested  that  the  warming  effects  were  probably 
larger'.  No  studies  (of  which  we  are  aware)  have  been  undertaken 
to  resolve  this  disagreement.  Studies  have,  however,  been 
reported  which  pertain  specifically  to  aerosol  effects.  It  is 
of  interest  that  there  now  seems  to  be  agreement  that  strato- 
spheric aerosols  do  cause  some  cooling  effect.  The  quantifi- 
cation of  this  effect,  however,  for  a specified  fleet,  remains 
highly  uncertain,  because  both  steady-state  burdens  and  effects 
per  unit  mass  involve  many  unknowns. 

A.  paper  of  specific  and  general  climatological  interest  is 
that  of  Mass  and  Schneider  (1977).  These  authors  in  a follow-up 
study  to  an  earlier  paper  (Schneider  and  Mass,  1975 ) report  on 
statistical  evidence  as  to  the  influence  of  sunspots  and  vol- 
canic dust  on  long-term  temperature  records.  The  authors 
conclude  that  ’'heir  earlier  (1975)  conclusions,  which  showed 
a dominant  effect  of  sunspots,  does  not  hold  up  or,  critical 
examination;  they  do  find,  however,  that  a volcanic  signal  can 
be  weakly  detected  in  the  temperature  records.  This  result 
thus  minimizes  questions  associated  with  numerous  contentions 
or  analyses  by  other  authors  who  have  studied  effects  of  vol- 
canic dust,  and  who  have  largely  or  totally  ignored  sunspot 


numbers  (e.g.,  Lamb,  1972;  Mitchell,  1961,  1970;  Oliver,  1976; 
Bradley  and  England,  1978).  It  Is  also  of  interest  that 
Livingston  (1978)  has  recently  reported  measurements  which 
indicate  a cooling  of  the  sun's  photosphere  with  increased  sun- 
spot activity,  rather  than  a warming,  as  required;  e.g.,  if  the 
"Little  Ice  Age"  is  to  be  explained  by  the  lack  of  sunspots 
during  the  Maunder  Minimum  (Schneider  and  Mass,  1975). 

A second  paper  requiring  comment  is  that  duo  to  Hamiil, 
Toon,  and  Kiang  (1977).  These  authors  have  examined  in  detail 
the  microphysical  processes  affecting  stratospheric  aerosols, 
specifically  those  composed  of  sulfuric  acid  and  water,  such  as 
would  result  from  S0n  effluent  from  aircraft . The  authors 
consider  the  altitude  region  in  which  nucleation  of  such  parti- 
cles takes  place  and  find  that  it  corresponds  to  that  region  of 
the  stratosphere  where  such  particles  are  generally  found. 
Various  features  of  the  aerosol  layer  can  be  explained  by  the 
analyses.  Of  particular  interest  is  their  observation  that  such 
droplets  are  unstable — i.e.,  they  will  evaporate--at  higher 
altitudes,  above  about  30  km,  and  that  such  particles  will  grow 
with  time  at  lower  altitudes,  depending  or.  the  ambient  H„,S  • 
vapor  pressure.  Computations  of  residence  times  and  injection 
ffficients  for  particulates,  such  as  those  by  Bauer  and 
Gardner  (Appendix  D in  Oliver  et  al . , 1977)  , assuming  such 
particulates  to  be  inert,  thus  cannot  hold  for  sulfuric  aci:; 
aerosols.  Hamiil  et  al.  (1977)  point  out  some  of  the  still 
outstanding  questions  about  their  modeling  work,  among  which  ! r 
the  role  and  importance  of  (NH^)0  SO^  at  various  altitudes. 
Analyses  of  such  particles  are  provided  by  Farlow  et  al.  (!'•  ’?); 
these  authors  report  finding  various  insoluble  materials  ir.  "he 
particles  and  conclude  that  sources  other  than  volcano* . c - 
tribute  significantly  tc  sucl  aerosols . Fhe  aeros  Ls , in 

neral,  wer<  f und  t be  mixture:  of  solids  and  Liqui  h r ther 
than  pure  sulfuric  acid. 


A 1-D  aerosol  model  based  on  this  and  additional  analysis 
has  been  developed  by  the  same  group  (Turco  et  al.,  1976). 

4.4  A SOVIET  VIEW 

Budyko  and  Karol  (1976)  in  a review  of  CIAP  research,  and 
Budyko  (1977)  in  an  article  on  present-day  climatic  changes, 
present  views  on  climate  modification  which  differ  from  that 
presented  in  CIAP  and  elsewhere.  Since  stratospheric  pollution 
problems  are  global  in  nature,  and  agreement  must  be  multi- 
national, it  seems  important  that  these  views  be  recognized. 

In  Budyko  and  Karol  (1976),  the  climatic  impact  of  aerosols 
and  water  vapor  as  presented  in  CIAP  are  reviewed.  A case  is 
drawn  from  CIAP  which  indicated  by  the  year  2000  that  a drop  ir 
the  average  air  temperature  of  0.l6°C  would  result  from  air- 
craft-generated aerosols.  They  note  a cooling  of  about 
0.1 °C  (including  NCu  compensation),  due  to  10  percent  ozone 
reduction,  and  a warming  of  0.1  to  0.2°C  from  water  vapor.  They 
then  state  "Thus,  the  total  effect  of  the  influence  of  the 
change  in  the  mass  of  the  ozone  and  water  vapor  on  thermal 
condition  may  prove  to  be  insignificant  even  with  high  levels 
of  aircraft  propulsion  effluents.”  [Note  that  this  interpre- 
tation differs  from  that  in  CIAP,  where  (p.  74,  Grobecker  et  al . , 
1974)  it  was  argued  that  NO.,  increases  compensated  for  ozone 
decreases,  and  temperature  changes  due  to  ozone  changes  were 
not  treated.] 

After  further  discussion  of  various  other  aspects  of  the 
CIAP  work,  Budyko  and  Karol  (1976)  note  that  the  effects  of 
sulfate  aerosol  obtained  by  CIAP  agree  reasonably  well  with 
their  own  independent  estimates.  They  then  note  the  possibility 
of  a "directed  impact"  on  climate  due  to  the  burning  of  sulfur 
in  the  lower  atmosphere,  and  argue  that  another  circumstance  is 
of  essential  importance  to  an  evaluation  of  the  conclusions  of 
the  'I A:  research . This  circumstance  is  the  expected  increase 
in  temperature  due  to  increases  it.  carbon  dioxide  fr  m combustion 


of  fossil  fuels,  and  they  point  out  that  the  increase  in  expected 
temperature  of  0.5  to  1.0°C  exceeds  the  cooling  due  to  strato- 
spheric aircraft.  This  warming  would  result  in  reduced  precipi- 
tation in  mid-latitudes,  with  "substantial  losses  in  the  economie 
of  many  countries."*  The  stratospheric  aviation  effect  may  then 
prove  to  be,  first,  insignificant  compared  to  other  changes,  and 
second,  the  aerosol  may  diminish  the  CO^-induced  tendency 
towards  warming.  They  point  out  that  the  second  effect  was 
noted  by  them  earlier,  and  give  references. 

In  the  later  paper,  Budyko  (1977)**  argues  for  the  coupling 
of  climate  theory  with  economic  aspects.  He  includes  plots 
which  show  that  as  the  mean  annual  temperature  anomaly  for  the 
northern  hemisphere  (and  for  75-85°N)  increased  (since  1883) 
the  precipitation  in  mid-latitudes  decreased;  decreases  in  in- 
tegrated precipitation  (as  measured  by  the  Caspian  sea  level) 
correlated  with  decreases  in  wheat  yield.  He  then  discusses 
climate  models,  effects  of  atmospheric  transparency  and  changes 
in  cloudiness  and  albedo,  examines  data,  and  concludes  that,  at 
least  for  periods  over  which  the  ice-albedo  feedback  is  small, 
climate  models  satisfactorily  explain  natural  changes  and  are 
thus  capable  of  predicting  anthropogenic  factors  on  climate. 

He  considers  C0o  and  aerosol  (tropospheric)  effects  on  climate, 
and  warns  that  warming  may  result  with  a decrease  in  rainfall 
in  mid-latitudes  (presumably  leading  to  reduced  food  production) . 
He  then  refers  to  his  earlier  works,  and  points  out  that  an 
increase  in  stratospheric  aerosol  could  be  effected  to  counter 
the  effects  of  CO^  warming. 

Budyko  and  Karol  also  criticize  the  CIAP  investigators  for  not 
including  in  their  evaluation  the  increased  rainfall  in  mid- 
latitudes with  global  cooling,  even  though  the  information, 
in  their  view,  was  available  to  the  CIAP  investigators. 

*# 

According  to  Tellus , the  manuscript  was  received  November  19, 
1975  and  in  final  form  June  4,  1976. 
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AN  ALTERNATIVE  LOOK  AT  THE  GENERAL  ELECTRIC  EMISSIONS  DATA 
(Lyon  et  dl.,  1S75;  Lyon,  1975) 

A. 1 BACKGROUND 

ns  discussed  in  the  main  text,  the  issue  of  the  appropriate- 
ness of  probe  sampling  techniques  (particularly  under  afterburn- 
ing conditions)  to  determine  NG  in  jet  engine  exhausts  was 
raised  by  the  work  of'  Davidson  and  Domal  (1973),  who  found, 
using  in-situ  spectroscopic  techniques,  that  probe  techniques 
may  give  severely  low  NO  values.  The  question  led  to  the  ref- 
erenced work  by  GE  investigators  who  studied  1.0  measurements 
by  probe  techniques  over  a wide  range  of  cor.uit:  ,ns  (military 
power  to  full  afterburning)  and  from  zero  to  at  ut  2“  nozzle 
diame  ers  iownstream.  The  and  J 7 ■ sngines  re  used.  The 
GE  work  did  not  employ  spectroscopic  techniques,  but  did  sample 
cor.oustior.  muses  at  high  and  low  temperature  conditions  and 
with  various  degrees  of  dilution.  GE  concluded  that  the  probe 
techniques,  d properly,  give  emission  index  values  indepen- 
dent of  ; osi  In  downstream,  and  hence,  implicitly,  the  probe 
techniques  must  be  accurate. 

The  GE  effort  was  broken  into  several  phases.  The  first 
phase  was  exploratory.  The  second  phase  involved  a great  deal 
of  data  collection,  but  some  difficulty  with  measuring  instru- 
ments was  noted,  and  analytical,  procedures  were  under  develop- 
ment . The  data  did  show  a disconcerting  tendency  to  give  higher 
NO  emiss'on  indices  as  gas  samples  were  taken  from  cooler  (more 
dilate)  regions,  but  material  balance  questions  were  evident  in 
Examples  of  these  effects  are  shown  in  Table  A-l. 
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TABLE  A-l . APPARENT  NOx  EMISSION  INDICES  vs.  DOWNSTREAM  POSITIONS 
J85-5 , MILITARY  POWER,  PRELIMINARY  TEST  RESULTS 


Station,  ft 

Metered  Fuel 

Flow,  pps 

Calculated 

Fuel  Flow,  pps 

E.I.  NO, 

0 

0.683 

0.687 

4.4 

3.75 

0.689 

0.614 

4.7 

7.5 

0.708 

0.654 

4.7 

15.0 

0.719 

0.598 

6.7 

30.0 

0.714 

0.280 

10.3 

Source:  Lyon  et  al . , 1975  (p.  124) 

Similar  trends  were  noted  at  other  conditions  and  for  both 
engines,  although  the  calculated  fuel  flow  and  metered  fuel 
flows  did  not  all  behave  in  the  fashion  shown  in  this  example. 

The  difficulty  was  ascribed  to  measurement  errors,  primarily 
in  CO^  content  at  large  dilutions,  and  in  some  cases  to  back- 
ground contaminant  levels.  A technique  was  thus  developed  in 
which  measured  NO  was  plotted  against  measured  CO-  values,  and 
emission  indices  deduced  from  the  slopes.  This  technique  was 
then  used,  together  with  improved  instrumentation  in  Phase  III. 
Phase  II  involved  measurements  at  four  stations.  Phase  III  in- 
volved measurements  only  at  the  exhaust  plane  and  about  24  noz- 
zle diameters  downstream.  Results  of  Phase  III  work  to  be  ana- 
lyzed here  are  shown  in  Figs.  A-l  to  A-4;  the  corresponding 
detailed  tabular  data  are  given  by  Lyon  (1975)  (Tables  4l,  45, 

49,  and  53  for  military  power  and  Tables  44,  48,  52,  and  56  for 
maximum  afterburning  power).  Tables  A-2  and  A-3  show  the  GE 
summarized  results. 

Figure  A-4  represents  an  extreme  case  in  which  a severely 
polluted  atmosphere  is  implied,  with  background  NCx,  hydrocarbon, 
and  CO-  contents  of  about  2,  5>  and  10  ppmv,  respectively.  The 
plots  indicate  a negative  hydrocarbon  emission  index;  the  entries 
are  dashed  in  Table  A-3. 
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FIGURE  A-2.  Emissions  concentrations  vs.  CO2  for  J85- 5 engine 
30  feet  aft  of  exhaust  nozzle,  at  max.  A/B  power. 
Final  data  by  Far  Plume  Method. 

Source:  Lyon  et  al . , 1975 


TABLE  A-2.  SUMMARY  OF  EMISSIONS  DATA  FROM  FINAL  J85-5  TESTS 


Note  that  the  slopes  of  the  lines  in  these  figures  are 
largely  determined  by  the  richer  (and  hotter)  gas  samples,  the 
data  at  highest  dilutions  tending  to  have  little  weight.  In 
order  to  see  whether  there  may  be  information  of  interest  in 
these  high  dilution  data,  an  alternative  procedure,  which  gives 
greatest  emphasis  to  dilute  gas  samples,  has  been  developed, 
and  is  described  below. 

A. 2 ALTERNATIVE  PROCEDURE 

In  the  approach  that  follows,  no  attempt  is  made  to  deter- 
mine a mass-weighted  emission  index  for  the  entire  stream. 
Rather,  the  question  is  asked  whether  the  various  data  on  Indi- 
vidual air  samples  imply  higher  emission  indices  as  the  parcels 
are  cooled  by  mixing  with  ambient  air. 


Note  that  if  fuel  and  exhaust  gas  compositions  are  known, 
the  NO  emission  index  and  the  air-to-fuel  ratio  for  that 
sample  can  be  calculated  simply  according  to 
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measured  NO  , CO,  and  unburned  hydrocarbon  contents,  in  ppmv 
(corrected,  if  necessary,  for  ambient  levels). 

At  dilute  conditions,  where  CO  and  HC  contents  are  negli- 
gible , 
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or,  more  precisely, 

E . I . NO 

x 

where  correction  is  made  for  background  levels  of  N0x  and  C0„ . 

This  is  the  basis  of  the  GE  procedure.  Unfortunately,  however, 
background  NO  and  C05  levels  were  not  measured  (except  indirectly 
by  noting  intercepts)  so  that  high  dilution  results,  which  may 
contain  the  desired  information  are  subject  to  severe  uncer- 
tainties. To  explore  these  effects,  certain  assumptions  were 
necessary  as  follows: 

1.  The  C0o  content  was  reduced  by  0.032  percent,  the  known 
minimum  level,  or,  in  some  cases,  where  the  instrument 
showed  a slightly  smaller  value  (e.g.,  0.02  percent) 

by  the  minimum  value  measured. 

2.  The  ambient  NO  was  assumed  (a)  to  be  zero,  or  (b)  for 

X 

the  downstream  cases,  a value  was  read  from  a plot  such 
as  Fig.  A-^,  reading  the  NO^  value  at  0.032  percent  CO^. 
(Ambient  N0x  values  on  unpolluted  air  are  < 1 ppbv, 
which  are  negligible  in  this  context.)  Plots  such  as 
Fig.  A-b  were  not  prepared  by  GE  for  the  cases  at  the 
nozzle  exits.  The  NO  levels  at  this  station  are, 
however,  still  normally  far  above  any  normal  background 
level . 

Calculations  using  these  procedures  were  made  for  military 
and  maximum  afterburning  power  for  the  J85  and  J79  engines. 
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A. 3 RESULTS 


Results  for  the  four  cases 
through  A-8.  Solid  points  are 
gible  background  NO  : crosses 

X 


ground  N0x- 

Note  that  this  method  of  plotting  permits  incorporation  of 
samples  taken  at  any  point  in  the  plurne. 

A. 4 DISCUSSION 

In  all  cases,  the  apparent  NO  emission  index  increases  at 

X 

high  dilution.  The  effect  takes  place  at  much  lower  dilutions 
at  maximum  afterburning  than  at  military  power. 

A correction  for  background  N0x  has,  as  would  be  expected, 
a powerful  effect  on  the  apparent  emission  index  at  high  dilu- 
tions. The  correct  value  of  background  NO  to  use  is  unknown, 
and  may  have  changed  during  the  course  of  the  run.  Thus,  there 
is  great  uncertainty  in  the  derived  value  of  the  NO  emission 

X 

index . 

In  the  afterburning  cases,  an  interesting  region  exists  in 
which  the  NO  apparent  emission  index  decreases  with  air-to- 
fuel  ratio.  This  may  imply  destruction  of  NO  with  time  by 

X 

reactions  with  CO,  H^,  or  hydrocarbons.  After  passing  through 
a minimum,  more  NO  may  be  generated  by  heating  of  excess  air. 
Probe  surface  reactions  may  also  enter. 

For  the  militarv  power  cases,  the  rise  in  apparent  emission 
index  at  air-to-fuel  ratios  greater  than  about  1000  to  1 is 
puzzling.  The  temperature  rise  at  dilutions  of  this  order  is 
small:  at  18,000  Btu/lb  of  fuel,  the  temperature  rise  at  1000 

to  1 is  about  72°F.  Surface  catalytic  effects  would  seem  im- 
probable. The  temperature  maintained  inside  the  line  is  higher 
(about  300°F)  but  still  doesn't  seem  high  enough  to  initiate 
surface  catalytic  reactions  (such  as  the  water-gas  shift  reaction, 
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followed  by  NO  reduction).  [It  would  be  of  interest,  if  it 
has  not  been  done,  to  run  known  mixtures  of  NO  (in  moist  air 
containing  CO  and  hydrocarbons)  through  a probe  tip  and  several 
hundred  feet  of  line  as  used  in  the  GE  tests.] 

The  GE  report,  in  general,  shows  slightly  higher  emission 
indices  for  the  downstream  location  than  for  the  nozzle  exit. 

In  the  J79  maximum  afterburning  case,  the  value  is  lower,  but 
the  method  of  plotting  shown  in  Fig.  A-8  would  not  lead  to  the 
same  conclusion. 

In  early  work  at  AEDC  [Davidson  and  Domal  (197  3)  ]>  the  error 
in  using  probe  techniques  was  greater  with  afterburning  than 
without.  It  is  debatable  whether  the  results  plotted  here  sup- 
port this  observation,  due  to  uncertainties  about  the  NO  cor- 
rection. (Note  suppressed  zeroes  in  Figs.  A-6  and  A-8.  Note 
also  that  a correction  of  2.2  pomv  in  Fig.  A-8  appears  to  be 
excessive . ) 

A. 5 FINAL  COMMENTS 

Before  any  interpretation  of  NO^  emissions  data  can  be 
accepted  which  shows  large  increases  in  emission  index  with 
dilution,  the  mechanisms  involved  would  need  to  be  established. 
The  work  reported  by  Benson  and  Samuelsen  (1977),  e.g.,  might 
suggest  N0x  reduction  on  the  probe  tip  by  unsaturated  hydro- 
carbons. However,  if  this  mechanism  were  operating,  one  would 
expect  low  N0x  emission  indices  in  samples  with  large  hydro- 
carbon contents.  A brief  perusal  of  the  afterburning  data 
(Table  52,  e.g.,  Lyon,  1975)  shows  no  such  effect.  Furthermore, 
the  good  linearity  of  N0x  vs.  CO,,  found  in  the  GE  plots,  and  the 
reasonable  agreement  (by  their  methods)  between  nozzle  exits 
and  30  (J85)  or  60  (J79)  ft  downstream  stations,  would  seem  to 
preclude  any  massive  effect  of  probe  tips,  etc. 
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APPENDIX  B 

ON  THE  USE  OF  Zr-95  DATA  FROM  CHINESE  ATMOSPHERIC 
THERMONUCLEAR  EXPLOSIONS  TO  STUDY  STRATOSPHERIC 
TRANSPORT  IN  A ON E - D I ME  NS  I ON AL  PARAMETERIZATION 

B.l  INTRODUCTION 

There  has  been  considerable  discussion  [Grobecker  e4:  ai . 

( 197 A ) , COMESA  (1975),  NAS  (1975  ),  NAS  (1976),  Oliver  et  al. 
(1977)*]  on  the  transport  of  materials  injected  into  the  lower 
stratosphere.  Particular  questions  arise  for  the  only  currently 
existing  SSTs,  Concorde  and  Tu-1'44,  which,  in  airline  service, 
will  inject  most  of  their  exhaust  products  in  the  15-18  krr. 
altitude,  40-60°N  latitude  region.  The  one-dimensicnal  para- 
meterizations  of  this  motion  by  different  workers  differ  sub- 
stantially . 

From  the  standpoint  of  altitude  and  latitude  of  injection, 
the  single  set  of  data  which  are  most  directly  comparable  to 
stratospheric  aircraft  flight,  in  particular,  the  Concorde  and 
Tu-l44  SSTs,  are  the  Zr-95  data  from  the  five  Chinese  3-Fit 
thermonuclear  explosions  that  injected  debris  at  18  km  alti- 
tude, 40°N,  90°E  - see  Telegadas  (1974  and  1976).  Here  we 
examine  the  integrated  stratospheric  burden  of  the  radioactive 
isotope  Zirconium-95  (Zr-95;  65  day  half-life)  in  terms  of  a 

If 

We  should  like  to  thank  H.S.  Johnston,  D.M.  Hunten,  L.  Machta, 
J.D.  Mahlman,  and  K.  Telegadas  for  helpful  comments.  The 
present  work  supersedes  that  reported  in  Appendix  C of  Bauer 
et  al.  (1977),  which  was  based  entirely  on  Model  I (see  Sec- 
tion B. 3* 1) • 
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one-dimensional  model  which  uses  an  effective  mean  eddy  diffu- 
sivity  K (assumed  to  be  constant  above  the  tropopause)  and  an 
effective  mean  tropopause  height  zq  (which  characterizes  an 
effective  height  of  injection  above  the  local  tropopause)  as 
parameters  to  characterize  atmospheric  transport  of  precipita- 
tion- scavengeable  materials. 

B . 2 EXPERIMENTAL  DATA 

The  U.S.  ERDA  (now  DOE)  "Airstream"  program  has  involved 
sampling  radioactive  debris  from  aircraft  flying  at  altitudes 
of  15  to  20  km  along  a flight  path  following  generally  along 
the  West  Coast  of  North  and  South  America,  from  72°N  down  to 
50°S  (more  recently  to  10°S  only).  These  flights,  at  regular 
intervals  of  4-6  months,  are  supplemented  by  occasional  balloon 
probes  which  are  flown  mainly  from  San  Angelo,  Texas  (32°N, 
101°W)  and  which  go  to  higher  altitudes  (27—35  km). 

The  Chinese  nuclear  explosions  are  conducted  at  Lop  Nor 
( 40°N , 90°E),  which  is  roughly  at  the  opposite  meridian  of 
longitude  to  the  Airstream  corridor;  in  general,  East-West 
transport  is  considered  to  be  fast  compared  to  North-South 
transport,  so  that  the  Airstream  profiles  of  radioactive  debris 
should  provide  an  effective  zonal  average  profile  of  a tracer 
of  atmospheric  motions  some  weeks  after  the  detonation.  This 
assumption,  made  by  Telegadas  in  deriving  an  estimate  of  the 
vertical,  zonal  averaged  distribution  of  radioactivity  with 
height,  introduces  some  undertainty  into  the  data,  especially 
at  early  times . * 

K 

Thus,  J.F.  Louis  (private  communication,  May  1976)  who 
analyzed  the  1967-1970  data  with  a two-dimensional  model, 
treated  the  effective  latitude  of  injection  of  the  1967 
debris  as  55°N  rather  than  40°N;  evidently,  the  cloud  from 
that  detonation  was  caught  in  the  Central  Asian  summer  high- 
pressure  cell  and  thus  took  a relatively  long  time  before  it 
was  zonally  spread. 


Our  primary  concern  in  this  work  is  with  the  rate  of  removal 
of  radioactive  debris  from  the  lowest  stratosphere  on  a global 
basis,  although  the  height  distribution  of  radioactivity  in 
the  stratosphere  is  of  some  interest. 

Among  other  parameters  (such  as  altitude  and  season),  this 
removal  rate  will  vary  with  the  latitude  of  injection,  so  that 
the  analysis  that  follows  applies  specifically  only  to  the 
latitude  under  study;  i.e.,  40°  or,  generally,  mid-latitudes. 

The  data  are  described  in  terms  of  altitude/latitude  pro- 
files presented  every  3 months  in  terms  of  a stratospheric 
burden  of  the  isotope  Zr-95 > decay-corrected  to  the  time  of 
detonation.  The  integrated  stratospheric  burden  of  Zr-95  for 
the  various  Chinese  3-Mt  thermonuclear  explosions  computed  by 
Teiegadas  (197^,  1976)  are  shown  in  Fig.  B-l.  The  numerical 
dab.'-  *'rom  which  Fig.  B-l  was  prepared  are  the  basic  experimental 
data  analysed  here.  Mote  that  for  a summer  injection  there  is 
no  significant  decrease  in  stratospheric  burden  until  midwinter, 
while  with  a winter  injection  of  tracer,  the  stratospheric  bur- 
den decreases  steadily. 

Figure  B-2  shows  the  vertical  Zr-95  activity  distribution 
observed  after  the  various  Chinese  nuclear  tests.  While  the 
data  are  subject  to  major  uncertainties  because  the  Airstream 
sampling  took  place  so  far  from  the  latitude  of  injection  and 
ften  quite  a long  time  after  the  actual  injection,  all  the 
tint  a suggest  an  injection  height  of  the  debris  in  the  17-19  km 
altitude  range.  The  low  value  observed  some  16  weeks  after  the 
October  1970  test  is  evidently  due  to  large-scale  cloud  motion 
•it  some  time  well  after  the  injection  (cf.  the  distribution 
estimated  one  week  after  that  same  test). 

In  a one-dimensional  model,  one  crucial  parameter  for  a 
quantitative  interpretation  of  transport  out  of  the  strato- 
s:  here  is  (z.-z  ),  the  height  cf  injection  above  the  local 
tropopause.  Inspection  of  Fig.  B-2  shows  that  the  present 
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FIGURE  B-l.  Stratospheric  Zr-95  burden  (decay-corrected) 
following  the  Chinese  nuclear  tests.  Ob- 
served burdens  are  indicated  by  the  short 
horizontal  lines. 

Source:  Telegadas,  1974,  1976 
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results  are  consistent  with  injection  heights  z^  of  17-19  km 
and  with  local  tropopause  heights  zq  of  10-12.5  km.  We  describe 
the  data  in  terms  of  a mean  injection  height  = 18  km  and  a 

mean  tropopause  height  zq  = 11  km,  corresponding  to  a net  in- 
jection height  (z-^-z  ) = 7 km  above  the  local  tropopause.  [The 
tropopause  height  over  Lop  Nor  (40°N,  90°E)  is  highly  variable, 
ranging  from  10  km  in  winter  up  to  14-15  km  in  summer — see, 
e.g.,  Crutcher  and  Davies  (1969)-] 

As  an  extreme  in  variability  of  the  effective  stratospheric 
injection  height,  we  also  consider  zq  = 14  km  which  gives 
(z^-zq)  = 4 km.  The  quantity  (z^-zq)  will  tend  to  be  preserved 
in  North-South  transport  since--while  the  local  tropopause 
height  changes  significantly  with  latitude — the  lines  of  con- 
stant mixing  ratio  for  stratospheric  injectants  tend  to  remain 
parallel  to  the  local  tropopause,  suggesting  North-South  trans- 
port parallel  to  the  sloping  tropopause  [e.g.,  NAS  (1975), 
p.  146;  Machta  and  Telegadas  (1973)-] 

B . 3 ANALYSIS  OF  THE  DATA 
B.3.1  Mathematical  Formulism 

For  a chemically  inert  tracer  of  mixing  ratio  f(z,t)  which 
is  injected  at  time  t = 0 at  some  altitude  z = z,  in  the  strato- 
sphere, the  transport  equation  may  be  written  in  terms  of  an 
effective  eddy  diffusion  coefficient  K(z)  as 

(3/3z)[n(z)K(z)3f/3z]  = n(z)  3f/3t  . (B.l) 

We  treat  the  stratosphere  as  a region  of  constant  and  rela- 
tively small  eddy  diffusivity  coefficient  K,  characterized  by 
a tropopause  at  height  z , below  which  transport  is  very  much 
faster;  the  eddy  diffusivity  profile  is  taken  as 
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K,  constant,  in  the  stratosphere;  i.e.,  for  z > 


(B.2) 


K.  , constant,  in  the  troposphere;  i.e.,  for  z < z 


where,  typically,  K^.  ~ 10  K.  This  simple  K-profile  was  used  by 
Machta  (197*1)  and  by  Crutzen  (197*0;  it  is  used  here  because  it 
permits  transport  out  of  the  lowest  stratosphere  to  be  described 
in  terms  of  only  two  parameters,  namely  K and  z . (The  value  of 
is  typically  not  critical,  as  long  as  it  is  sufficiently 

is 

greater  than  K.) 

Most  current  models  for  K(z)  show  a monotonic  increase  with 
height  above  some  minimum  in  the  10-  to  20-km  altitude  range  that 
is  associated  with  the  tropopause;  Lindzen  (1971)  points  out  that 
this  is  a general  characteristic  associated  with  the  decrease  in 
density  with  increasing  height . The  present  model  does  not  show 
this  increase  in  K at  higher  altitudes.  However,  any  model  sat- 
isfying Lindzen' s criterion  requires  at  least  three  parameters 
in  the  region  of  interest  here  [perhaps  z , the  lowest  value  of 
K above  zq,  and  the  slope  dK/dz  or  d(log  K)dz  above  z ],  and  we 
prefer  to  analyze  the  Zr-95  data  in  terms  of  this  simplest  pos- 
sible model  before  going  to  the  additional  complexity  of  a three- 
parameter  model  (which  we  plan  to  do  in  future  work). 

We  also  assume  an  isothermal  stratosphere,  so  that  in  the 
region  of  interest 


, , -(z-z  )/H 

n ( z ; = n e o 
o 


(B . 3) 


where  n(z)  = total  particle  number  density  at  height  z,  n = 
n(zQ),  and  the  "scale  height”  H = kT/Mg  is  taken  as  6.3  km  (cor- 
responding to  T = 216  K) . In  the  region  from  10  km  to  20  km, 
which  is  of  direct  interest  here,  the  assumption  of  constant 
temperature  is  quite  well  satisfied. 


i 


With  the  assumption  of  Eq . (B.3),  the  differential  equation 
(B.l)  becomes 

( 3/3z  ) {K ( z ) 9f/9 z } - { K ( z ) /K } 3 f / 3 z = 9f/3t  . (B.H) 

Clearly,  this  equation  simplifies  somewhat  further  with  the 
assumption  of  Eq . (B.2)  for  the  eddy  diffusivity  profile,  K(z). 

We  initially  ignore  sedimentation  of  the  aerosols  carrying 
the  Zr-95  tracer  (postponing  discussion  to  Section  B.3. 3 below), 
and  now  we  wish  to  solve  the  diffusion  Eq . (B.4)  with  the  con- 
dition (B.2)  for  K(z),  for  the  following  initial  condition 


f (2,0)  * q 5 (z  - z,  ) 


(B.5) 


which  corresponds  to  instantaneous  point  injection  at  height 
z1 ; the  initial  injection  is  not  precisely  a delta-function 
injection,  so  there  is  a slight  time  correction,  as  discussed 
in  Section  B . 3 • 2 . 


Main 
posphere 
mations , 
namely , 


emphasis  is  clearly  on  the  stratosphere;  in  the  tro- 
i.e.,  for  z < z , we  have  made  three  different  apprc 
of  successively  increasing  accuracy  and  complexity, 


• "Infinite  Model"  I:  We  solve  the  transport  equation 
(B.l)  for  constant  K (equal  to  the  stratospheric  value) 
for  all  z.  This  is  a very  simple  model,  but  does  not 
describe  the  situation  in  the  troposphere  or  the  lower 
boundary  at  all  accurately. 

• "Semi-infinite  Model"  S:  Once  again,  we  solve  the  trans- 
port equation  (B.l)  for  constant  K in  both  the  strato- 
sphere and  the  troposphere,  but  now  we  assume  an  absorb- 
ing surface  at  some  height  z.,  (which  may  be  taken  as  on 
the  ground,  at  the  tropopause,  or  at  some  intermediate 
altitude,  possibly  corresponding  to  the  altitude  of  rain 
cut ) . 
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• "Two-step  Model"  T:  We  assume  different  rates  of  trans- 
port in  the  stratosphere  and  troposphere  [i.e.,  f K 
in  Eq.  (B.2)]  and  zero  mixing  ratio  at  the  ground. 

The  solution  of  Eq . (B.l)  with  the  initial  condition  (B.5) 
for  the  infinite  model  I can  be  shown  to  be — see,  e.g.,  Martin 
(1976)  — 

f-j-U.t)  = q(4TTKt)-1/2exp  j -(Kt/4)  (1/H)  - (z-s^'/KtJ  " J.  . (B.6) 

With  the  more  satisfactory  treatment  of  the  lower  boundary 
by  model  S,  it  can  be  shown  [Feller  (1968),  Martin  (1976)]  that 
the  solution  for  z > z0  is 


' r , 

T ^ • 


a - 


exp  [-(z-z2 ) (z, -z0 )/Kt ] } 


(B.7) 


For  the  "two-step"  model  T,  there  is  an  additional  require- 
ment to  matching  stratospheric  and  tropospheric  solutions  at  the 
tropopause,  z = z_.  It  is  evident  that  f(z,t)  must  be  continuous 
at  the  tropopause,  and  further,  since  there  is  no  source  or  sink 
at  the  tropoDause,  the  net  flux  transoorted  across  z = z must 
be  continuous.  By  integrating  Eq . (B.l)  with  respect  tc  z,  this 
condition  is  seen  to  be 


or 


r,(z)  K(z)3f/3z  = continuous  at  z = z 

o 


(B.  : 


..(Etrar .)  (Etrat. , 

a d I 


j_  =K(Trop-)3f(TroP-)/3z  i 


z=z 


vo. Mb) 


It  may  be  shown  that  the  solution  of  the  two-step  model 
characterized  by  Eq . (B.4),  with  the  initial  condition  of  Eq . 
(B.5),  with  the  matching  condition  (E.8b)  at  the  tropopause  z , 
and  with  a sink  at  the  ground  so  that 


fm(  2*0  , t ) = 0 


rT(z,t)  = & exp  (z-J1)/2K 


/ 2 1 \ 

d<  exp-Kt  < + — W(<fz,2  )W(<,z1 ,z^ ! 

UK*  ‘ ” 


, (B . 10  ) 


where 


W(<,:,s  ) = 


K„-K 

<TK).cos<mz^  + ■ 1 s 1 r.  < ~ z 


. sir.<(z-z  ' + <K  sinx^z  cosk(z-z  ) 

CJ  0 1 O 0 


^ .(-hr)  sinKTzo  + ‘tV^Vc . + <2-Sl 


n x^z  (B.ll) 

- U 


stratospheric  K _ K_ 
tropospheric  K Kf 


0 < d < 


(B.12) 


(B. 13) 


The  function  f^Czjt)  cannot  be  expressed  in  terms  of  tabu- 
lated functions,  but  the  total  stratospheric  burden  of  tracer, 
which  is  the  quantity  that  is  analyzed  here,  can  be  expressed 
as  a single  integral  that  can  be  evaluated  numerically--see 
Eq.  (B.15)  below. 

If 

A derivation  of  this  result,  using  transform  techniques 
appropriate  for  the  boundary  conditions,  has  been  given  by 
W.  Wasylkiwskyj , IDA  Inter-Office  Memorandum  of  1 July  1977, 
and  addendum  of  2 September  1977,  "A  one-dimensional  model 
for  studying  diffusion  of  injectants  in  the  stratosphere." 
Copies  of  this  memorandum  are  available  on  request.  It  may 
be  verified  by  direct  substitution  that  Eq . (B.10)  is  indeed 
a solution  of  the  problem. 


For  comparison  with  the  results  of  Fig.  B-l,  which  gives  the 
total  measured  stratospheric  burden  of  Zr-95  as  a function  of 
time,  we  need  the  total  calculated  stratospheric  burden  of  tracer. 


B,  (t) 


■f- 


(z)  fj(z,t)  dz;  (J=I,S,T) 


(B.14) 


Numerical  calculations  of  B.(t)  have  been  made  using  the 
two-step  model  T of  Eq.  (B.10).  Significant  qualitative  insight* 
can  be  gained  with  the  simple  model  I of  Eq . (B.6)  which  gives 


BI(t ) 


-f- 


-(zT-z  )/H 

(z)  f x ( z , t ) dz  * nQq  e 1 U(Xq)  , (B.15) 


where 


xq  - (1/2) (Kt) 


1/2 


[(1/H)  - (Zl-zo)/Kt)]  ; (B . l6a ) 


11/2(1  + erf  x ) for  x^  < 0 
o o 

(B. 16b ) 

1/2(1  - erf  x ) for  x^  > 0 
o o 


The  overall  relationship  between  time  t,  x , and  U(xQ)  is 
given  in  Table  B-l. 


There  are  quantitative  defects  in  model  I,  which  gives  high 
values  of  K because  the  relatively  fast  transport  in  the 
troposphere  cannot  be  simulated  in  any  other  way  by  this 
model . 
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TABLE  B-l. 


Time  t 


xq  AND  U(x  ) FROM  EQ.  (B.9)  AS  FUNCTIONS  OF  TIME  (FROM  MODEL  I) 


Smal  I 


H(z1-zq)/K 


large,  0 

negative 

1 0.5 


Large 

(Kt)1/2/2H 

( 1/2 ) ( 1 -erf  x )-*0  as  t-*» 


The  physical  meaning  of  Eqs.  (B.15  and  B.16)  is  relatively 
clear-cut.  The  quantity  qn  e*  “1  o'  in  Eq . (B.15)  is  the 
total  number  of  molecules  per  unit  area  injected  at  t = 0 at 
2 = xQ  is  a time-dependent  variable,  and  U(x  ) is  the  frac- 
tion of  injected  material  that  remains  in  the  stratosphere 

after  time  t,  where  x = x (t). 

’ o o 

Prom  Table  B-l  we  see  that  a simple  way  to  determine  K is 
from  the  time  tT.  at  which  U(x  ) or  the  stratospheric  burden 

-L  > "2  O 

in  the  approximation  I has  half  its  initial  value,  because 


t 

u T 


H/K 


(B.17) 


This  result  shows  qualitatively  how  the  effective  stratospheric 
K-value  depends  on  the  decay  time  of  the  stratospheric  burden 
and  also  on  the  effective  injection  height  above  the  tropopause. 
Mote,  in  particular,  that  a given  decay  rate  (l/tT  , ) may  be 
obtained  by  rapid  transport  (K  large)  through  a large  distance 
(z^-z  ),  i.e.,  a low  tropopause;  or  by  slow  transport  (K  small) 
and  a high  tropopause,  i.e.,  slow  transport  through  a short 
distance  only. 

*e  that  the  net  transport  rate  in  the  stratosphere  is  prc- 
• . • K(z)  n(z)3f/3z.  In  comparing  the  infinite  model  I 
••  - tep  model  T,  we  note  that  in  the  lower  stratosphere 
?f-/3z  ! because  of  f^  = 0 at  the  tropopause  while 
' • ■*  early  time  fT(z1)  = f^Cz^)].  Thus,  for 

!•  -he  stratosphere,  model  I gives  significantly 


faster  1 ’anspi  rt  than  does  model  T,  and  conversely,  to  explain 
given  set  of  data  with  model  T requires  a smaller  K-value 
chan  does  model  I. 

For  the  two-step  case  T,  an  expression  for  the  stratospheric 
burden  B^Ct)  in  the  form  of  a single  integral  which  can  readily 
be  evaluated  numerically  is  given  below: 


b„.  ( L ) « e 

T h 


r .u  e"Ksl(K'  f <m2) 

■ 1 X 

JQ  ,...2 


|(-  ^ oy  > It  S,IVo  (Wo  * h®  * Kd3l,*^0COS*(VV 

( , ~ v 7777 

I -^TT  hill*  , * K„l(,GSKrnZ  I dl  Sill  K„,Z 

\ fH  1 o To/  1 o 


(B.18) 


Here,  k„  is  defined  in  Eq . (B.12)  and  d is  defined  in  Eq 

i 

(3  13).  Note  chat  the  value  of  kt  passes  through  zero  as  < 
c-s  from  0 to  ®;  when 


<2  + 


u r. 


(B . 19a) 


we  must  make  the  substitution 


:os<T:o 


sin  <_z. 


C0ShKm2 


isinnic-z. 


(B.l9b) 


in  the  integrand  of  Eq . ^B.l8). 

It  is  of  interest  to  compare  the  solution  f^  of  Eq.  (B.10) 
with  the  much  simpler  solution  fe  of  Eq . (B.7),  bearing  in  mind 


E-15 


P 


( B . 20a ) 


lim  f„  = fQ(z0  = z );  i.e.  , a sink  at  the  trcpause, 
a n 1 ^ d 0 


because  K<<K. 

t 

lim  f_  = fo(z0  = o) ; i.e.,  a sink  at  the  ground  if  (3.20b) 

d - 1 ^ 


K = 


K 


t 


Some  representative  results  for  the  burden  B0(t)  and  E^(t)  are 
shown  in  Table  B-2  for  various  values  of  the  parameters.  We 
see  that  for  short  times  (about  1 — 4 months)  the  a-deper.dence  of 
BT(t)  is  very  small,  but  for  the  times  of  interest  here,  which 
are  typically  8-l6  months,  there  is  a significant  variation  be- 
tween the  values  of  B,p(t)  for  d ~ 0.03  - 0.1  (which  are  considere 
to  be  representative),  and  the  limit  d ■*  0 . The  present  cal- 
culations are  reported  for  d = 0.03;  the  choice  of  d affects 
Bg(t)  and  thus  the  rate  of  transport  out  of  the  stratosphere 
for  a given  value  of  K,  but  the  numerical  value  of  the  eddy  dif- 
fusivity  in  the  troposphere,  = K/d,  is  not  important  for  the 
subsequent  discussion. 


TABLE  B-2.  STRATOSPHERIC  BURDEN  Bj(t)  FOR  THE  TWO-STEP  MODEL  AS  A FUNCTION 


OF  TIME,  FOR  VARIOUS  VALUES  OF  d = K 


/ K 

stratosphere  troposphere 


d Kstrat^Ktrop 

Time,  months 

1 

4 

8 

12 

16 

oa 

2.515 

1.993 

1.382 

1.017 

0.782 

0.01 

2.516 

2.003 

1.392 

1.026 

0.790 

0.03 

2.516 

2.020 

1.412 

1.043 

0.804 

0.1 

2.517 

2.066 

1.480 

1.105 

0.857 

0.3 

2.518 

2.131 

1.614 

1.258 

1.003 

lb 

2.520 

. 1 R . 

2.217 

1.807 

1.518 

1.303 

1 o 

We  show  values  of  B(t)/10  q for  injection  height  z^  = 18  km,  tropopause 
height  z0  = 11  km,  stratospheric  eddy  diffusivity  K = 10^  cm^/sec. 


aThis  is  actually  case  S,  z9  = z = 11  km  - see  Eq.  (B.20a) 

h c O 

This  is  case  S,  z2  = 0 - see  Eq.  (B.20b) 
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B . 3 . 2 Initial  Condition 


The  initial  Zr-95  profile  as  a function  of  height  used  here 
is  taken  from  Telegadas  (1974),  Fig.  4 5.  It  can  be  fitted  quite 
well  to  a Gaussian  profile 

f(z,tQ)  = fQ  exp  I - (z-z^  )2/cq2  J-  (B.21) 

with  a =2.15  km,  which  can  be  reconciled  with  the  delta-function 
o ' 

inj ection  of  Eq.  (B.5)  by  the  equivalence 


= 4 K t 


'calc 


t , + t 

obs  o 


(B. 22a) 
(B . 22b ) 


where  t , = time  from  a delta  function  injection,  and 

calc 

t , = time  after  the  nuclear  explosion, 

obs 

Representative  values  of  t are  shown  in  Table  B-3  as  functions 
of  K.  F r K i 3 x 10J  cmVsec  the  t(  values  are  so  small  com- 
j ir<  It  the  i Sf  rvat  Lor  1 Lines  shown  in  Fig.  B-l  that  the  equi- 
valence of  E'q.  (B.22)  is  indeed  satisfactory . 


TABLE  B-3.  TIME  CORRECTION  t AS  A FUNCTION  OF  K [SEE  EQ.  (B.22a).] 


K,  cm  /s 
■>3 


10“ 

3 x 10 
104 
3 x 104 


3 


tn,  months 

4.4 

1.5 
0.44 
0.15 


B . 3 . 3 Correction  for  Sedimentation 

Zr-95*  like  most  other  radioaet iv<  is  topes, 
small  particles  whose  mean  radius  lies  in  the  range 


carried  on 


0.02  ym  to 
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0.15  ym  (see  Drevlnsky  and  Pecci,  1965;  Telegadas  and  List, 
1969).  We  assume  that  the  density  of  the  particles  Is 
p = 2 g/cm  and  that  the  mean  radius  R = 0.1  ym,  as  did 
Telegadas  and  List  (1969).  Spherical  particles  of  this  density 

_ O 

and  radius  have  a sedimentation  speed  v = 3-6  x 10  J cm/sec 

s 

at  an  altitude  of  18  km  (see,  e.g.,  Junge  et  al.,  1961),  which 
corresponds  to  a net  settling  of  about  100  m per  month.  One 
way  to  correct  for  sedimentation  is  to  assume  that  the  mean  in- 
jection height  falls  with  the  effective  sedimentation  speed  v . 
Thus,  the  effective  injection  height  above  the  tropopause  is 


since  the  sedimentation  speed  depends  on  radius  and  on  ambient 
height  (or  density).  A plot  of  z^(t)  is  given  in  Fig.  B-3,  for 
p = 2 g/cm',  R = 0.1  ym. 

With  this  correction,  we  replace  xq  of  Eq.  (13)  by 

xc  = 1/2  (Kt)1/2  [(1/H)  - (zx  o(t)/Kt)]  , (B.24) 

and  thus  U(x  ) by  U(x  ). 

o o 

It  has  been  pointed  out  (by  H.S.  Johnston,  private  communi- 
cation, 1975;  also  Poppoff,  Whitten,  and  Turco,  1978)  that  the 
residence  time  of  aerosols  is  very  sensitive  to  particle  size 
and  could  be  affected  by  aerosol  growth  and  coagulation  processe 
In  keeping  with  the  use  of  a relatively  simple  model,  we  make 
the  present  simple  estimate  of  the  effects  of  sedimentation. 
Mote,  however,  that  the  effect  of  sedimentation  on  the  present 
estimate  of  stratospheric  K-values  is  particularly  significant 
for  relatively  small  stratospheric  diffusion  distances  (z^-z  ), 
i.e.,  for  high  values  of  the  assumed  mean  tropopause.  Overall, 
the  correction  for  sedimentation  introduces  significant  uncer- 
tainties, and  it  may  be  desirable  to  reexamine  the  present  esti- 
mate on  the  basis  of  any  additional  experimental  data. 


Z (0  , km 


20 

(”Z  (t  =0)  = 18  km;  p-  2 g/cm^;  R - 0.  1 pm 


FIGURE  B-3.  Effective  injection  height  for  aerosols  as  a 
function  of  time,  7(t). 


B.3.4  Data  Reduction  for  (K,zq)  Parameterization 

To  reduce  the  data,  we  normalize  the  tracer  concentrations 
shown  in  Pig.  B-l  to  give 


A _ tracer  burden  (corrected  for  decay 

exp'  ' ~ maximum  burden  observed  for  the  event 


(B. 25) 


This  quantity  is  shown  in  Fig.  B-4,  where  we  also  show  calcu- 
lated values  of  the  normalized  tracer  burden,  A , (t),  for 

C3.1C 

the  following  set  of  parameters: 

zq  = 11  km  - height  of  tropopause 
z^  = 18  km  - height  of  injection 

from  Fig.  B-3  to  correct  for  sedimentation  from 
an  initial  18-km  injection. 


We  see  from  Fig.  B-4  that  experimental  data  for  all  five 

events  extend  to  A s 0.6  but  not  to  A ^ 0 . 5 , and  thus  we 

exp  exp  ’ 

define  the  quantity  tg  g by  the  relation 


Aexp(t0.6)  ’ °'6 


(B.26) 


and  use  tQ  g rather  than  t-^^,  which  cannot  be  determined  from 
all  five  events)  to  compute  K for  a given  value  of  zq. 

Note  in  Fig.  B-4  (or  Fig.  B-l)  that  the  stratospheric  bur- 
den due  to  the  October  1970  shot  falls  rather  abruptly  some  7 
months  after  the  event.  While  this  decrease  may  be  spurious 
(since  the  later  values  of  burden  are  lower-bound  values  due 
to  contamination  from  the  French  tests  of  summer  1971),  yet  it 
may  be  correlated  with  the  fall  in  mean  altitude  in  the  profile 
some  4 months  after  the  event  shown  in  Fig.  B-2.  In  the  analysis, 
we  shall  use  the  three  events  of  June  1967  (slow  decay-summer), 
September  1969  (intermediate  decay-fall),  and  December  1968 
(fast  decay-winter)  as  representative  of  the  data.  The  charac- 
teristic times  tn  c for  these  events  are  listed  in  Table  B-4. 


FIGURE  B-4.  Decay  of  stratospheric  tracer  burden  as  a function  of  time. 

Comparison  of  results  from  the  Chinese  tests  (from  Telegadas, 
1974,  1976)  with  the  predictions  of  the  "T"  model  with  and 
without  sedimentation  for  various  eddy  diffusivity  (K)  co- 
efficients. Note  that  the  final  data  points  marked  with  an 
arrow  are  lower-bound  values,  as  the  radioactivity  data  from 
the  Chinese  tests  were  obscured  by  new  fallout  from  later 
(French)  tests.  Tropopause  height  taken  as  11  km,  d = 0.03. 
Note  that  the  data  of  10/70  fall  abruptly  at  *;  this  corre- 
lates with  the  earlier  motion  down  of  the  peak.  Also,  the 
data  of  12/68  have  been  extrapolated  to  the  12-month  point  P 
used  to  determine  (12  mos.)  for  the  analysis  of  Eq.  (27) 


TABLE  B-4.  EXPERIMENTAL  PARAMETERS  FOR  VARIOUS  TESTS 
USED  IN  THE  DATA  ANALYSIS 

Test  of 

Parameter  12/68  10/70  9/65  6/73  6/67 

tQ  g (month)  5.7  8.2  9.0  9.45  11.7 

A0xp  (12  months)  (0.305)  (0.555)  0.59 

In  Fig.  B-5,  we  show  K as  a function  of  t„,  c for  z =11, 

611  U . O O 

14  km,  Zj  = 18  km,  and  as  z,  (t)  from  Fig.  B-3,  correcting  for 
sedimentation  as  is  suggested  in  Section  B.3-3- 

The  qualitative  character  of  the  results  is  plausible: 

a.  Increasing  z decreases  (z.-z  ) and  thus  [cf.  Eq . (B.17)’ 

O 1 o 

decreases  K for  a given  time  delay;* 

b.  Sedimentation  accelerates  transport;  i.e.,  a smaller 

value  of  K is  needed  if  sedimentation  is  considered 

than  if  it  is  neglected.  In  particular,  for  a 14- km. 

tropopause,  which  only  requires  transit  through 

(z,-z  ) = 4 km  in  the  stratosphere,  very  little  diffusion 
1 o 

is  needed  to  supplement  the  sedimentation  as  predicted 
by  the  present  model. 

c.  The  experimental  results  in  summertime  call  for  a much 
smaller  value  of  K than  do  the  winter  results. 

A different  way  to  analyze  the  data  is  to  obtain  an  annual 
average  of  activity;  i.e.,  to  ask  for  the  activity  12  months 
after  the  event  and  to  match  this  with  the  model  as  a function 
of  K,  rather  than  by  asking  for  the  time  tg  g of  Eq.  (B.26). 

We  see  from  Figs.  3-1  and  B-4  that  for  the  fast  event  of 
December  1968,  there  are  no  data  as  late  as  12  months  after 
If 

Note  that  the  quantitative  relation  of  Eq . (B.17)  that 
for  a given  removal  rate,  the  ratio  of  K to  (zi-zQ)  is 
constant,  is  not  satisfied  for  the  present  model  T. 


( 0.6  l,  month* 


FIGURE  B-5.  Time  for  stratospheric  tracer  burden  tc  decrease  to 
60  percent  of  Its  peak  value,  tg  g,  computed  with  and 

without  sedimentation,  shown  as  a function  of  eddy 
diffusivlty  K.  Calculations  are  shown  for  tropopause 
height  Zq  x 11  and  14  km.  The  corresponding  decay 
times  calculated  from  the  various  nuclear  explosions 
are  also  shown. 


the  event,  but  by  extrapolation  we  obtain  the  point  P in  Pig.  B-4 
which  is  used  here;  the  values  used  here  are  listed  in 

Table  B-4. 

Explicitly,  we  obtain  K from  the  relation 

= 12  months)  = A , [t  , =12  months  + t (K);  (K)]  (B.27) 

exp  obs  calc  calc  o 

In  Table  B-5  we  show  effective  mean  K-values  in  the  lower 

stratosphere  obtained  from  t_  c [see  Eq . (B.26)]  and  from  A 

0.6  M exp 

(12  months)  [see  Eq.  (B.27)]. 

The  following  points  should  be  noted: 

a.  The  K-values  corrected  for  sedimentation  are  considered 
preferable  to  those  not  so  corrected  (but  note  the 
limitations  of  the  correction,  especially  for  zq  = 14  km). 

b.  The  two  different  procedures  of  determining  K,  of 

Eq . (B.2 6)  and  Eq.  (B.27),  give  rather  similar  numeri- 
cal values;  the  differences  can  be  understood  on  phy- 
sical grounds. 

c.  The  significant  difference  between  K-values  for  summer 
and  winter  injection  can  readily  be  understood  in  terms 
of  the  experimental  results  of  Pig.  B-l. 

c.  Calculations  have  also  been  made  for  varying  injection 
heights,  = 16,  20  km,  with  reasonable  consequences, 
in  that  larger  values  of  require  smaller  values 

of  K to  compensate,  etc. 

B . 4 APPLICATION  TO  STEADY-STATE  INJECTION 

Given  the  effective  stratospheric  K-values  deduced  here, 
it  is  possible  to  compute  atmospheric  residence  times  and  in- 
jection coefficients  for  the  continuous  injection  of  tracers 
(injected  at  18  km  altitude  and  40°N  latitude),  by  using  the 
present  model  for  a K-profile  [cf.  Eq . (B.2)]. 
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TABLE  B-5.  EFFECTIVE  K-VALUES  COMPUTED  FROM  EXPERIMENTAL  DATA 

(IN  10^  cm^/sec) 


Experimental  Situation 


June  1967 
Slow-Summer 


September  1969 
Intermediate-Fall 


December  1968 
Fast-Winter 


Tropopause  at  11  km 

From  A(12  mo. ) , i .e. , 
12-month  average 


With  sedimentation 

4.0 

4.5 

10.5 

Without  sedimentation 

5.5 

6.1 

14.0 

From  tQ  6 - decreases  to 

0.6  of  maximum 

With  sedimentation 

4.3 

6.2 

10.3 

Without  sedimentation 

6.1 

7.9 

12.0 

Tropopause  at  14  km 

From  A(  12  mo. ) 

With  sedimentation 

0.77 

0.95 

4.1 

Without  sedimentation 

2.2 

2.6 

6.9 

From  t0>6 

With  sedimentation 

0.5 

1.35 

3.0 

Without  sedimentation 

1.8 

2.9 

5.2 

NOTE:  We  use  model  "T",  with  d = 

^strat/Ktrop. 

= 0.03 
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The  atmospheric  residence  time  tR  is  defined  as  atmospheric 
burden/flux,  while  the  injection  coefficient  a is  defined  for 
a steady  state  by  McElroy  et  al.  (1974) — see  also  Hunten  (1975)-- 
as  follows: 

a = X/Q  (B.280 


where  X - molecular  mixing  ratio  at  and  above  the  injection 
height  z-L  (which  is  constant  in  a steady  state),  and  Q = source 
strength . 


For  the  case  of  an  isothermal  stratosphere  with  scale  height 
H (=  6.3  km,  corresponding  to  T = 2l6°K),  and  an  isothermal 
troposphere  with  scale  height  H , and  for  the  K-profile  of 
Eq.  (B.2)  we  have  the  following  expressions  for  gaseous  tracers: 


(B.29) 


(B. 30) 


For  aerosols,  the  treatment  is  more  complex  [see  Bauer, 
Gardner,  Oliver  (1977)]-  Here  we  confine  the  discussion  to 
gaseous  injections,  which  are  of  greatest  interest. 

Given  these  formulae,  we  can  compute  injection  coefficients 
and  residence  times  for  various  values  of  the  spatial  and  tem- 
poral mean  eddy  diffusivity  coefficient  K (between  18  km  and 
the  tropopause)  as  derived  from  the  experimental  data.  Before 
doing  so,  however,  some  discussion  of  the  appropriate  mean  on  a 
time-weighted  basis  is  necessary. 

Let  us  reexamine  how  material  leaves  the  stratosphere.  From 
Fig.  B-l,  it  appears  that  material  injected  in  summer  does  not 
leave  at  all  until  the  following  winter.  However,  the 
month)  estimates  of  Eq . (B.27)  do  provide  annual  mean  K-value 
for  a summer  injection.  Here  we  shall  use  an  annual  average 


K 


summer 


(B.  3D 


= 1/2  (K 


+ K 


winter 


) 


based  on  the  A (12-month)  estimates,  taking  into  account  the 
effects  of  sedimentation.  The  results  obtained  in  this  way  are 
compared  with  those  of  some  other  studies  in  Table  B-6.  Refer- 
ring to  Fig.  B-2,  we  see  that  an  annual  mean  tropopause  height 
zq  = 11  km  is  most  appropriate;  we  show  results  also  for  a high 
tropopause,  z0  = lA  km,  for  which  the  short  diffusion  distance 
implies  that  K is  quite  small  and  the  effects  of  sedimentation 
are  extremely  important. 

In  Fig.  B-6  we  compare  the  present  estimates  of  K with  these 
of  several  current  models  listed  in  Table  B-6.  Note  the  effects 
of  sedimentation  on  an  annual-average  basis,  and  in  particular 
the  wide  variation  between  effective  K-values  for  summer  and 
winter  injection. 


B.5.  A FINAL  COMMENT 

The  most  striking  characteristic  of  the  radioactive  fallout 
data  of  Telegadas  (197A  and  1976)  is  that  they  display  a very 
strong  seasonal  variation,  with  transfer  to  the  troposphere 
being  much  more  rapid  in  late  fall  and  winter  than  in  summer. 
This  effect  is  not  included  in  most  one-dimensional  models,  in 
which  an  annual  average  eddy  diffusivity  profile  is  used,  and 
may  have  resulted  in  an  overestimate  by  such  models  of  the  cal- 
culated ozone  reduction  following  the  massive  Soviet  atmospheric 
nuclear  tests  of  fall  1961  and  fall  196?  [see  Foley  and 
Ruderman  (1973),  Bauer  and  Gilmore  (1975),  and  p.  159  of  MAS 
(1975)]. 
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TABLE  B-6.  COMPARISON  OF  RESULTS  FROM  DIFFERENT  MODELS 
FOR  18-km  INJECTION 


Source 

KT 

a, 

cm2-sec 

tR, 

ZL_ 

This  analysis  (two-step  case  "T",  d = 

0.03) 

11  km  tropopause 

With  sedimentation  ("K  = 7.25  x 103 

0 

cm  /sec) 

2. 

4 

2.0 

Without  sedimentation  (l<  = 9.75  x 

103  cm2/sec) 

1. 

8 

1.5 

14  km  tropopause 

With  sedimentation  ("K  = 2.44  x 103 

o 

cm  /sec) 

5. 

0 

3.5 

Without  sedimentation  fK  = 4.55  x 

103  cm2/sec) 

2. 

5 

1.9 

Other  K-Profiles* 

‘ 

Chang  (1974) 

1. 

7 

1.4 

Chang  (1976) 

1. 

8 

1.6 

Hunten  (1974) 

4. 

6 

3.6 

Hunten  (1974)  (+2) ( 18  km. 

7. 

1 

4.6 

latitude  adjusted  to  20  km) 

• , 

Crutzen  (1974) 

1. 

8 

1.7 

Crutzen-Isaksen  (1975) 

1. 

9 

1.7 

Wofsy  (1975) 

3. 

2 

2.6 

Crutzen  et  al . (1978) 

2. 

3 

2.0 

★ 

The  more  recent  (1975-78)  eddy  diffusivity  profiles 

used 

here  are 

shown 

in  Fig.  B-6;  the  older  ones,  which  were  used  during  CIAP,  are  listed  in 
Oliver  et  al.  (1977).  All  estimates  are  made  for  18-km  injection,  except 
for  the  "Hunten  (1974)  (+2)"  estimate  which  corresponds  to  a 20-km  injec- 
tion as  an  adjustment  for  the  latitude  variation  of  mean  tropopause  height, 
the  procedure  recommended  in  NAS  (1975),  p.  118.  Conceivably,  other  pro- 
files should  be  similarly  adjusted,  but  we  have  seen  no  such  recommendation 
for  the  other  profiles. 
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FIGURE  B-6.  The  present  results  for  an  11  km  tropopause  compared  with  some 
other,  currently  used,  eddy  diffusivity  profiles.  Note  the 
effect  of  sedimentation  (see  Section  B.3.3)  on  an  annual  av- 
erage bases,  and  also  the  large  summer/winter  variation. 
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